Primordial Black Holes in 
Non-linear perturbation theory 



Juan Carlos Hidalgo Cuellar 

Astronomy Unit 
School of Mathematical Sciences 
Queen Mary College 
University Of London 
United Kingdom 



A dissertation submitted in candidature for the degree of Doctor of Philosophy 

in the University of London 



Friday 1st May 2009 



Abstract 



The thesis begins with a study of the origin of non-linear cosmological fluctuations. In partic- 
ular, a class of models of multiple field inflation are considered, with specific reference to those 
cases in which the non-Gaussian correlation functions are large. The analysis shows that pertur- 
bations from an almost massless auxiliary field generically produce large values of the non-linear 
parameter /nl- 

Next, the effects of including non-Gaussian correlation functions in the statistics of cosmo- 
logical structure are explored. For this purpose, a non-Gaussian probability distribution function 
(PDF) for the curvature perturbation 7Z is required. Such a PDF is derived from first principles in 
the context of quantum field theory, with n-point correlation functions as the only input. Under 
reasonable power-spectrum conditions, an explicit expression for the PDF is presented, with cor- 
rections to the Gaussian distribution from the three-point correlation function (1Z1Z1Z) . 

The method developed for the derivation of the non-Gaussian PDF is then used to explore 
two important problems in the physics of primordial black holes (PBHs). First, the non-Gaussian 
probability is used to compute corrections to the number of PBHs generated from the primordial 
curvature fluctuations. Particular characteristics of such corrections are explored for a variety of 
inflationary models. The non-Gaussian corrections explored consist exclusively of non-vanishing 
three-point correlation functions. 

The second application concerns new cosmological observables. The formation of PBHs is 
known to depend on two main physical characteristics: the strength of the gravitational field 
produced by the initial curvature inhomogeneity and the pressure gradient at the edge of the 
curvature configuration. The latter has so far been ignored in the estimation of the probability 
of PBH formation. We account for this by using two parameters to describe the profile: The 
amplitude of the inhomogeneity and its second radial derivative, both evaluated at the centre of 
the configuration. The method developed to derive the non-Gaussian PDF is modified to find the 
joint probability of these two parameters. We discuss the implications of the derived probability 
for the fraction of mass in the universe in the form of PBHs. 
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Chapter 1 

Introduction 



Cosmology is at the forefront of modern physics. Over the last two decades, it has moved 
from a predominantly theoretical discipline to a sound observational science. Today's ex- 
periments are capable of observing tiny fluctuations of a faint signal coming from the Big 
Bang, emitted about thirteen billion years ago. The observations of primordial inhomo- 
geneities are a unique probe of the physical conditions in the early universe. The inflationary 
paradigm indicates that the inhomogeneities are the result of quantum fluctuations of the 
matter dominating the universe in its first moments. In this widely accepted picture the ob- 
served inhomogeneities fix the normalisation of an inflationary potential setting the energy 
scales for inflation to about 10 16 GeV, the GUT scale. This is 10 6 times more than the energy 
of particles released by supernovae. A similar ratio arises for the energy scales to be tested 
by the large hadron collider (LHC). These numbers show how the geometry of the universe 
and its inhomogeneities constitute a unique probe of high energy physics. 

Several observational parameters have been defined in cosmology in order to determine 
the physical conditions of the early universe. The density and nature of matter observed 
today, the distribution and mean amplitude of initial inhomogeneities, and most recently 
the non-Gaussianity of primordial fluctuations are among these parameters. The latter has 
received considerable attention from cosmologists but the analysis of the latest observations 
has not yet provided conclusive evidence for departures from Gaussian statistics. A great deal 
of effort is under way to reduce the detection thresholds of the non-Gaussian parameters. 
Even if non-Gaussianity remains undetected by future experiments, we can still constrain 
theoretical models that are known to develop large non-Gaussianity. 

The main objective of the present work is to study how non-Gaussian statistics, inherited 
from inflation, can modify the probability of primordial black hole formation. The class of 
models of inflation that motivate this study and the development of statistical tools to address 
this question are complementary projects, and both are included in the present thesis. In the 
rest of this chapter we provide a brief description of the state of the art in cosmology, with 
special attention to the open questions that motivate this thesis. 
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1.1 Cosmological observations and the Big Bang 



It has been more than four decades since 



Penzias & Wilson! ll 196511 managed to identify, 



for the first time, the cosmic microwave background (CMB) radiation. This was detected, 
almost by accident, while calibrating a large reflector at the Bell Laboratories. The uniform 
and isotropic radiation observed corresponds to the most perfect black-body radiation ever 
measured, peakin g at A = 1.9 mm, with a red- shifted temperature of T CMB = 2.725 Kelvin 
jjaffeetallEoOlll . 

The detection of the CMB gave decisive support to the Big Bang theory. The standard Big 
Bang model considers a universe dominated by uniform and isotropic matter. Its dynamics 
is governed by gravity, with equations prescribed by the theory of general relativity. (Grav- 
ity is the only long-range force to be considered since the universe is electrically neutral.) 
The conditions of isotropy and homogeneity, in this context, imply that the spacetime admit- 
ting these pr operties is necessarily a Friedmann-Robertson- Walker universe (FRW) (see e.g. 
Waldl il984h . 



1.1.1 Basic dynamics of the universe 

We write the FRW metric in the form of the line-element in spherical coordinates. 



ds 2 



-dt 2 + 



a(t) 
1 — nr' 



(dr 2 + r 2 [d9 2 + sin 2 6d<f?]) , 



(1.1) 



where t and k are the coordinate time and the uniform curvature of the spatial sections 
respectively. The usual spherical coordinates in the spatial hypersurfaces are r, 9 and <\>. 
Finally, a(t) is the scale factor, with present value a = 1. The Einstein equations of gen- 
eral relativity provide the dynamical relation between the matter and spacetime variables. 
Assuming homogeneous and isotropic matter, with density p and isotropic pressure p, the 
Einstein equations show that the evolution of the scale factor is given by 



H 2 



K 



(1.2) 



where an over-dot is the coordinate time derivative and H is the Hubble parameter, a measure 
of the expansion rate. I ts present value is H = lOO/i kms _1 Mpc _1 , with h = 0.71 ± 0.08 
IIFreedman et all 1200 ill . This last equation is known as the Friedmann equation. We use 
throughout units where c = h = 8irG = 1. 

The matter contents of our universe has several components and the fraction of each 
component relative to the critical density is called the density parameter fij = /3H 2 . If 
we denote the sum of all matter components as fix, the Friedmann equation can be written 
simply as 



fi K (t) + n T (t) = i, 



(1.3) 
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where Q K = K/(aH) 2 is the curvature density parameter. When the matter density is equal 
to the critical density 3H 2 , then Ox = 1 and the universe is flat at all times. Observations 
tell us that we live in a nearly flat universe (\Q K \ < 1CT 2 ), so we assume Q K = hereafter. 
The energy density is dominated by two main components, a cold dark matter component 
(^cdm — 0.23) and another component referred as dark energy ~ 0.72). The nature of 
both these components is a crucial question in cosmology and has motivated a lot of research. 
We will return to this point and to an analysis of the Einstein equations later in this work. 



1.1.2 The Big Bang model 



The hot Big Bang model is now accepted as the standard model describing the evolution of 
the universe. This model characterises, with impressive accuracy, the evolution after the first 
second. At this time, the universe was a primordial fireball with high enough temperature 
and pressure to dissociate any nuclei. The formation of nuclei was only possible once the 
cosmic expansion reduced the average kinetic energy sufficiently. The formation of the first 
elements took place at temperatures of around T ~ 0.1 MeV, when the universe was around 
1 s old. This process involves conditions that cannot be replicated elsewhere (cf. stellar nu- 
cleosynthesis). Within the current observational limitations, the Big Bang prediction for the 
present abundance of light elements is confirmed remarkably by the present measurements. 

Big Bang nucleosynthesis halted once matter had cooled down enough, due to the cos- 
mic expansion. The electrical neutrality of the matter was reached at a more recent event: 
the so called 'recombination' process refers to the time when each electron was captured 
by a nucleus forming the first neutral atoms. Subsequently, at a temperature of around 
T ps 0.1 eV (ps 10 3 Kelvin), CMB photons decoupled from ordinary matter and have since 
travelled freely. These same photons reach us in the form of microwave radiation. The 
surface of emission of these primordial photons is called the last- scattering surface. CMB 
observations constitute irrefutable proof that the universe was homogeneous at early epochs 
and dominated by radiation when T > 10 3 Kelvin. 

The current temperature of the CMB radiation (Tcmb = 2.725 Kelvin) is measured with 
such precision because its fluctuations are ti ny. The first observational ev i dence for the CMB 



1992 



Bennett et all 119961. The results 



anisotropics came from the COBE satellite [Smoot, 

of this experiment showed that the temperature fluctuati ons have a mean amp litu de 8T/T 
■0~ 5 . The amplitude of such deviations was predicted by Peebles & Yu [1970] and lzeldovich 
lll972ll in terms of the matter density pertu rbation 8p/p ~ 10~ 5 . T hese inhomogeneities are 



related through the Sachs-Wolfe formula [[Sachs & WolfeL LL967Q . This prescribes that for 
inhomogeneities of comoving size A, 



6T ~ l ( IT U 2 A 



(1.4) 



where we have defined 5 P = Sp/p, and where a subscript LS indicates an evaluation at the 
last-scattering surface. 
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More recent expe riments, such as BOOMERANG INetterfield et all 120 021. MAXIMA 



IIHanany et all 120001 and WMAP IIHinshaw et al. 



2007 



Komatsu et al. 



2008], managed to 



measure the acoustic oscillations in the radiation plasma due to the small-scale density varia- 
tions in the early universe. Measurements of acoustic oscillations in the CMB demonstrated 
the flatness of the universe to 1% precision (i.e. \Q K \ < 10~ 2 ). They were also used to rule 
out cosmic strings as a significant contributor to struc ture formation an d suggested 'cosmo- 
logical inflation' as the theory of structure formation liJaffe et all 1200 ill . 



1.2 Cosmological inflation 
1.2.1 Motivation and achievements 

The observations mentioned above provided strong arguments in favor of the Big Bang 
model but also showed the necessity of a larger theoretical framework due to the follow- 
ing problems: 

1 . Horizon problem. In the Big Bang model, the distance light could have traveled up to 
the time of last- scattering c/ls is of order 180 Mpc. This is called the particle horizon 
and determines the radius of causally connected regions at that time. The particle hori- 
zon today is much larger, with radius d ~ 6000 Mpc. Therefore, the measurements 
of CMB radiation at angular scales larger than one degree include regions that were 
causally disconnected at the time of the photon decoupling. The temperature at such 
scales is observed to be uniform up to one part in 10 5 . This means that causally disjoint 
patches of the universe in the past had the same thermal history. In the context of the 
hot Big Bang model there is no plausible explanation for this fact. 

2. Flatness problem. The density of matter components in the universe is diluted with 
time due to the cosmic expansion. Conversely, if there was an initial curvature com- 
ponent k, then this would rapidly dominate the matter contents. This is easily derived 
from Eq. (11.31) . which can be written in the form 

^-i = ^— = n n . (1.5) 

The product aH decreases with time in a radiation or matter dominated universe. If 
the universe is initially flat, then it remains flat for subsequent times, but observations 
show that \Q K \ < 10~ 2 today, and the Friedmann evolution demands an even smaller 
curvature in the past. For example, at nucleosynthesis, when the universe was around 
1 s old, we require \Q K \ < 10~ 16 to be consistent with the present value. Such a 
small value requires an extreme fine-tuning of initial conditions f2 T , for which a causal 
explanation would be desirable. 

A solution to these problems is provided by the inflationary paradigm, which we will 
study in detail in Chapter [2] The main feature of this theory is that it changes the behaviour 
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of the comoving cosmological horizon by considering an accelerated expansion of the uni- 
verse at early times, i.e., at times prior to nucleosynthesis. In terms of the scale factor, this 
condition demands 



a > 



dt 



1 



aH 



< 0. 



(1.6) 



The shrinking of the cosmological horizon represents a 'reverse' evolution of spacetime 
which avoids the fine-tuning of initial conditions demanding homogeneity and flatness. If 
we consider an inhomogeneous patch of the universe when inflation starts, at an initial time 
U, the cosmic accelerated expansion brings all initial inhomogeneities out of the comoving 
cosmological horizon. If inflation lasts long enough, then after the inflationary period we are 
left with a much larger region composed of small patches of size of the cosmological horizon 
which are out of causal contact but with common physical characteristics. The number of 
e-folds of expansion required for the listed problems to be solved is 



N = In 



a(t 



end i 



a(ti) 



> 60. 



(1.7) 



This number is required to guarantee that the comovi ng scale of the current size of the uni- 
verse exited the horizon at the beginning of inflation IILiddle & Lv th. 2000]. This indicates 
that inflation must last longer than 60 e-folds. Arguably, it was Guthl 1 198 lh who first brought 
these ideas together. 

The theory of inflation has received important contributions from particle physics. In par- 
ticula r, the theory of particle creation from vacuum fluctuations [|HawkingLll982l : IStarobinskyL 
198211 gave inflation its strongest argument: the vacuum fluctuations generated during infla- 
tion are redshifted to superhorizon scales by the action of the inflationary mechanism. At the 
end of inflation, the thermalisation of the inflaton false vacuum reheats the universe and the 
standard hot Big Bang phase begins. In this transition, the vacuum fluctuations of the infla- 
ton field are transformed into matter density perturbations with a prescribed amplitude. From 
this transition onwards, the modes re-enter the expanding comoving horizon. Thus, initial 
conditions of cosmological perturbations in the hot Big Ba ng are set by inflation. The ob 



served mean amplitude of the temperature inhomogeneities [|Smooti 1 1992c iNetterfield et al 



2002 



Spergel et all 1200711 sets the energy scale at which the initial vacuum fluctuations were 



generated by tracing back the evolution of fluctuations described above. This simple expla- 
nation of the origin of the temperature fluctuations constitutes a decisive argument in favour 
of the inflationary scenario. It represents the greatest advantage of inflation over many other 
alternative extensions of the standard Big Bang scenario. 

In summary, the requirements for a period of inflation are: (1) a mechanism to generate 
an accelerated expansion maintained for at least 60 e-folds of expansion; (2) a way of ac- 
counting for the transition to the subsequent FRW stages of evolution, thereby providing the 
suitable initial conditions for the Big Bang scenario; (3) quantum fluctuations of the infla- 
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tionary field, generated at observable scales such that the matter density fluctuations of size 
A meet the relation (aHX)5 p ~ 1CT 5 and this product is almost invariant over the observed 
scales. 

In practice, measurements of CMB anisotropics, combined wi t h measurements of back - 



ground parameters inferred from supernovae surveys IIAstier et all 120061 : 



Riess et al 



2007|], 



indicate that the root-mean- square (RMS) amplitude of temperature fluctuations is 



5T 
T 



2 x 10" 5 , 



(1.8) 



RMS 



at the pivot scale with comoving size Acmb = 150 Mpc customarily used in CMB studies. 
Observations also indicate that this value does not vary significantly over the range of ob- 
served scales. In other words the mean amplitude is almost scale-invariant for angular scales 
larger than one degree. In Chapter [2] we show how this relates to the curvature perturbation £ 
and discuss its basic properties. In particular, we will show that, in the cases which concern 
us, £ is constant for scales larger than the particle horizon. 



1.2.2 An embarrassment of richness 

The required amount of inflation and the corresponding amplitude of the curvature pertur- 
bations determine the kind of matter and energy scale necessary to satisfy the conditions 
for accelerated expansion. These prerequisites have been met by several models of inflation 
which may or may not be motivated by more fundamental theories of physics. One of the 
main problems faced by the inflationary paradigm is that of richness. There are many mod- 
els that meet the dynamical requirements. Most of them invoke one or more scalar fields 
{(pi} with dynamics governed by a potential V(0j). There are a plethora of models, each of 
which corresponds to particular realisation of this potential, which satisfy the observational 
constraints up to the level of the observed inhomogeneities. Consequently, many of the 
models cannot be distinguished at the level of linear perturbation theory. This demands the 
formulation and experimental determination of new parameters that provide complementary 
information about the early universe. An important constraint on the inflationary models can 
be obtained by considering the statistical deviations from a Gaussian field of fluctuations. 
This idea has opened a new window in the study of the early universe, namely the nonlinear 
extension of perturbation theory and its non-Gaussian statistics. 



1.3 Non-Gaussianity 

By non-Gaussianity in cosmology we refer to the small deviations of observed fluctuations 
from the random field of linear, Gaussian, curvature perturbations Cj(i, x). C(^? x ) is me 
curvature perturbation in the comoving gauge, that is, as measured by an observer which 
sees no net-momentum flux. The mathematical expression for x) in terms of the matter 
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density perturbation is provided in Chapter [2l 

Among the parameters of nonlinearity, the nonlinear coupling /nl is the most useful 
observable for describing non-Gaussianity. Its definition comes from the second order ex- 
pansion of curvature perturbations in real space, which can be written as 



C(x) = Cl + \& 

where d refers to the Gaussian perturbation with variance £ 



x 



(1.9) 



tmss( x ) and C2 is the 



second order perturbation parametrised by the nonlinear parameter / NL in the following way 



C 2 (x) = -^nl(Ci(x) 2 -Crms(x) 2 ). 
Note that the perturbative expansion of £ implies also the rough definition 



(1.10) 



5C 2 (x) 



NL 



6Ci 2 (x) 



(1.11) 



which gives an intuitive notion of this parameter. Histori cally, non-Gaussia nity as a test of 
the accuracy of perturbation theory was first suggested by 
of /nl used here was first introduced by lSalopek & Bond 



Alle n et all I1987h . The definition 



ll 199011 in terms of the Newtonian 



or Bardeen pot ential <£r (defined in Chapter [2b . Their initial definition ha s been preserved 



by convention llGangui et all 1 19941 : 



Verde et al. 



2000 



Komatsu & Spergell . 1200 ill , which is 



why the transformation to the curvature perturbation ( 2 involves the numerical factor —5/6. 
In the context of perturbation theory, the study of dynamical equations at second order yields 
important information independent of the parameters of linear perturbations. Thus, in the 
nonlinear regime, we can discriminate different models of inflation which are degenerate at 
linear order. This fact has motivated the search for non-Gaussianity in the CMB and large- 
scale structure. 

Statistically, the lowest order effect of including a non-Gaussian contribution is a non- 
vanishing correlator of three copies of the curvature field (. The three-point function in 
Fourier space is given by the bispectrum B, defined by 



(C(k 1 )C(k 2 )C(k 3 )) = (27T) 3 B c (k 1 ,k 2 ,k 3 )5^(k 1 + k 2 + k 



3), 



(1.12) 



where 5^ is the three-dimensional Dirac delta function. 

The bispectrum is directly related to the parameter /nl and for each mode k — |k| . More- 
over, being a function of three momenta, the ^-dependence of the bispectrum also provides 
valuable information which could help us to understand the physics of the early universe. 

The nonlinear parameters have been investigated through the analysis of hi gher order cor- 
relati ons in the CMB anisotropics observed mostly by the WMAP satellite jSpergel et al. . 
200711 . After five years of collectin g data, WMAP observations give the limits —151 < 
/nl" 1 < 253 I Komatsu et al. . I2OO8 I for an equilateral triangulation of the momenta and 
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-4 < /j|^ al < 80 ISmithet all 1200911 for a local triangulation. The triangulation of the bis- 
pectrum is a characteristic which arises due to the following: The momentum conservation in 
the three point correlation is guaranteed by the delta function in Eq. (11.121) . which demands 
that the sum of the three vectors is zero. In consequence the three momenta represent the 
sides of a triangle in A;-space. Two main triangulations can be distinguished: the equilat- 
eral triangulation and the isosceles or l ocal triangulation, w hich are characteristic shapes of 
different models of inflation (see e.g. iBabich et al.l [2004J]). An experimental detection of 
/nl would greatly narrow the range of cosmological models which meet the observational 
bounds. In the near future, space telescopes, and in particular the PLANCK satellite, are ex- 
pected to tighten these bou nds considerably. Specifica 



•P 
1, 



ly, any signal with J fm.\ <; 5 should 
be observed by PLANCK IlKomatsu & SpergelL 120011 : iLiguori et all boodl . This raises the 
exciting possibility of looking for particular signatures of inflationary models. 

Another attractive observational prospect for non-Gaussianity is to look at the implica- 
tions of considering primordial non-Gaussian fl uctuations in the s t udy of the statistics of 
galaxies and other la rge-scale structures (LSS) llVerde et all \200(x iMatarrese et all 



200C 



Lo Verde et all l2008|1 . Such observations probe inhomogeneities at scales smaller than those 



observed in the CMB. 

The effects of non-Gaussianity in the LSS can be classified into two categories, which 
provide distinct observational methods for detecting non-Gaussianity. The first is the bis- 
pectrum of galaxies, potentiall y determined by computing the three-point correlation func- 



tion from redshift catalogues llVerde et all 12001 ; 



Scoccimarro et al 



200411 . The second is 



the non-Gaussian correlations in the probability distribution function (PDF) which leads to 
modi fications in t he number of galaxies and o ther structures with respect to the Gaussian 
case I Verde et aL . 



20001 : IMatarrese et all [2000] . 



Both methods involve delicate issues, crucial for the correct interpretation of observa- 
tions. Most important is the fact that the inhomogeneities that collapse to form galaxies 
evolve in a nonlinear fashion at late times. This is because the primordial fluctuations en- 
ter the horizon much before they form virialised structures. Consequently, the nonlinear 
evolution of fluctuations may blur the primordial non-Gaussianity of the initial statistics. 

Another important problem is that there is no single way of constructing a non-Gaussian 
PDF from theoretical models, i.e., several non-Gaussian PDFs can be constructed with a 
common varian ce and skewness. This well known problem has been expressed pithily by 
n2006n : "We know what a dog is, but, what is a no-dog? A no-dog can be anything". 



Heavens 



The effects on, say, the integrated number of galaxies may change substantially with every 
realisation of the PDF. This complicates the interpretation of non-Gaussian signatures. 

In Chapter [3j a formalism is presented to attack this problem. We construct the PDF 
of the curvature perturbations with a direct input from its higher-order correlations. The 
formalism is then applied to compute the modification which a non-Gaussian distribution of 
fluctuations brings to the abundance of primordial black holes. 
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1.4 Primordial black holes 



1.4.1 Standard picture 



The idea that large amplitude matter overdensities in the universe could have collapsed 



throug h self-gravity to form primo rdial black holes ( PBH s) was first put forward b y 
[| 1966TI and then independently by iHawkingl 11197 ill and ICarr & Hawking! II 197411 more than 
three decades ago. They suggested that at early times large- amplitude overdensities would 
overcome internal pressure forces and collapse to form black holes. The standard picture of 
PBH formation from initial inhomogeneities prescribes that an overdense region with size 
r\ will overcome pressure and collapse to form a black hole if its size is bigger than the 
associated Jeans length 



Zel'Dovich & Novik 



rj = 4.7T 



■w 



5 + 9w 



■d 



Hi 



(1.13) 



where the particle horizon dn is of order of the Hubble radius th = 1/H. Here we assume 
an equation of statep = wp, where w is constant. For the case of radiation-domination, for 
example, w = 1/3. 

The size of the initial inhomogeneity must also be smaller than the separate universe scale 



H 



f{w), 



(1.14) 



where the function f(w) has been derived by lHarada & Card [1200511 . and is of order unity. 
Thus, rj < r- x < ru, both limits being of order the Hubble radius. Consequently the mass of 
a PBH is close to the Hubble horizon mass. This g ives a simp le formula for the mass of a 
PBH forming at time t during radiation domination [jCarrl . Il975|1 : 



Mi 



PBH 



M, 



H 



10 



15 



10 



-23 



(1.15) 



The PBH mass spectrum depends mainly on two characteristics of the early universe: the 
equation of state w, which determines how large the amplitude of initial inhomogeneities 
should be to halt the background expansion and recollapse, and the nature of th e initial den- 



sity fluctuations, which determines how likely such amplitudes are. ICarrl ||1975|l determined 
the threshold amplitude 5 th = (<5 p ) t h required for the density perturbation to collapse to a 
PBH to be §th ~ w. In this case, one needs perturbations to the FRW metric with mean 
amplitude of order unity to form a significant number of PBHs. 

The special characteristic of PBHs is that they can form at very early epochs and have 
very small masses. The smalles t PBHs would have formed at the end of the inflationary 



expansion [ 
ILvthet al.. 


uarr & 
20061: 


Liasey,. 
Zaballa et al. 


I, even 
.120071 



The mass of the horizon at the end of inflation is 
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BZaballa et all 1200711 



M, 



10 17 g 



10 7 GeVY 



RH 



(1.16) 



where the reheating temperature Trh depends sensitively on the model of inflation con- 
sidered. I n the canonical slow-r oll inflationary model this temperature can be well above 
10 10 GeV ||Kolb & Turned. LL990J]. Taking on account the production of dark matter candi- 
date particles in supersymetric models, this temperature could be dropped by several orders 
of magni tude, however, leptogene sis does not alow the reheating scale to be smaller than 
10 9 GeV IIBuchmuller et all 1200511 . This in turn means that PBHs could have been produced 
with masses much smaller than 10 11 g. On the other hand, PBHs that formed at 1 s have 
masses of order 1O 5 M which is already in the range of masses of black holes at the centre 
of galaxies. 

The small masse s of PBHs promp ted the investigation of their quantum properties. The 
well known result of 



Hawking 



1974|1 shows that black holes radiate with a temperature 



T ~ 10" 



M 

Mr 



Kelvin 



(1.17) 



and evaporate entirely on a time scale 



L evap — 



10 



(SI 



M 

Me 



(1.18) 



where M P) is the sola r mass. With the age of the universe estimated as 1.37 ± 0.015 x 10 10 y 
llSpergel et all 1200711 . we can predict that PBHs with mass M cvit = 5 x 10 14 g are evaporating 
now. PBHs are also the only type of black holes for which the effect of Hawking evaporation 
could be observed. In deed, the black holes ev aporating now would be producing photons 
with energy 100 MeV [|Page & Hawking. 1 1 97 ql . The observed 7-ray back ground radiation at 
this e nergy implies that the density parameter of such PBHs must satisfy llPage & HawkingL 
193 



fi PBH (M ~ 10 15 g) < 10 



(1.19) 



This bound remains the tightest constraint to the abundance of PBHs. Additional cosmolog- 
ical bounds to the mass fraction of PBHs are reviewed in Chapter|4] 

The mass fraction of the universe turning into PBHs of mass M at the time of their 
formation is denoted by Ppbh (M) . This is equivalent to the probability of formation of PBHs 



of mass M. In a rough calcu 



ation, 0pp;u(M) is given by the Press-Schechter formalism 



jPress & Schechten . ll974l : ICamll975|l as the integral of the PDF over all amplitudes 5 P above 
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the threshold 5 



th- 



(3?bh(M) = 2 / F(S p )dS p , 
>s th 



(1.20) 



where the factor two has been added to account for the half volume of the universe that is 
necessarily underdense. With this factor the Press-Schechter f ormula gives a good fit to the 
results of N-body simulations for the case of galactic haloes [|PeeblesL 1 198Qf| . For the case 
of PBHs, an upper limit of integration is formally required. This is the amplitude of an 
inhomogeneity for which the total mass would form a separate closed universe. However, 
the contribution of higher values to the probability is almost negligible and we do not include 
an upper limit he re. For the case of a Gaussian PDF with variance E P (M) this integral is 
approximated by [ICarrl . 1 197511 



Ppbk(M) « 5th exp 



51 



2£2(M) 



(1.21) 



This equation demonstrates the sensitive dependence of the probability of PBH formation 
with 5 th . The above integral is expected to be small due to the exponential dependence on 
the threshold value 5 t h- Ppbu is also known to be small because it is related to the current 
density parameter f2 PBH of PBHs formed at time t and with mass M by 



^PBH — /^PBH 



a(t) 



^ 10 6 /?PBH 



t 

Is 



-1/2 



10 18 # 



PBH 



M 
10 15 g 



-1/2 



(1.22) 



where = 8 x 10~ 5 . The factor a~ l arises because PBHs form mostly during the radiation- 
dominated era but PBH density scales as a~ 3 , while radiation scales as a -4 . From this 
relation we see that any limit on f2pBH places a direct constraint on /3pbh- For example, from 
the bound in Eq. (11.191 ), we infer that /3pbk{M = 10 15 g) can only have a small value of order 

io- 26 . 



1.4.2 Shortcomings 



The simple picture of PBH formation described above has several shortcomings 

1 . In the radiation era the inhomogeneities forming PBHs must have a large amplitude 
when they enter the horizon and they must be bigger than the horizon for a consid- 
erable period of their evolution. As we will show in Chapter [2l the inhomogeneities 
at superhorizon scales are best described in terms of curvature perturbations because 
they are constant in this regime. The curvature perturbatio n has already been used 
in the more recent numerical simulations of PBH forma tion llShibata & Sasakil Il999i : 

1999j 



Niemeyer & Jedamzik . 



cent works on the subject IIYokoyamal . 



Polnarev & Muse 



199 



9; 



„, 



20071 ] . Here, as in several other re 



Green et al 



. 120041 : IZaballa etal . 



2007 
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Josan et all 1200911 . we compute the probability of formation of PBHs from the statis- 
tics of the curvature perturbations. This has the advantage of relating the formation of 
PBHs directly to the initial perturbation spectrum. Additionally, it avoids the gauge 
anomaly associated to the matter density fluctuation. 

2. In the calculation of the probability of PBH formation, one could argue that the Press- 
Schechter formula in Eq. (11.201) is only an empirical a pproximation. Alte rnative ap- 
proaches have therefore considered the theory of peaks [|Green et all l2004f| . However, 
this does not render significant corrections to the Press-Schechter result. Moreover, 
the Press-Schechter formula can be used to c alculate the probabilit ies of large-scale 
structure formation from non-Gaussian PDFs I Matarrese et all 2000 1. Indeed, the lat- 
est numerica l simu lations confirm that it is a good approximation even in this case 
IIGrossi et all 1200911 . This justifies our choice of the Press-Schechter formalism to ex- 
plore new aspects of the probability of PBH formation. 

3. A severe oversimplification of the usual calculation of the probability of PBH forma- 
tion is the assumption of Gaussianity. The exponential decay of the Gaussian PDF is 
preserved after its integration in the Press-Schechter formula (11.211) . The fact that the 
mass fraction involves an integration over the tail of the normal distribution, where the 
probability density is small, leads us to consider that a slight variation on the profile 
of the PDF might modify this picture significantly. Indeed, non-Gaussian probability 
distributions have been co nsid ered in s tudies of the probability of PBH formation by 



Bullock & Primackl 1199711 and 



Ivanov 



[| 1998(1 • The discrepancy in their results and the 
large departures from the Gaussian case make this problem worth revisiting. One main 
objective of this thesis is to derive the modifications that non-Gaussian PDFs bring to 
the probability of PBH formation in the most general cases. We explore for the first 
time the modifications that a non-Gaussian PDF may bring for the bounds on the ampli- 
tude of fluctuations and the higher order statistics parameter / NL on the cosmological 
scales relevant to PBH formation. 



4. The last important problem in the calculation of /3pbh is the determination of the pre- 
cise value of the threshold amplitude 5 th or £ th for the density or the curvature inho- 
mogeneity. This approximation of Ppbh prompted several studies of PBH formation to 
determine the precise value of the threshold amplitude. Early numerical simulations 
of gravitational collapse, however, already showed that t his value depends sens itively 



on the shape and profile of the initial configuration 6„(x) IINadezhin et all 1 197811 . This 



dependence indicates that the lower limit of the integral (11.201) is not uniquely pre- 
scribed for all configurations collapsing to form PBHs. The problem then is how to 
differentiate profiles of initial inhomogeneities in the calculation of the probability of 
PBH formation. This is another problem we address in this thesis. We calculate the 
probability of PBH formation by taking into account the radial profiles of initial cur- 
vature inhomogeneities. This represents a first attempt to incorporate profiles into the 



1.4: Primordial black holes 



24 



calculation of /? PB h and allow for a more precise estimation of the probability of PBH 
formation. 



1.4.3 Alternative mechanisms of PBH formation 



The formation of PBHs is not limited to the collapse of overdensities. PBHs may also form 
at the phase transitions expected in the early universe. Let us here briefly review other known 
mechanisms of PBH formation. 

• PBHs may form at early phase transitions where the equation of state is soft for a small 
period of time. In such transitions, the effective pressure in the universe is reduced due 
to the the formation of non-relativistic particles. Hydrodynamical simulations show 
that at such a phase transition the value of <5 t h is reduced below the value pertaining 
to the radiation era. This mechanism enhances the proba bility of PBH formation at a 



mass sea 



JedamzikLll997h . 



e of th e order of the horizon mass at that time [Kh lopov & PolnarevL LL980; 



Loops of cosmic strings can collapse to form PBHs. Cosmic strings are topological 
defects formed at the phase transitions in the very early universe. Closed loops can 
be formed from string self-intersection. The scale of a loop will be larger than the 
Schwarzschild radius by a factor (Cr//) -1 , where \i is the string mass per unit length, 
a free parameter in the theory. In the cosmic string scenario, these loops are respon- 
sible for the formation of cosmological structures if (G/-0 is of order 10~ 6 . In this 
scenario, there is always a small probability that particular configurations, in which 
all the loop dimensions lie within its Schwarzschild radius, can colla pse to form blac k 
holes. This mechanism has bee n discussed by many autho r s (see e.g. iHawkingl 11198911 : 
Polnarev & Zembowiczl ll 199 ill : loarriga & Sakellariadoul Il993k However, WMAP 
and observations of galaxy distributions show that cosmi c strings can a t mos t con- 
tribute to 10% of the temperature anisotropy in the CMB [|Wyman et all l2005|l . The 
mass per unit length is less constr ained by the observational limits on primordial black 
holes jCaldwell & Casperl Il996|l . Because the /i parameter is scale-invariant and its 
most stringent limit comes from CMB observations, we can say that the formation of 
PBHs from cosmic string loops is subdominant with respect to the standard picture of 
collapse of overdensities. 

One can also consider closed domain walls which form black holes. Domain walls 
are hypothetical topological defects of higher order. In a phase transition of second 
order, such a s might be associated with in flation, sufficiently large domain walls may 
be produced llCrawford & Schramm!. 1 1 98211 . This lea ds to the formation of PBHs in the 
lower end of the range of masses [|Rubin et al. 



• Recently, a mechanism to form PBHs as the result of warping cosmic necklaces has 
been suggested. These topological defects arise in the process of symmetry breaking 
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in the framework of quantum strings IIMatsudal 1200611 . 



In all these mechanisms the PBHs have mass of order the horizon mass at phase transitions 
in the early universe. They are also expected to produce PBHs with a Gaussian distribu- 
tion. Here we are interested mostly in PBHs with a non-Gaussian distribution in order to 
produce constraints on models of inflation, so we do not study these alternative formation 
mechanisms. 



1.5 Thesis outline 

Chapter [2] presents a study of non-Gaussianity from inflationary scalar perturbations. It first 
introduces the relevant definitions and the main tools used in the study of inflationary pertur- 
bations. It then focuses on the derivation of non-Gaussian correlation functions. Specifically, 
the three-point correlation is studied in models where an auxiliary scalar field during infla- 
tion is responsible for the generation of non-adiabatic fluctuations. The cases in which the 
non-adiabatic fluctuations may generate large values of / NL is considered in detail. 

The method used to derive the non-Gaussian correlators requires the solution of the 
Klein-Gordon equation beyond linear order. This equation is solved considering a pertur- 
bative expansion of the nonlinear terms without taking on account the metric back-reaction. 
For the cases in which analytic solutions are possible, the derivation of the three-point cor- 
relation is presented. Finally, the observational limits on /nl are used to constrain models of 
inflation which include a curvaton field, a special case of an isocurvature field. 

Chapter [3] discusses the decomposition of the curvature perturbation 1Z into harmonics. 
This is a technical step, which is necessary in order to write down a path integral for the 
PDF F(1Z). We present the calculation for the Gaussian case first, in order to clearly explain 
our method with a minimum of technical details. This is followed by the equivalent calcula- 
tion including non-Gaussian corrections which follow from a non-zero three-point function. 
Finally we calculate the probability F[TZ(k)], which will be used to derive a non-Gaussian 
probability of PBH formation. 

In Chapter |4] we compute the mass fraction /3pbh resulting from a non-Gaussian PDF 
of primordial curvature fluctuations 1Z. We restrict ourselves to the case in which the non- 
Gaussian PDF corresponds to a constant value of / NL . It is first shown ho w to reconcile the 



discrepancy b e tween two previous studies of non-Gaussian PBH formation [|Bullock & Primackl . 



1997 



ep 

; I 



IvanovL 1 199811 . We then calculate the modifications to the observational bounds to 
/?pbh when a large value of /nl is included. 

Chapter [5] explores the probability of finding non-trivial spatial profiles for the perturba- 
tions that form PBHs. The numerical simulations show that the usual assumption of homo- 
geneous spherically symmetric perturbations collapsing to PBHs is not appropriate. Chapter 
|5] provides a probabilistic analysis of the radial profiles of spherical cosmological inhomo- 
geneities that collapse to form PBHs. Based on the methods used to construct non-Gaussian 
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PDFs, we derive the probability distribution for the central amplitude of 1Z and for the sec- 
ond radial derivative d 2 TZ/dr 2 at the centre of the spherically symmetric inhomogeneity used 
to describe the radial profiles explored in studies of gravitational collapse. We then consider 
the joint probability of both parameters to compute the correction to Ppbh- The results show 
how much the probability of PBH formation can be reduced if we do not include all possible 
configurations forming PBHs. 

Chapter[6]is the summary and conclusion of this thesis. We also describe future research 
which may follow. The key achievement of this thesis is to combine for the first time the 
study of two crucial probes of the early universe. The effects of nonlinear non-Gaussian 
inhomogeneities and primordial black hole formation. 



Chapter 2 

Non-Gaussian curvature perturbations 



2.1 Outline 

Observations of cosmological structure and CMB parameters are best interpreted in the con- 
text of cosmological perturbation theory. This is a useful tool to connect observations with 
models of inflation derived or motivated by high energy physics theories for which there is 
no other available test. Surprisingly enough, the simplest inflationary model, consisting of 
a single scalar field slowly rolling down a quadratic potential, motivated mainly by its sim- 
plicity, has passed all observational tests. The future of cosmology relies on the extension of 
experimental tests and predictions for new cosmological parameters, mostly beyond linear 
order. This is crucial if we want to achieve a better understanding of the physics dominating 
the early universe. 

This is enough motivation to study the nonlinear regime of cosmological inhomogeneities. 
Among the observable effects, the non-Gaussianity of perturbations has been widely studied 
in inflationary models. Non-Gaussianity is an important observational test as it might elimi- 
nate models of inflation even for a null detection. Our goal in this chapter is to compute the 
nonlinear correlations of a general isocurvature field which is valid for all models. 

We first introduce the theory of perturbations and then focus on the situation in which 
the curvature perturbation is generated by the quantum fluctuations of an isocurvature scalar 
field. The isocurvature or entropy perturbations are transformed into curvature inhomo- 
geneities at the end of a period of inflation or shortly after it. We will show that only the 
presence of entropy fluctuations can affect the evolution of curvature fluctuations on super- 
horizon scales. 

At linear order, we will show under which conditions the observed power spectrum of 
curvature fluctuations can be attributed to the action of the isocurvature field. Subsequently 
we present a method of deriving such correlations from the solutions to the Klein-Gordon 
equation of the isocurvature field. For specific cases we are able to derive an explicit expres- 
sion for the nonlinear parameter / NL . The prospects of observationally testing the predictions 
for the models of structure formation presented here are also briefly discussed. 

The introductory sections of this chapter present a review of the elements of the stan- 
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dard inflationary scenario, including the linear perturbation theory. We present the relevant 
definitions and conventions to be used, with particular attention to those results of linear 
perturbation theory which will be used in this and subsequent chapters. From Section [231 
onwards, we focus on the description of the non-Gaussian correlators of an auxiliary isocur- 
vature field x- The expressions for the curvature perturbation three-point correlators and the 
/nl values are presented in the last section of this chapter. 



2.2 Linear perturbations 



In cosmological perturbation theory, the universe is described to a lowest order by a homo- 
geneous, isotropic background spacetime. The large-scale inhomogeneities and anisotropies 
observed in the real universe result from the growth of dens ity fluctuations, the amplitudes 
of which are small in the early stages of the universe. (See Peebles! [1 198QTI for a textbook 
description of the development of perturbation theory.) 

In the framework of perturbation theory, the homogeneous background spacetime is ac- 
counted for by an ansatz metric. The most useful ansatz in this case is the Friedmann- 
Robertson-Walker (FRW) metric: 



-1 

-Hj 



(2.1) 



where the conformal time 77 is given in differential form by 

dt 



drj 



a(t)' 



(2.2) 



a is the scale factor and 7^ is the metric of the three-dimensional space. In our notation Greek 
indices have values 0, 1, 2, 3, while Latin ones have value s 1, 2, 3. We assume t hroughout a 
flat space, relying on the observational limit \Q K \ < 1CT 2 IIKomatsu et all 1200811 . The FRW 
metric describes the isotropic space-time expanding at a uniform rate. The expansion rate is 
conventionally characterised by the Hubble parameter 



H 



din a 1 din a 1 



df 



a dr] 



H, 



(2.3) 



where H is defined with respect to coordinate time t and H with respect to conformal time 
77. 
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2.2.1 Metric perturbations 

In perturbation theory, observed anisotropics and inhomogeneities are considered as depar- 
tures from the metric (12.11) . For a perturbed metric, the metric tensor can be split as 

9v> = 9$ + ( 2 - 4 ) 

where gf} is the homogeneous FRW background and 5g^ v encodes the perturbed quantities. 
First order scalar perturbations of the metric are expressed in terms of the functions </?, B, ip 
and E, which are defined by 

- 2aV( 7 7: x ) ! 
a 2 B(r),yt) jh 

2a 2 (ip(rj,x)j ij + E tij (T},x)), 

where the index (s) denotes scalar modes. The vector constructed from the scalar B is 
necessarily curl-free, i.e. Biyi = 0. The pure vector contributions to the metric perturbations 
are 

8$ = -<?S u 5g$ = 2a 2 F {ijj) , 

where we demand Suj] ^ 0. The symmetric derivative of the function F 4 is the vector 
contribution to g^ . To distinguish scalar and vector contributions, the vector part is forced 
to be divergence-free, i.e., ^ Si j = 0. (The decomposition of a vector field into curl- and 
divergence-free parts is formally known as Helmholtz's theorem.) The tensor contribution 
to the perturbation quantities is 5 gf) = a 2 hij. This is constructed as a transverse, traceless 
tensor, which guarantees that it cannot be constructed from scalar or vector perturbations. 

The perturbation functions if, B, ip and E, represent four degrees of freedom. The diver- 
genceless vectors Si and Fj each have two degrees of freedom and the transverse traceless 
tensor hij has two more. We therefore have 10 degrees of freedom in total. The contravari- 
ant metric tensor of the perturbed metric is constructed, to first order, from the condition, 
9^9^ — &a- Fi nai ly> the line element of the metric is 

ds 2 = a 2 (r]){ - (1 + 2cp)dri 2 + 2(B ti - S i )dr ] dx i 

+ [(1 + 2^) 7ij + 2E ij + 2Fij + h^} dx l dx j }. (2.5) 

In the present work we will study the nonlinear perturbations as the quantum fluctuations 
of scalar matter fields. We will establish the correspondence between scalar matter fluctua- 
tions and scalar perturbations in the metric at first and second order in perturbation hierarchy. 
We will then derive statistical parameters of nonlinearity. 

In contrast to the scalar metric fluctuations, the vector and tensor perturbations in the 
metric are not sourced by scalar matter perturbations at first order. In the standard picture, 



Sg {s) 
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they a re only related at second or higher order in perturbation theory (see e.g. iLu et al 
|2008|]), therefore their contribution to the statistical parameters of nonlinearity are sub- 
dominant and henceforth we neglect their contributions to the perturbations in the metric. 



2.2.2 Gauge freedom 

In general relativity, the mathematical relations between physical quantities are manifestly 
independent of the coordinate choice. However, there is no covariant way of splitting back- 
ground and perturbed variables. There is always an unphysical coordinate or gauge de- 
pendence associated with perturbed spacetimes. T his issue of gauge ambiguity was disre 



garded in the initial works of perturbation theory [|Lifsh itz. 



1946 



Lifshitz & Khalatnikov 



19631 . This cou 



way by iBardeen 



d lead to erroneous results which were eventually resolved in a systematic 
[| 1980(1 - The importance of determining the gauge changes that equations 
and perturbations undergo leads us to look at this problem in detail. In the following we 
adopt a 'passive ' appro ach to gauge transformations (For a recent review of these results, see 



Malik & Wands 



Q2008I]). Let us consider the general coordinate transformation, 



(2.6) 



where £° = £°(?7, x l ) is a scalar that determines the choice of constant-?/ hyper surf aces. The 
scalar £ and the divergence-free vector are also functions of the original coordinates within 
these hypers urf aces. 

The principle of relativity states that any physically meaningful measurement must be 
invariant for all observers, in particular, for observers with different coordinate systems. One 
of these invariants is the line element ds 2 , where coordinates enter via differentials. Such 
differentials and the scale factors in both coordinate systems are related in the following 
way: 



d v =dfj-t ' dfj-Q dx\ 

dx % = d? - (£'; + c) dfj - + cj) 

a(rj) =a(f)) -£°a'(fj). 



(2.7) 



Where a ' is the derivative with respect to conformal time. To first order in the metric per- 
turbations and coordinate transformations, the perturbed line element, Eq. (12.51) is written in 
the 'shifted' coordinate system as 



ds 2 = a 2 (f,){ - (i + 2(</? - ne - e)) dv 2 +2[(B+e- o.< -Si+ g ^ 

+ [ (1 + 2(^ - 7^°)) Ty + 2 (E - + 2F id - 2&J + h tj \ dx j dx l } , (2.8) 
where, as before, H = a'/a is the Hubble parameter in terms of conformal time. This metric 
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can also be written using the initial definitions in terms of the 'shifted' coordinates: 

ds 2 = a 2 {fi){ - (1 + 2y)df} 2 + 2(B ti - S^dfjdx 1 (2.9) 
+ [(1 + 24>)% + 2Eij + 2F id + hij] dx { dx j }. (2.10) 

This shows that the coordinate transformation Eq. (12.71) induces a transformation of the met- 
ric perturbations. Comparing Eqs. (12.51) and (12.101) . the change is given to first order by 

$=i[,-ng>, (2.ii) 
<p=<p-ne-t ', (2.i2) 

B=B + f?-£', (2.13) 
E=E-£. (2.14) 

It must be stressed that the gauge transformations are, in effect, a change of the correspon- 
dence between the perturbed spacetime and the unperturbed background spacetime. 

A first exercise concerning gauge transformations is to find those quantities which re- 
main invariant after a gauge transformation. To first order in perturbation variables, gauge- 
invariant quantities are linear c ombinations of the gauge-dependent quantities presented 



above. For scalar perturbations iBardeenl lll980ll shows that only two independent gauge 



invariant quantities can be configured purely from the metric perturbations: 

$ B =cp + H(B-E') + (B-E')', (2.15) 
■$ B =-il ) -H{B-E'). (2.16) 

Any other gauge invariants in the metric are linear combinations of these two quantitie s 
because the gauge freedom allows only two arbitrary scalar functions £° and £ IIMalikl. 1200 ill . 



The Bardeen invariants will be useful in relating curvature perturbations in different gauges, 
as we will show below. 

2.2.3 Perturbations of the matter sector 

Before displaying conservation equations for the curvature perturbations, we will discuss the 
perturbations of the matter sector. For a perfect fluid, that is, a fluid with no heat conduction 
or viscosity, the stress-energy tensor is 

Tt = (p + p)u»u u +p5Z, (2.17) 

where the 4- velocity is defined with respect to proper time r as 

u» = ^ (2.18) 

dr 

and is subject to the normalisation u^u^ — —I. Anisotropic stresses would be encoded in a 
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stress tensor 11^, but are absent for perfect fluids and for scalar fields minimally coupled to 
gravity. These are precisely the kinds of matter considered here, so we ignore the tensor 
in the subsequent analyses. 

Using the normalisation = — 1, the perturbed velocity has components 

u° = -(l-<p), + (2.19) 

a a 

u = -a(l + ip), Ui = -a(vi + v ti + B ti - Si), (2.20) 

where the spatial parts are written in terms of the gradient of a scalar v yi and a (solenoidal) 
vector Uj. The perturbed energy-momentum tensor is: 

r °= -(po + Sp), (2.21) 
T? = (po +Po)(B,i + v { + Vi - Si), T = -(p +p )(y> i + v% (2.22) 
T) = (p„ + 5)6 l jy (2.23) 

where p and p represent the uniform pressure and matter density. In general our scalar 
stress-energy components can be written as f(r), x l ) = fo(r]) + 5f(r), x % ), with the subscript 
denoting the background homogeneous part. As in the case of metric perturbations, coordi- 
nate transformations will affect the matter perturbations. This means that the matter density, 
velocity and pressure perturbations are gauge dependent. Under the transformation Eq. (12.61) . 
perturbed scalar functions of the form are thus transformed as 

Sf = Sf- io£°. (2.24) 

The vector perturbations are derived either from a potential, which will transform with the 
shift or from a pure divergence-free vector, whose transformation depends on £\ In 
particular, the velocity potential v transforms as: 

v = v + £, (2.25) 

and the vector function v % is transformed as 

v* = v l + C'. (2.26) 

2.2.4 Physical quantities and scales 

Before addressing the characteristics and governing equations of the perturbed spacetime, 
let us define the physical scales and the quantities that determine of the size and age of the 
universe. 

The time-like 4- vector field 

Np = -a{l + <p)S°, (2.27) 
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defines the direction perpendicular to the hypersurfaces of constant time. In consequence, 
this vector field defines a coordinate system. This vector is unitary (N U N U = — 1) and the 
contravariant vector N u = g Ufl N^ has components 

N° = -(l-ip), N^-iS'-B' 1 ). (2.28) 
a a ' 

The expansion rate of the spatial hypersurfaces with respect to the proper time of observers 
with 4- velocity is 9 = N^.^. Considering only scalar perturbations, this is given by 

6 = 3^(1 -<p) + 3 V - -V 2 (B - E') , (2.29) 

or a a 

where the operator V denotes the usual three-dimensional gradient. 

By looking at the relation between proper and coordinate time, dr = (1 + ip) dt, we 
extract the expansion with respect to coordinate time from the above expression by writing 

(t) = (l + ^)0 = 3# + 3^ + W (t) , (2.30) 

where the Hubble parameter, H = a/ a is the background uniform expansion rate with re- 
spect to the coordinate time. The shear scalar in coordinate time is o"( t ) = (E — B /a). 

For the sake of completeness, we include here the definition of some useful scales. A 
comoving observer is one moving with the expansion of the universe, i.e., one who measures 
zero net momentum density. The distance of a comoving point from our location (taken to 
be at the origin of coordinates), is given by r(t) = a(t)x, where x is the comoving distance. 

The Hubble radius r H = H^ 1 provides a good estimate for the distance light has travelled 
since the Big Bang. Formally, the integral 

r* ds 

V- / -pr, (2-31) 
Jo a i s ) 

defines the comoving distance travelled by a free photon since t — and until time t. This 
is important because no information could have travelled further than 77. This define the 
'comoving particle horizon'. In the above integral rj can be taken also as the conformal time. 
In a matter dominated universe 77 oc a 1 / 2 , while in radiation domination 77 oc a. In a de Sitter 
inflationary universe 

f ds f da H^ 1 
V — I ~ TT = / TTT = TIT- (2-32) 



a(s) J Ha 2 a(t) 

This shows that in an inflationary phase 77 — > —00 as the universe approaches the initial 
singularity a = 0, and increases monotonically towards 0. This leads us to consider the mag- 
nitude 1 7] I when we use the conformal time in our calculations for inflation. The maximum 
distance light travels from time t = to us is simply the comoving horizon times the scale 
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factor, 

/•tnow 

d H = a (tnow) / —tt- (2.33) 

Jo a ( s ) 

This is called the particle horizon, i.e., the radius of the region which is in causal contact 
with us. This equation can be applied to find the horizon radius at times different from £ now . 
Considering the current dark energy domination and a cold dark matter component (with 
density parameter f2 m ), the horizon size does not coincide exactly with the Hubble scale. 
However, an approximate solution to the integral Eq. (12.331) shows that, 

d H (t aow ) « 2H^ w 1 + °-°^l nnm ~ 3.5H^, (2.34) 



where the last expression assumes Vt m = 0.25 [|Hu et al.Lll998|1 . 



2.2.5 Particular gauges 

Let us now focus on the expressions for the curvature perturbation on three useful choices 
of time slicing and threading. These are the uniform curvature gauge, the uniform density 
gauge and the comoving gauge. 

For the first case the spatial hypersurfaces present an unperturbed 3-metric, which means 
ip = E = 0, in other words, curvature perturbations of the three-metric are set to zero. 
We distinguish the quantities written in this gauge with a subscript n indicating a constant 
curvature. So for a general coordinate system we require the following transformations: 

^Lconst = yr, £« = E - (2.35) 

In this case, the scalar perturbation becomes 

8f K = V- f^. (2.36) 



In particular, for scalar fields, this is the gauge-invariant Sasaki-Mukhanov variable [|Sasakil . 



1986l : lMukhanovLll988h . explicitly, 



50 K = 50-0^. (2.37) 

In the uniform density gauge, the requirement 5p = for constant-time hypersurfaces 
implies 



$ P = J- (2-38) 
The gauge-invariant curvature perturbation on these hypersurfaces is denoted by ( and de- 
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fined as 



C = i'Sp = ip - H 



Sp 



(2.39) 



In this case, there is another degree of freedom and one can pick either B, E or v to be zero. 
The gauge invariance is made explicit when the curvature perturbation is written in terms of 
the Bardeen variables. We will return to this point once we have defined the curvature in the 
comoving gauge. 

The comoving gauge is subject to the condition that the spatial coordinates comove with 
the fluid, that is, for a constant-time slice the 3-velocity of the fluid vanishes, v 1 — 0. The 
threading is chosen so that the constant-?] hypersurfaces are orthogonal to the 4-velocity w M , 
which demands v+B = 0. An immediate consequence of this choice of gauge is that the total 
3-momentum vanishes on constant-time hypersurfaces. For this reason several authors call 
this gauge the zero-momentum gauge. Using Eqs. (12.251) and (12.131) . the chosen conditions 
imply that 



(2.40) 



with £(x % ) the residual coordinate gauge freedom. This quantity is not specified at this stage 
because it is not required for the determination of the scalar quantities like curvature, expan- 
sion and shear. For arbitrary coordinates, the scalar perturbations in the comoving orthogonal 
gauge are given by 

<p m = <p+-[(v + B)a]' 1 ^ m = ^ + H{v + B), E m = E+ [ vd V -i. (2.41) 
a I 



The scalars <p m and i]) m defined in this way are gauge-invariant. The density perturbation in 
the comoving gauge is also given in gauge-invariant form by 



5p m = 5p + p'(v + B). 



(2.42) 



Some authors use the gauge-inva riant dens i ty perturbation in the comoving gauge by defining 
the combination A = 5p mi / p | Bardeenl . ll980l : lKodama & Sasaki 1 19841. 



The curvature pe rturbation tjj in the com oving gauge was first used by LLukashl [[1980] and 
first denoted as 1Z by lLiddle & Lythl [119930. It is mathematically defined as 



K = ip m = ip + H(v + B). 



(2.43) 



In the next section we will find that, through the Einstein equations and gauge-invariant 
quantities ^ B and defined in Eqs. (12.161) and (12.151) . one can establish an equivalence 
at large scales between the curvature perturbation to linear order in the uniform density 
gauge and the same perturbation defined in the comoving gauge. Taking this equivalence for 
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granted, in the meantime, allows us to relate 71 and 5p directly. Indeed, if we consider the 
gauge transformation (12.391) from an initial fiat hypersurface, then 



n = c = - 




p 



at scales beyond the cosmological horizon (as will be made explicit below). 

This last transformation shows the way of avoiding the gauge anomaly. One can always 
change the gauge (or frame of reference) and establish the equivalence between the perturba- 
tions of any two gauges as long as a particular gauge is chosen at the start and all quantities 
are initially defined in this gauge. 



2.3 Evolution of perturbations and conserved quantities 

Just as the perturbed scalars in the metric are gauge dependent, so are the evolution equa- 
tions for these quantities and they must be treated carefully in order to avoid spurious gauge 
modes. 

The equations governing the dynamics of space-time are found by varying the action X 
with respect to the metric and matter components. The action is defined as 

7-1 



1= Cy/=g<Fx, (2.45) 

J — oo 

where £ is the Lagrangian density of the matter and the gravitational field. The gravitational 
Lagrangian of general relativity is, 

C = -R/2, (2.46) 

where R is the Ricci scalar. The Lagrangian density for a classical matter field minimally 
coupled to gravity is 

C = K — V — R/2, (2.47) 

with K the kinetic energy and V the potential energy. For the matter sector of the Lagrangian 
we can define the energy-momentum tensor as 

dC 

T„„ = -2— + g,X. (2.48) 

In particular, the energy-momentum tensor for a perfect fluid with density p and isotropic 
pressure p and 4-velocity m m is given by Eq. (12.171) . 

Let us now look at the Lagrangian density of a single scalar field 0, minimally coupled 
to gravity. Its kinetic energy is K = —l/2g^ v d l _ l (J)d u (f), so from Eq. (12.471) we find a canonic 
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action 



C M = ~ [g^d^d^ + 2V{<j>)} . (2.49) 

Using the definition Eq. (12.481) . it is easy to show that the scalar field energy-momentum 
tensor is 

T; = g a »<t>,^, a - 5; (v{<j>) + lg a % a <f>,^ ■ (2.50) 
The comparison between the last expression with Eq. (12.171) . T ji u for a perfect fluid, show s 



that we can define the density, isotropic pressure and velocity as IITabensky & Taubl . 



197311 



\g» u < 



„r-ri, p = -<T<1>,^ + v, p = -gTM,u-V. (2.51) 



This identification provides an easy way to quantify the energy density of a scalar field and 
its perturbations. The Lagrangian density for two real fields is 

Cm = -^GT<MV) " l^X^u) - U(<f>, x), (2.52) 

which features the joint potential U (0, x) of the participating fields, each minimally coupled 
to gravity. This case is what concerns us in the rest of the chapter and we shall focus on its 
dynamical equations. 

2.3.1 Background equations 

The Einstein equations are found by varying the action (12.471) with respect to the metric. 
Under regular conditions, with no variations of the fields at the boundaries, the equations are 
found by applying the operator 



6 x 5 



Sg "5{d,g) 



(2.53) 



to the Lagrangian. The Einstein equations dictate the dynamics relating the local spacetime 
curvature to the local energy-momentum. In the adopted natural units, 

= T^, (2.54) 

where the left-hand side is the Einstein tensor, defined as 

G^v = Rfiv — -g^R- (2.55) 

The Einstein equations can be split into components that are parallel or orthogonal to the 
time-like field at any order in perturbation expansion. The two independent equations 
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obtained at the background level are the Friedmann and acceleration equations: 

H 2 = l -a 2 p 0) (2.56) 

H' = -]:a 2 (3p + p ). (2.57) 
o 

Additionally, the Bianchi identities G^.^ = imply the local conservation of energy and 
momentum, 

= 0, (2.58) 

where ; denotes a co variant derivative with respect to the metric g^. 

For the background quantities, the energy-momentum conservation equations provide an 
expression for the expansion in terms of the matter fields. For the case of a single fluid in the 
background FRW universe, the equation T u . u = gives 

p' = -3H(p + Po). (2.59) 

Note that the isotropy assumption means there is no net background momentum and thus no 
other conservation equation at zeroth order. Moreover, Eq. (12.591) can also be obtained as a 
combination of the Einstein equations (12.561) and (12.571) . 

The homogeneous Einstein equations can be solved for the variables a(t) , p(t) , p(t) when 
an equation of state for the matter components is provided. This is dictated by the micro- 
physics of the matter. In particular, the equation of state 

p = wp, (2.60) 

describes most of the relevant cases of the post-inflationary cosmology. For the case of pure 
radiation w = 1/3, while for pressureless dust w = 0. For fluids with such an equation of 
state, the solutions to the Einstein equations are 

/ \ -3(w+l) / \ 2/(3iu+5) 

M-) • ■ < 2 - 61 > 

Pi V a i/ fl i \ViJ 

with initial conditions p = p\ and a = a x at r] = rji. 

When more than one fluid is present, we account for the contribution of each component 
to the total matter density by defining a dimensionless density parameter 

%) = -g^r, (2-62) 
where the factor 3H 2 /a 2 is the critical density. In the flat universe that concerns us, the 
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curvature contribution Q K is zero and the sum of all matter contributions is unity, i.e., 



(note that the dark energy that dominates the expansion in the present stage of our universe 
should also be included in this sum. This component is customarily denoted by With 
these definitions Eq. (12.561) can be written in the form 



For the case in which the matter is dominated by a single scalar field 0, energy density 
conservation leads to the Klein-Gordon equation, 




(2.65) 



which can also be derived from the variation of the action with respect to <p. The potential 
V is assumed to be an explicit function of the field alone. In the case where there is more 
than one field, a Klein-Gordon equation is obtained for every scalar field, with an interaction 
potential U . Note that the Klein-Gordon equation is valid at all orders in the perturbation 
expansion. We will rely on this important fact to derive the contribution of nonlinear pertur- 
bations to the non-Gaussianity of the primordial fluctuations, the ultimate objective of this 
chapter. 

The Klein-Gordon equation for a homogeneous scalar field O f° r a FRW background 
metric is 



The solutions of this equation will be explored in the context of inflation in Section 1241 The 
inflationary behaviour is guaranteed when the dominating scalar field meets the so called 
slow-roll conditions. The dynamics of inflation will be discussed in more depth in the fol- 
lowing sections. In the meantime we note that, as mentioned in Section [L2l the perturbations 
produced during a period of inflation exit the cosmological horizon due to the shrinking of 
the latter scale. In a super-horizon regime, under suitable conditions, the curvature inhomo- 
geneities are time-invariant. 




(2.63) 




(2.64) 



0o + 3#0 O + — = 0. 



(2.66) 
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2.3.2 Dynamics of perturbations 

At linear order, the scalar metric perturbations are related to matter perturbations via the 
Einstein equations. The density and momentum constraints are 

3H(H<p - i//) + V 2 [if> - Ha] = - ^a 2 5p, (2.67) 
W - Htp\ =\a 2 (p + Po) [v + B], (2.68) 
and two evolution equations for the scalar metric perturbations 

i)" + 2Hif>' - H<p' - (2W + H 2 )p = - ^a 2 Sp, (2.69) 
a + 2Ha + ^-<p=Q. (2.70) 

The energy-momentum conservation equations for the perturbed spacetime are related to 
the ones above via the Bianchi identities. Specifically, the evolution for the energy density 
perturbation is 

5p' + 3H(5p + dp) = -(po + po) [3^ + V 2 (v + E')] , (2.71) 
and the momentum conservation equation is 

[(Po + Po)(v + B)}' + 5p= -(p + p )[cp + 4H(v + B)}. (2.72) 

Instead of solving the equations at first order, let us show how the dynamical linear equations 
encode two important implications for cosmological perturbations. A convenient way to find 
conserved quantities is to work in Fourier space (the Fourier transformation is here denoted 
by J 7 .) In this case, a generic coordinate-dependent perturbation fit, x) is decomposed into 
harmonic functions of time: 



CO 



/(t ' x)= (2^/ ex P(- ik " x ) /*(*)> Le -' F[f(t^)} = fkit). (2.73) 

Each function f k (t) is referred as a perturbation mode and labelled by its comoving wavenum- 
ber k and has an associated scale Afc = a(t)/k. The only characteristic scale of the unper- 
turbed universe is the Hubble scale or cosmological horizon as defined in Eq. (12.331) . When 
perturbation modes lie well outside the cosmological horizon, the ratio 

e = d H /\ k = k/a{t)H{t) (2.74) 



is much smaller than one. Since the spatial derivatives V are transformed to k/a, this shows 
that we can neglect gradients terms in the equations compared to the time derivative which 
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scales as H, i.e., for a given function f(t, x) with Fourier transform fk(t), 

k 



-/*(*) 



(2.75) 



This simplifies the equations considerably and shows under which conditions the curvature 
perturbations are conserved in the superhorizon regime. 

Let us use this approximation to look at the equivalence of the curvature perturbation 
on large scales in different gauges. From the Einstein equations we can rewrite the gauge 
invariant curvature perturbation TZ in terms of the curvature perturbation variables. Making 
use of the momentum constraint Eq. (12.681) . we note that 



v + B 



Hip-ip' 
H'-H 2 ' 



(2.76) 



We can insert this in the definition of TZ in Eq. (12.431) and write the latter in terms of the 
Bardeen gauge-invariant quantities. That is, 



n 



B 



13 > 



H'-H< 



(2.77) 



which represents an alternative form of Eq. (12.681) . 

We can follow a similar procedure for Eq. (12.671) and write the uniform density curvature 
perturbation £ in terms of Bardeen invariants: 



H'-H 2 
H 



V' B + H§ 



B 



H'-H 2 
H 



B 



1 

m 



v 2 ^ 



B- 



(2.78) 



Note that, because the curvature perturbations in the last two equations are written in terms 
of gauge-invariant quantities, TZ and ( are manifestly gauge-invariant themselves. Moreover, 
the combination of these equations leads to the gauge-invariant generalisation of the Poisson 
equation, 



V 2 * B = 3 («' - H 2 ) (TZ - = y fa 



(2.79) 



As before, if gradients are discarded, both TZ and £ coincide. This result is important in view 
of the consequent correspondence (12.441) . which is used extensively throughout this thesis. 

A second important feature is the evolution of ( on superhorizon scales. The energy 
conservation Eq. (12.711) can be written in coordinate time as 



S'p + 3H (Sp + Sp) = - (p + p ) 3^ + V 2 (- + E) 



(2.80) 



This leads to an evolution equation for the perturbed energy density when we include Eq. (12.591) : 



[p + 5p\ +3(H + i))\p + 6p + po + 6p] = -V 



E 



(2.81) 
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which in view of the definition of the expansion 9^, Eq. (12.301) . gives to first order 

P + 0(t) [P + P}= V 2 [- + E] + V 2 (a (t) ) [p + p ] + 0(5 2 ). (2.82) 



We emphasise that, in our notation, the density and pressure with no subscript represent the 
sum of the background function and its perturbation, i.e., p = p + Sp. The 0(5 2 ) term 
indicates that the above equation is valid to first order in perturbation expansion. 

Since we are discarding spatial gradients in the equations of motion on super-horizon 
scales, Eq. (12.811) provides an evolution equation for the curvature perturbation: 

i>\Po + Po] = ~-H(p + p), (2.83) 



i = -\-Je- + (2.84) 



or 



3p + Po 3 VG°o + Po) 2 

If we work in a gauge where the time slices are uniform-density hypersurfaces, i.e., the 
uniform-density gauge, we may set 

Sp -> 0, 5p -> Sp Sp , ip -> ip$p = (■ 

All perturbations, including the pressure perturbation in this gauge are independent of the 
density perturbation. From thermodynamics we know that in fluids where entropy is constant 
the pressure is a function of the density. In this category fall the barotropic fluids, defined 
as those fluids in which the pressure is only a function of the density p and vice-versa. In 
general the pressure of a thermodynamic system (in our case the universe) is a function of 
both the density and the entropy, 

5p = c 2 s 5p\ s + 5p\ p , (2.85) 



where c 2 = dp/ dp\ s is the adiabatic sound speed in the system (see e.g. IChristopherson & Malik 



||2009|l for a careful treatment of thermodynamics of fluids in cosmology). If one defines the 
entropy perturbation 5s from the identity 5p\ p = p5s, then one has 

5s = S -E - 5 -P. (2.86) 

p p 

In view of this, Eq. (12.841) for the curvature perturbation reduces to 

Po ( 5p Sp \ Hpo 
C = tt I 7 = ; Ss. (2.87) 



This result is important in the description of the evolution of perturbations after inflation. It 
indicates that for inhomogeneities with characteristic scales much larger than the size of the 
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cosmological horizon, the curvature can on ly be modified if the matter content of the universe 
has a non-adiabatic or entropy component jWands et all 12000(1 . It is also remarkable that this 



argument requires the conservation of the energy-momentum tensor, and not necessarily the 
Einstein equations. This means that the described property is valid for any theory of gravity 
in which energy is conserved. 

This result has been extended beyond linear orde r in the perturbation e xpansion. On 
superhorizon scales, the expression at second order is IIMalik & Wandsl 120041. 



(2 



Po 



3(i>o + Po) 2 



5p : 



2\p 



(Po + Po) 



Sptlp - 2 (p + po) Ci 



Ci 



(2.88) 



where numerical indices indicate the order of each quantity in the perturbation expansion. 
This result shows that, as in the case of the linear (, the evolution at second order depends 
y on the entropy perturbation and its derivatives. This has also been proved to all orders 



on 
by 

i.e., p/p = const, and the sin gle field inflationary case , for which the conserva tion of ( had 



Lyth et all I12005ll . This result generalises the special cases of a constant equation of state, 



1999 



Salopek&BondLll990h . 



been previously been verified UShibata & Sasakil 

Eqs. (12.871) and (12.881) have motivated several studies searching for significant growth 
of £ on superhorizon scales during and after inflation. In particular, theories of multi-field 
inflation have been proposed to generate the the curvature perturb a tion and, at the same 
time, an observab 
Enqvist & Nurmil 



e sign ature of non-Gaussianity IIMollerach , 



1990 



Lvth&WandsL 12002 : 



200511 . We will now study the effects of considering an auxiliary field, 



with special attention to those models where a large non-Gaussian contribution arises. We 
start with a short description of inflation in the next section. 



2.4 Inflation 

One can define cosmological inflation as the epoch when the scale factor of the universe is 
accelerating^: 

a > 0. (2.89) 

This condition can be written in terms of a more physical quantity. The period of inflation 
can be considered as an epoch in which the comoving Hubble horizon decreases with time: 

Jt ( H ~^ < °- (2 ' 90) 

Within general relativity, the above conditions on the time dependence of the scale factor give 
conditions on the matter content through the Einstein equations. In particular, Eq. (12.571) can 



'This does not include the late time acceleration at the current epoch which is attrib uted to dark energy. Whereas 
inflation-like scalar fields may be responsible for such behaviour (see, e.g. jMartir] 1200 8ll ). in this thesis we are not con- 
cerned with the dynamics of the late universe 
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be used to write the condition (12.901) as 

3p + po < 0. (2.91) 

Demanding a positive energy density T° = p is a sensible physical condition. The condi- 
tion T 00 > is known as the weak energy condition. In view of this, the above equation 
demands that the dominant matter must have negative pressure during a period of inflation. 
The simplest matter field with this property is a scalar field, which is composed of spin-0 
particles. The concept of a scalar field is prevalent in particle physics where scalars such as 
the Higgs scalar are essential in the construction of the standard model. Although no scalar 
particle has so far been observed, they play a fundamental role in cosmology, as they possess 
the unusual feature of their potential energy dominating over their kinetic energy. 

As indicated by Eq. (12.511) . the Lagrangian definition of the energy-momentum tensor 
requires that the energy density and pressure for a homogeneous scalar field be 

Po = \<p 2 + V^), PQ = l -^-V{<P). (2.92) 

This shows that in order to meet condition (12.911) we require the potential V(4>) to dominate 
over the 'kinetic' term. This can be dynamically achieved with a sufficiently flat potential 
provided the field is displaced away from its minimum. Such physical conditions are con- 
trolled by two parameters: 

Ml (V'V V" 
e SR = -f (—J , r/sR = M P — . (2.93) 

These are called the slow-roll parameters and, during inflation, they are subject to the slow- 
roll and friction-dominated conditions 

e S R < 1, Vsr < 1. (2.94) 

The first condition, the 'slow-roll' condition, ensures the slow rolling of the field down its 
potential. The second, a 'friction-domination' condition constrains the potential to be very 
flat one for the period of inflation. This condition is imposed to allow for an extended pe- 
riod of inflation (which should last for over 60 e-folds of expansion), required to recover 
a sufficiently flat and homogeneous universe in the observed scales. When both slow-roll 
parameters are much smaller than one, the dynamics of the single field <p guarantees an ac- 
celerated expansion with a shrinking comoving Hubble horizon. 

The scalar field satisfying these properties is called the inflaton. The equations dictating 
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the background dynamics of a universe dominated by the infiaton are 



r2 v{4>) 



H z = -f^, (2.95) 
3H<j)= -V'(<f>). (2.96) 

The homogeneous Klein-Gordon equation (12.661) . reduces to Eq. (12.961) in the slow-roll 
regime. The exact solution to these equations requires the specification of the potential as 
a function of the homogeneous scalar field 4>(t). However, the first equation already shows 
that, if we assume a constant V, 



a(t) oc exp(v / V/3i). (2.97) 

This illustrates explicitly the exponential growth of the scale factor in the inflationary regime. 
Additionally, the field depends only linearly on time to lowest order, as expected for a slow 
rolling field. An expansion in powers of the slow-roll parameters shows that any deviation 
from this behaviour sh ould be orders of magnitude smaller than the form expressed here 
jstewart & Lyth . 1993 1- In the following we focus on the study of the tiny inhomogeneities 
produced by quantum fluctuations of the inflationary field. This aspect is crucial in under- 
standing the origin of the observed structure in the universe. 

2.4.1 Inflationary field power spectrum 

The mean amplitude of matter or curvature perturbations is inferred from their power spec- 
trum. This is constructed through the quantizatio n of real perturbation fields, as prescribed 



by quantum field theory IIBirrell & Daviesl 1 198411 . Following a perturbative expansion, we 
split the scalar field as 

<f>{t,x) = M*)+ W,x). (2.98) 

Such an expansion separates the general Klein-Gordon equation (12.651) into its homogeneous 
part (12.661) and the perturbation equation 

Scj) + 3H5(f) - + m 2 6 J<f> = 0, (2.99) 

where the effective mass of the field fluctuation is defined as 

m 2 si = ml + V& L , (2.100) 

and V/vl = V—mi(j) 2 /2 represents the nonlinear part of the potential. Note that here we have 
neglected the perturbations of the metric which enter the Klein-Gordon equation through the 
operator , v v . In this case, Eq. ( 12.991) is a valid approximation because we consider </> to be 
subject to the slow-roll conditions. 
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In order to compute the power spectrum of the field perturbations, we need to consider 
the free field 54>t, or the fi eld in the so-called interaction picture of quantum field theory 
llPeskin & Schroedeii 1 199511 . This field is the solution to the Klein-Gordon equation (12.991) 
in the absence of the nonlinear term, i.e. 



^ I + 3HS(f> I -^^- + mlS ( f )I 



0. 



(2.101) 



In the quantum framework we Fourier decompose this field as 



(x,t) 



dk 3 



e lk '*a k 5(f) Ik (t) 



-ikx t 



(2.102) 



where * denotes complex conjugation and f the hermitian adjoint operator. a k and a\ are 
operators satisfying the usual canonical commutation relations, 



[a k , a)} = 5 {3) (k - p), [a k , a p ] = [ 



a ki a p\ 



0. 



(2.103) 



and 5^(k) is the three-dimensional Dirac delta function. Thus the field fluctuations in 
Fourier space are solutions of the linear equation 



V + 3H5<P Ik + — + m 



0. 



(2.104) 



For simplicity, we only consider de Sitter inflation, where H is constant, and the scale factor 
is a = —l/(Hrj). As mentioned above, this is a good approximation in slow-roll inflation. 
In terms of conformal time, the previous equation can be written as 



50" +37r 1 ^' 



Ik 



k 



0. 



(2.105) 



The solution to this equation involves the set of Bessel complex functions. After proper nor 
malisation, i.e., taking the Bunch-Davies vacuum for the de Sitter spacetime [|Bunch & Davies . 
1988D, one finds 



where H^ 1 ^ is the Hankel function of the first kind and of order 



(2.106) 



v 



The star indicates that we are evaluating H just after Hubble horizon exi 



& i.e., 



(2.107) 



when k 



2 In the perturbed KG equation the mass term is negligible because of the slow-roll condition V" /V <C V , which 
is equivalent to m$ <C H. Well before horizon exit the harmonic flat spacetime equation is recovered in the solution 
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The two-point function of the field fluctuations reads 

(S^M^iM) = (2vr) 3 5( 3 )(k 1 + k 2 ) 1 ^(|A; 1 r / |) 3 |H( 1) (IMI)r. (2-108) 

where the angled brackets indicate the expectation value, in this case, of two modes of the 
perturbation field, evaluated at time 77. In terms of the two-point function, the power spectrum 
P^ is defined as 

(SMvWkM) = C27r) 3 5^ (ki + k 2 ) P^, h), (2.109) 
and the dimensionless power spectrum in this case is 

V^k) = f- 2 P,(v,k) = ^(\k V \f |H«(|^|)| 2 . (2.110) 

By expanding the Hankel function in Eq. (12.1081) on large scales, i.e. for \krj\ <C 1, one 
obtains 

n(ri, k) = 2- 2 - F( %^ Hi (\k v \) 2 " ~ (1 + (\k V \r , (2-1H) 



where 



3 



ws 2-"~5«- <2 ' 112) 

This is negligible in the massless case which corresponds to a de Sitter inflationary phase. 
Note that, in this case, the approximation of a linearised potential is guaranteed by the slow- 
roll approximation. Higher order terms in the Klein-Gordon equation are suppressed by 
powers of the slow-roll parameters. In this way, interaction between Fourier modes are 
absent, i.e., the vacuum fluctuations of different Fourier modes are decoupled and the field 
fluctuations are Gaussian. 



2.4.2 Observables 

In the standard picture of inflation, the perturbations of the inflaton field are stretched out 
of the horizon and subsequently transferred into curvature perturbations which survive after 
the universe has reheated. Observationally, the temperature inhomogeneities in the CMB are 
related to the mean am plitude of the curvature perturbations through the Sachs- Wolfe effect 



||Sachs & WolfeL Il967|1 . Thus, the power spectrum of curvature perturbations is observed to 



Eq. J2.106t . On the other hand, there is no need to compute the correlation at times well after horizon exit. The relation 
J2.44b and the fact that 1Z is conserved well outside the horizon indicate that the required field power spectrum can be 
evaluated a few Hubble times after horizon crossing. 
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be IIKomatsu et all 1200811 



V 



■R 



(2.45±g:8g) x io- 



95% CL). 



(2.113) 



at a pivot scale k = 7.5a H « 0.001 /i Mpc -1 . The power spectrum is also probed at 
other scales in the CMB with various filter functions and also with the power spectrum of 
galaxies and clusters. The rate of change of the observed value of V with k is parametrised 
by the spectral index n s . This is defined by 



n. — 1 = 



dlogPrc 
dlog/c 



(2.114) 



where 1 is subtracted by convention due to the fact that the matt er density power spec trum 



has the form P p oc k Tls . Observationally, the outp ut from WMAP IIKomatsu et al 



distance measurements from typ e la supernovae IIRiess et al. L|2007|; 



the baryon acoustic oscillations [|Percival et all l2007|l constrain the spectral index to be 



Astier et al 



2008], the 



200611 and 



-0.256 < l-n s < 0.025 (95% CL.). 



(2.115) 



for 0.001 Mpc" 1 < k < 0.1 Mpc" 1 . 

In single-field inflation, the power spectrum of curvature perturbations can be calculated 
from V<f, and Eq. (12.441) . which for a single scalar field can be written as 



K 



H 
— ( 





(2.116) 



This relation, used at linear order in Eqs. (12.1091) and (12.1 101) . shows that 



H 



2 



VlZ= {j) V * ~ ^ + 0(U})) (2 ' 117) 

In terms of the slow-roll parameters, using Eqs. (12.951) and (12.961 ) and evaluating the previous 
expression at horizon crossing, we have 

_ 2 8 K 3 _8K 

Note that the evaluation of the power spectrum at horizon exit is justified by the fact that 1Z is 
constant on super-horizon scales. In this regime we can define the root-mean-square (RMS) 
value of 1Z (or Q as 

H H 2 

■ft-RMs = Crms = -t<$0rms ~ — r ■ (2.119) 

(p 2n(p 



This quantity will play an important role throughout this thesis. 
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We can also write the spectral index as llStewart & Lythl 1 199311 



1 + 2?7sr - 6e SR , 



(2.120) 



to lowest order in slow-roll exp ansion. When re, = 1, the power spectrum is called scale 



invariant or Harrison-Zeldovich [|Harrisonl . 



1970 : IZeldovichLll972h . When n s ^ 1, the spec- 



trum is described as tilted and, for re. f > 1, i t is called 'blue' because the power is enhanced at 
small wavelengths IIMollerach et all 1 199411 . With the precision reached by the latest probes 
of cosmological structure, it has been possible to constrain the field param eters and discard 
some models of inflation a lterna tive to the simplest picture presented above l Alabidi & Lyt' 



2006al : lAlabidi & Lidsey . 



h. 



200811 . Current observations of the CMB and large-scale struc- 



ture, however, are compatible with the predictions of many other models of inflation. It is 
therefore crucial to study additional observables which provide further insight into the char- 
acteristics of the early universe. At the level of scalar perturbations, the most convenient 
observables for discriminating between these models are the tensorial perturbations mean 
amplitude and spectral index, the running of the spectral index and the non-Gaussianity of 
perturbations. In this thesis we focus on the effects of the latter. 



2.4.3 The 5N formalism 

We now present a formalism to account for the contributions of multiple fields to the curva- 
ture perturbation at all orders. An important feature of the perturbed curvature on superhori- 
zon scales is that we can calculate its magnitude by considering the change in the number of 
e-folds of expansion of the relevant patch of universe with respect to a uniform background 
expansion. This in turn allows us to compute the amplitude of the curvature fluctuations 
from the matter fluctuations. The idea behind this technique is to consider C as a pertur- 



bation in the local expansion [ Starob inskvL 



1996 



Sasaki & Tanakal . 



1998 



1985 



Lythetalll2005h . i.e. 



: balopek" 



& Bond . 1 19901 : ISasaki & Stewart . 



(2.121) 



where 5N is the perturbed expansion of the uniform-density hypersurfaces with respect to 
spatially flat hypersurfaces. 

We now describe the elements of the above formalism. The number of e-folds of expan- 
sion between two moments in proper time T\ and r 2 is given in the homogeneous background 
by 



N 



-0„ dr 



t-2 



-9(t) dt 



Hdt, 



(2.122) 



where 9 refers to the homogeneous expansion, that is, 9 considered to lowest order in 
Eq. (I2.29I ). In the perturbed metric, the number of e-foldings along an integral curve of 
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the 4-velocity, i.e., along a comoving worldline between T\ and r 2 , is 

M = Sr 2 ^ 6dT = it 2 + " ^ ^ = it + ^ ^ (2 ' 123) 

The last equality holds on superhorizon scales. It is clear that the difference J\f — N will 
provide the change in the curvature perturbation from an initial hypersurface at time ti and 
a final one at time t 2 , i.e., 

5N = Af - N = Aip. (2.124) 

In particular, we can choose to integrate the expansion starting from an initial uniform- 
curvature hypersurface at time t\, that is, ip(ti) = 0. Then we can find the amplitude of 
the curvature perturbation at the later time t 2 by choosing a trajectory with an endpoint t 2 
fixed on a comoving or uniform-density hypersurface. Denoting the difference between the 
background and the perturbed expansion as 56, we can write 

SN= I" 59 {t) dr) = i> 2 . (2.125) 
Jti 

In particular, when we consider the endpoint embedded in a uniform-density hypersurface, 
then Eq. (12.1211) is recovered. 

On large scales, where spatial gradients can be neglected, the local physical quantities 
like density an d expansion rate obey the same evolution e quations as in a homogeneous 
FRW universe [ Wands et al. . 200ol : Sasaki & Stewart! . 1996 1. By using homogeneous FRW 



universes to describe the evolution of local patches, we can evaluate the perturbed expansion 
in different parts of our universe with particular initial values for the fields during inflation. 
This is known as the 'separate universe' approach and means that, when we neglect the 
decaying mode for the field perturbations on superhorizon scales, we can consider the local 
integrated expansion as a function of the local field values on the initial hypersurface. In 
particular, one can expand Eq. (12.1211) as 

C = SNfafa)) = V (2-126) 

~ 0(pi 

where the initial time t\ again corresponds to some initial spatially-flat hypersurface. We can 
use this formula to construct the curvature power spectrum in multi-field inflation: 



This formalism can be extended to establish the equivalence between nonlinear matter and 
metric perturbations. This is done by first assuming Eq. (12.1211) as the definition of the 
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curvature perturbation and then using Eq. (12.441) to integrate (. This gives IILyth et all 1200511 
C = -A — V ln f — J • ( 2 - 128 ) 



3(w + l) \po 

This nonlinear curvature perturbation can be written as a function of the initial field fluctua- 
tions, evaluated again at an initial flat hypersurface labelled by the time t\. We use a Taylor 
expansion 

< = Ea 5 ^(«>)+5flps'*'W*w + - (2 - 129) 

i J 

where the leading order term coincides with the expansion (I2.126I) . This last expansion 
greatly simplifies the derivation of the higher-order curvature correlations from the scalar 
field bi spectrum. 



2.5 Non-Gaussianity from isocurvature fields 

The perturbed energy conservation equations show the conditions under which the curvature 
perturbation £ may vary over time in a regime in which the perturbation modes lie well out- 
side the horizon. Specifically, Eq. (I2.87I) shows that the evolution of ( is directly related to 
the presence of an entropy perturbation 5s. This quantity can be the intrinsic non-adiabatic 
pressure of a single field, or the difference in the density perturbations of any two fields 
which contribute to the curvature perturbation. This effect motivates the study of models of 
inflatio n in which the observed inhomogeneities are the result of non-adiabatic field fluctu 



ations JMoller ach 



199 



0: 



Linde & MukhanovL 1 19971 : Enqvist & SlothL 12002 : iLyth & Wands . 



20021 : iMoroi & TakahashiL 1200 111 . Here we are interested specifically in models in which a 
highly nonlinear £ can be generated from the aforementioned entropy perturbation. Out of 
the multiple stages in which this field may have influenced (, we focus on the ca se of an aux- 
iliary scalar field during inflation generally referred as the isocurvature field x [|Zaldarriaga . 



20041 : Enqvist & Nurmi . 



2005 



Enqvist et all l2005bh 



A first approximation to nonlinear ( involves the first and second order perturbations in 
real space: 



C(t,x) = Ci(t,x) + -C 2 (t,x). 



(2.130) 



The second order perturbati on is conventionally writte n in terms of the first order perturba- 
tion and the parameter / NL [|Komatsu & Spergell . l200l|l . This gives 



C(t,x) = Ci(t,x) - -/ NL (d 2 (t,x) - (Cftt,x))) , 



(2.131) 
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which in Fourier space is written as a convolution 



C(k) = Ci(k)-§/NL([Ci*Ci](k)-<C?(*)». 



(2.132) 



Note that Eq. (12.1321) shows explicitly the superposition of modes that characterise non- 
Gaussian statistics. 

Statistically, non-Gaussianity refers to the non-vanishing higher order moments of the 
quantity in question. In quantum mechanics this corresponds to the n-point correlation func- 
tions with n > 3. To lowest order, the bispectrum B^(ki, k 2 , fc 3 ) is defined by the expectation 
value of the product of three copies of the curvature field: 

(C(kOC(k 2 )C(k 3 )) = (2n) 3 B c (k 1 ,k 2 ,k 3 )5^(k 1 ,k 2 M). (2.133) 
The amplitude of /nl is given in terms of the bispectrum by substituting Eq. (12.1321) in this: 

(2.134) 



6 U t kf gc 



In this section we present a method to compute the nonlinear ( by determining nonlinear 
solutions to the Klein-Gordon equation of the field fluctuation. Then with the aid of the SN 
formalism we will construct the three-point correlator of (. Special attention will be paid to 
the cases in which a large /nl can be obtained. 

2.5.1 Two-field inflation 

Here we consider a spacetime inflating by the action of a potential which depends on two 
minimally coupled fields. In this case, in addition to the canonical inflationary field, we 
consider a second field \ described by the action 



dx 4 



-g 



(2.135) 



where W(x) plays the role of a nonlinear potential. The joint Lagrangian density of the two 
scalar fields in this model is given by Eq. (12.521) . The model in question demands that there 
is no contribution of x to the background matter content. This is guaranteed by the following 
condition on the potential of the two-field Lagrangian defined in Eq. (12.521) : 



U{<f>,x)^V{<f>) 



W( X ) + 7 :m'x 2 ^V( ( f ) ). 



(2.136) 



To analyse the dynamics of this auxiliary field, let us consider a background homogeneous 
part and a coordinate-dependent field fluctuation, 



= Xo(t) + £x(*,x). 



(2.137) 
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The Klein-Gordon equation derived by varying the action (12.1351) with respect to X yields 
for the background field 

Xo + 3H X o + m 2 xX o + W'( X o) = 0, (2.138) 
and for the field fluctuation 

S X + SH5 X + m 2 Sx S X + £ T^W** = °" (2 " 139) 

n=3 ^ 

Here a superscript (n) denotes the n-th derivative and rrig is the effective mass for the field 
fluctuation defined as 



m 



ml + W". (2.140) 



Following the same steps as with the inflaton power spectrum, we recover a solution similar 
to Eq. (l2~T06l) . Specifically, 

V X (V, k) = ^{\Hf |HW(N)| = 2" 2 " r ^ ~ ^ ^ (l^|) 2a , (2.141) 
with a and /i defined through the equation 

/ Q TTL 2 

«-V^ 2 "i"^ >0 - (2 - 142) 

Because of this, the spectrum of X is blue. The power spectrum V x vanishes on large scales 
for large values of a, i.e. when mj x > H*, so we only consider the case 

m Sx <H if . (2.143) 

Indeed, we assume m x <C H*. As shown below, the interesting values of W" are those 
which are large enough to generate a large nonlinear coupling, but sufficiently small to satisfy 
condition (12.1431) . which also implies 

W"<Hl (2.144) 

For example, for 

W = ^X 3 + ^X\ (2.145) 
this condition turns into conditions on M and A: 

M<^, A 1/2 < — . (2.146) 
Xo Xo 
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Typically we have a nonzero expectation value Xo- Working in a perturbative expansion 
then requires that the quantum fluctuations 8x do not exceed this expectation value, i.e., we 
require \ > H*. This imposes additional constraints on the parameters M and A in the 
specific model of Eq. (12.1451) . as described below. 



2.5.2 Non-Gaussianity and nonlinear evolution 



The expectation value of the product of three fields was first computed by lMaldacenal (1200311 
for the case of single-field inflation, including slow-roll scalar field and gravitational interac- 
tions. In our case, since the field gives a negligible contribution to the energy density of the 
Universe, we will neglect its coupling to gravity. 

The three-point correlation function of the field perturbations Sx can be computed using 
the expression given by Maldacena: 



(Sx(t, xi)6x{t, x 2 )flx(*, x 3 )) 



dt'dHjit'), S X i(t, xi)5xi(t, *2)Sxi(t, x 3 )]>, 

(2.147) 



where Hj is the interaction Hamiltonian written in terms of the field perturbation in the 
interaction picture, i.e. in our case 



J n=3 



M/( n ) 
dx 3 a(t') 3 — -5x1- 



(2.148) 



Note that on the left-hand side of Eq. (12.1471) the expectation value is taken with respect to 
the vacuum of the interacting theory, while on the right-hand side it is taken with respect to 
the vacuum of the free theory. A generalisation o f this expression t o higher order correlators, 
including loo p corrections, has been provided by lweinberg [ 2005 1. 

Recently, iMussoJ |2006|] showed that the general expression for the correlators of m field 
fluctuations can also be derived by solving perturbatively the field equation of motion and 
using the expression for the two-point function of the free field fluctuation. Here we make 
use of this formalism to show that Eq. (12.1471) can be derived by solving Eq. ( 12.1391 ) pertur- 
batively. 

As mentioned before, the evolution equation for 5xi is given in identical form to the 
inflaton case, providing <p — ► X m Eq. (12.1051) . On the other hand, the evolution equation for 
the full nonlinear Sx can be given perturbatively providing | (5x — &Xi)/&x\ ^ 1- Rewriting 
Eq. (12.1391) and re-expressing the nonlinear term in terms of 5xi, we obtain 



5 X + SH5 X 



n=3 



n 



n-l 



(2.149) 
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This equation can be then solved by Green's method and its solution is 



n=3 V 1 >' ^ J 



where G R (t, x; t', y') is the retarded Green's function 



G R (t, x; y) = iQ(t - t') [S X i(t, x), S X i(f, y)] 



where 



G(x) 



5{z) dz 



(2.150) 



(2.151) 



(2.152) 



is the Heaviside step function. By using the fact that the free-field perturbation is Gaussian, 
i.e. 

(8 X i(xi, t)6 X i{x2, t)5xi(*3, *)) = 0, (2.153) 
one can rewrite the three-point correlation function of 5x making use of Eq. (12.1501) : 

(5x(t,xi)5 X (t,X2)5x(t,x S )) = - / ^ / dt ' a (t') 3 W^ x 

n=3 ^ U ''J 

[6xi(t,Xi),5 X i(t',y)} (5xr 1 (t',y)5x/(t,x 2 )<5x/(t,x 3 )) + {Perms.} (2.154) 



Using this equation and the result jMussoLfeOOql . 



n [<Jxj(*,Xi),<Jxr(* ',y)](8Xi y)6xi(t, x 2 )5x/(t, x 3 )) + {perms} 
-( [5xW> y)> S Xi(t, x-i)6xi(t, ^)6xi(t, xa)]>, 



(2.155) 



one can derive Maldacena's formula (12.1471) for the x-Hamiltonian (12.1481) after some ma- 
nipulation. This shows that one can obtain Maldacena's formula either from considering 
perturbations in the action or from the perturbative expansion of the equations of motion. In 
the case of slow-roll inflation, this equivalence was shown by Seery et al. | 2008| . Equation 
(12.1471) can be generalised to higher correlation functions: 

{S X (t, *i)6x(t, x 2 ) . . . 5 X (t, x m )> = (2.156) 



i / dt'{[H I (t'),5xi(t,x 1 )5 X i(t,yi2)---Sxi(t,^ 

J — OO 
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2.5.3 Field bispectrum 



Here we are interested in the three-point function of the scalar field, also called the field 
bispectrum F{k\, k 2 , k 3 ), defined as 



(2.157) 



For simplicity, we assume that the nonlinear potential of the scalar field is dominated by the 
cubic interaction W" and that this is approximately constant. This reduces the expansion in 
Eq. (12.1481) to the first term only. Substituting this in Eq. (12.1471) and using Eqs. (12.1501) and 
(12.1551) . the field bispectrum becomes 

F(r);k u k 2 ,k 3 ) = (2.158) 

W"V /2 ^-Re {-i^n^OMI) £ drfy/Wl[llM{\ktf\ 



During de Sitter inflation, the ratio k/aH can be written as the product \krj\ and we will use 
this to parametrise the various stages of evolution of the bispectrum. 

It is instructive to first evaluate the bispectrum (12.1581) during inflation (r) < ?7 rcll ) for the 
case of a massless field fluctuation (i.e., ms x = a = 0, u = 3/2" ) . In this case, the integra l 
in Eq. (12.1581) can be evaluated analytically I Bernardeau & Uzan . 20031 : Zaldarriaga . 2004 1 
using the expressions for the Hankel functions, 



H$(z) 



71 Z 6 ' 2 



1Z 



y3/2 ■ 



(2.159) 



One find 



F(r);k 1 ,k 2 ,k 3 



W"'Hl 

4EU 3 



^k 3 t +k t J2k i k j + ~ (^ + j + ]n(\k t r]\)J J2 k i 

i<j ' ' i 



(2.160) 



where k t = fcj. The integral in Eq. (12.1581) has been evaluated for three different stages: 
when the modes are well inside the Hubble radius, around the Hubble radius and outside 
the Hubble radius. The first stage does not give any contribution to the integral because the 
Hankel functions oscillate rapidly for \kiT)\ 3> 1 and the fields can be taken as free in the 
asymptotic past. At late time, the integral is dominated by the modes that are around the 
Hubble scale or larger, \kr)\ < 1. In particular, we will show below that the time-dependent 
term in Eq. (12.1601) . with the typical local momentum dependence, essentially comes from 
the nonlinear and classical super-Hubble evolution of the field perturbation. The finite part, 



3 Our result coincides with the one found in lZaldarriagal l2004ll . modulo the overall sign and the factor 1/3 inside the 
parentheses. 
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with the non-trivial momentum dependence, comes from integrating over times correspond- 
ing to Hubble-crossing. 

Let us now consider the case with ms x ^ 0, and 1/2 < fi < 3/2. The integral inside 
Eq. (12.1581) cannot be integrated analytically in this case. One can only evaluate it over the 
period when all modes are well outside the Hubble radius because, for small arguments, 
Hankel functions can be written in terms of Bessel functions (which for real arguments are 
real). Decomposing Eq. (12.1581) in this fashion, using the Bessel function expansion for small 
arguments and retaining the dominant real part, then yields 

F{ V , k u k 2 , k 3 ) =F( V *, k u k 2 , k 3 ) + — — ± (2.161) 

12 [iir z a 



2/i 



,3-30 



2/i 



(v/v*r a 



n*(*i 



, 2n 



The first term on the right-hand side is the non-Gaussianity at Hubble exit and it can only 
be computed numerically. The second term represents the non-Gaussianity generated at late 
times during inflation, when the large-scale local term dominates over the finite Hubble- 
crossing term. 

Now we will show that the large-scale local contribution to the non-Gaussianity of Eqs. 
(12.1601) and (12.1611) can be derived by solving the equation of motion of the field perturbation 
on large scales. The nonlinear evolution of the field fluctuation is derived by taking the large- 
scale limit \krj\ —> in Eq. (12.1051) for the free-field 5xi and Eq. (12.1491 ) for the nonlinear 
5x- This leads to the equations 

5x Ik + 3H5 Xnc + m] x 5xi* = 0, (2.162) 
$Xk + 3H5x k + m 2 Sx 5xk = S k , (2. 163) 

where the source on the left-hand side is given in terms of the linear solution, 

Sk = -^(5 X i*5 X i) k , (2.164) 

and where * denotes the convolution operation. The growing solution of the homogeneous 
Eq. (1271621) is 

SXiM-SxM(\krj\)i^. (2.165) 



Only in the massless limit is this constant. One can find the solution of the inhomogeneous 
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equation by using the method of variation of parameters, which yields 

(1^1)"^ [" dv'aWQkviyi-^r,') 



s_ 1 

SxM =^Xk(^)(|^|) 5 fl - 



(2.166) 



Using the source (12.1641) with Eq. (12.1651) and integrating over conformal time, one ob- 
tains 

SxM = SxMdWT + ^f(v/v*)[Sxi(v)*S X i(v)}u , (2.167) 
where 

f | (l-x 3 )+\n(x) for a = 3/2, 

fix) = { i / „ „ , \ , , , (2-168) 
1 ^172) (l - f^ 3 ~ 3a ) ~ ** 1/2 < A* < 3/2. 

We can now use Eq. (12.1671) to compute the bispectrum from its definition (12.1571) . In terms 
of the power spectrum of the field perturbation, this is 

F{ V ; h, k 2 , k 3 ) = F{ V , h, k 2 , k 3 ) + ^f{v/v*) Yl P ^ k ^ P ^ < 2 - 169 ) 

* i<3 

which correctly reproduces the large-scale contribution to the non-Gaussianity in Eqs. (12.1601) 
and (12.1611) . Note that, in the limit fi — > 3/2 and a — ► 0, the large-scale expression for a 
massive field converges to the massless case, as can be checked by taking this limit in the 
lower expression on the right-hand side of Eq. (12.1681) and using 

2/1 i , a , /n/ 2 



l + - + 0(a 2 ), (2.170) 
3(/i-i) 3 

x~ a = exp [-a ln(z)] -> 1 - a ln(z) + C(a 2 ). (2.171) 

In summary, at late times the non-Gaussianity of the field is dominated by the nonlinear 
evolution on large scales and thus the bispectrum of the field perturbation is of the local 
form, i.e. proportional to the product of two power spectra (See Sec. II .31 for the definition 
of triangulation of the bispectrum). In the massless case, in which [i = 3/2, the power 
spectrum of the field fluctuation is constant while the bispectrum grows as ln(a). In the 
massive case, however, the power spectrum decays as a~ 2a and the bispectrum decays as 
a~ 3a , thus 'growing' as a a with respect to the product of two power spectra. This relative 
growth is important for the non-Gaussianity in the curvaton mechanism, as can be seen from 
Eq. (12.1341) . The curvaton case will be discussed below in more detail. 

Before making contact with curvature perturbations, let us extend the solution (12.1671) 
and consider the evolution of the non-adiabatic perturbations in a radiation-dominated era. 
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This applies once inflation has ceased and the inflaton field has been thermalised, but before 
the isocurvature perturbation is converted into an adiabatic one. 

During the radiation dominated era, H = (2t) _1 and the field perturbation evolves on 
large scales according to 



6Xlk + 2t 6Xl]c + m Sx 5xik = °' 
3 • 



(2.172) 
(2.173) 



The growing mode of the h omogeneous equation is a Bessel function of the first kind J 
iLanglois & VemizziEool 



7r J\/i{rn 5x t) 
2 5 /4r(3/4) (msxt) 1 / 4 



(2.174) 



For m$ x t <C 1, well before decay, the growing mode is constant. Indeed, by solving 
Eq. (12.1731) at lowest order in mg x t, we finqj 



rch j 



7T Ji/4(m Sx t) 



2 5 /4r(3/4) (m Sx ty/ 4 



10mj x 



{m Sx t) 2 [5xi{t)*5xi{t)} k . (2.175) 



This result shows that the isocurvature fluctuation continues its nonlinear evolution through- 
out the radiation era. We will use this result in the context of a curvaton field to account for 
the consequences of considering a nonlinear source in the bispectrum of the field and then 
compute the corresponding non-Gaussian parameter /nl- 



2.5.4 The curvaton 

The cu rvaton is an alternative inflationary mechanism to generate the matter density fluctu 
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ation s ILM ollerach 



1990 



Moroi & Takahashi. 



iL booih 



Lind e & MukhanovL 1 19971 : Enqvist & SlothL 120021 : iLyth & Wands , 



200 111 . This is achieved without appealing to the perturbations in 



the original inflaton field. Instead, an auxiliary 'curvaton' field, subdominant during infla- 
tion, generates isocurvaure fluctuations which transform into adiabatic ones after the infla- 
tionary phase, during the decaying oscillations of the curvaton field. During inflation, the 
isocurvature field presents a negligible contribution to the energy density. After inflation the 
field still plays no significant role in the background evolution as long as its mass m x is neg- 
ligible compared to the Hubble parameter. However, once m x fa H 2 , the curvaton field starts 
oscillating at the bottom of its potential. At this stage the potential can be approximated as 
quadratic. The energy density of the field decays during the oscillations like a non-relativistic 
component (p x oc 1/a 3 ). The curvaton then contributes significantly to the energy density of 
the universe and the curvaton fluctuations are transformed into adiabatic matter fluctuations 



4 This result is in agreement with lEnqvist & Nurmil l2005ll where the computation considered a general nonlinear poten- 
tial up to order O ((ms x t/2) 8 ). 
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(for the simplest version of this mechanism see lBartolo & Liddld Q2002|]). 

We can write the number of e-folds of expansion in terms of the value of the field x- 



N(t dec ,t in ) = Iln(-^ 



(2.176) 



where 



2 2 



PXdo 



2 2 



(2.177) 



are the energy densities of the curvaton at the onset of the oscillations (on a fiat slicing) and 
at the moment of decay (on a uniform density slicing), respectively. 

The non-Gaussianities generated by the curvaton mechanism have been studied in sev- 
eral papers. In the following we consider and combine all the possible effects, including the 
intrinsic non-Gaussianity of the curvaton field fluctuation and the nonlinear relation between 
the curvature perturbation £ and the curvaton fluctuation. The intrinsic non-Gaussianity of 
the curvaton can be generated inside and outside the Hubble radius due to its nonlinear po- 
tential. In particular, as discus sed later in this section, we take into account the nonlinear 



evolution du ring both inflation JIBemard eau & Uzan . 



ation epoch IIEnqvist & Nurmil 12005 



Lvt 



h, 



2004]. 



2003 



ZaldarriagaL l2004f| and the radi 



In order to compute ( and its n-point functions, one can follow two equivalent procedures: 

1 . Expand the perturbed number of e-folds 5N on an initial flat slice at the onset of the 
oscillations (t = t in ), in terms of the field fluctuations Sx(t- m ), and then use Eqs. (12.1691) 
and (12.1751) to introduce the three-point correlators of the field fluctuations. 

2. Expand the perturbed number of e-folds 5N on an initial flat slice at Hubble cross- 
ing (t = £*) in terms of the field fluctuations Sx(tm) and then take into account the 
nonlinear relation between the field fluctuation at t = t in and the one at t — t* using 
Eqs. (I2T671) and (I2T751) . 

We will follow the latter procedure, which has been used in the rest of the literature on the 
curvaton. We write, as in Eq. (12.1291) 



C = N x J X (t* 



In 

2 ' 



x*x 



.Sx 2 (U) 



(2.178) 



where N is given by Eq. (12.1761) . 

In general, as we have seen in the previous section, Xm is a nonlinear function of the field 
value at Hubble exit, and we parameterise this dependence by the function g(x{t*)) an d its 
derivatives: 



9 =Xo{U. 



(n) 



(2.179) 
(2.180) 



n=l 
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Using -t^- = g'-^j and Eq. (12.1771) . one can differentiate N in Eq. (12.1761) to obtain 

2^ 



3 ^ 



(2.181) 



where the prime denotes here a derivative with respect to x* and 



c = i- dlnpx 



3p> 



(2.182) 



dec 



where p is the unperturbed total energy density. Here we have taken a uniform p dcc assuming 
that the radiation is unperturbed. Also, to arrive at the last equality we used (p Xdcc / PxiJ 1 ^ 3 = 

[(Pdcc-Px dcc )/(Pin-Pxin)] 1/4 - 

To compute the power spectrum we neglect the evolution during the radiation dominated 
era, and from Eq. (12.1671) we obtain g'/g = (\kr]\) a , which yields 



V ( (t,k) = ^C 2 V x (t,k). 



(2.183) 



The spectral index is given by iLyth & Wands! [ 1 200211 as 



H 

n s -l = 2— = 2a, 



where Eq. (12.1421) implies 

,2 



(2.184) 
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a 



3H? 



< 1. 



(2.185) 



If a is not too small, th e spectrum can be extremely blue and this is ruled out by observations 
Komatsu et al. . 2008 1. Differentiating N in Eq. (12.1761) once more yields 



N 



,x*x* 
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2C 
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2-C 



(2.186) 



If we neglect the non-Gaussianity of the field fluctuations at Hubble crossing, which are 
subdominant with respect to the ones accumulated during the super-Hubble evolution, and 
the definition (12.1341) gives 



NL 



5 

4C 



1 



99 
9' 2 
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3 



5C 

IT 



(2.187) 



We have arrived to a well known result o btained without assuming a dominant contribu- 
tion of W'" in the perturbation equations [|Bartolo et all 120041 : ILyth & RodriguezL 12005b 



Sasaki et al 



20061 



By comparing the above equation with the nonlinear evolution given by Eqs. (12.1671) and 
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(12.1751) . and stopping the nonlinear evolution of perturbations when the field starts oscillating 
at t ~ l/mg x , the non-Gaussianity in the curvature perturbation becomes 



f _ 5 



XQin^) 1 



(2.188) 

3 o 



where the function / has been defined in Eq. (12.1681) . At late times, / can be approximated 
by 

, , f -AN for aAN < 1, 



where we have used Eq. (12.1851) and AN = Nt^—N* ~ 60 is the number of e-folds between 
Hubble crossing and the end of inflation. 

Note that our result, Eq. (12.1881) . can also be obtained by computing the field bispectrum 
from the isocurvature field evolution during radiation domination Eq. (12.1751) . In this context, 
the field bispectrum is given by 

F(t; h, k 2 , k 3 ) = F(t rch ; k h k 2 , k 3 ) - — V P( V , ki)P(r), kj), (2.190) 

which clearly evolves with time before the decay of the curvaton. This non-Gaussianity is 
transferred to the curvature correlation by expanding the perturbed number of e-folds 8N on 
an initial flat slice taken at the onset of the oscillations (t = t- m ). Then we can write £ = 5N 
in terms of the field fluctuations Sx(t- m ) as in Eq. (12.1291) . With the aid of Eqs. (12.1691) and 
( 12.1901 ), we then replace the three-point correlators of the field fluctuations in the curvature 
perturbation bispectrum. Our result in Eq. (12.1881) is thus recovered. 



2.6 Model discrimination through observations 

In its simplest version, the curvaton model proposes an isocurvature field whose quantum 
fluctuations reproduce the spectrum of curvature fluctuations we observe in the CMB. Fixing 
the required perturbation amplitude and spectral index imposes important restrictions on the 
possible values of X- Additionally, the non-Gaussianity of the curvaton could constrain the 
parameter space further. 

As an example, let us consider a curvaton with bare mass m x and nonlinear poten- 

tial W = |f x 3 , such that mj x ~ W" = xo *M and W" = M. Note that the negligible mass 
in this case allows for a scale-invariant power spectrum in the field fluctuations and conse- 
quently the curvature perturbations. For simplicity we take xo osc — Xo *> which is consistent 
with neglecting the nonlinear potential in Eq. (12.1381) ). since this would only contribute a 
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term of order 0(a 2 ). When the nonlinear coupling is very small, Eq. (12.1881) reduces to 

/ NL = c ( 1 " i aAiV ) " 3 " y a ~ AiY_1 ' (2 - 191) 

and the intrinsic non-Gaussianity of the curvaton field gives a negligible contribution to the 
total non-Gaussianity in the curvature perturbation. However, when the nonlinear coupling 
is important, the nonlinear parameter is 

/nl = ^ f 1 - -exp(aAN)\ - - - — , AiV" 1 < a. (2.192) 

In this case the intrinsic non-Gaussianity of the curvaton can be the main source of non- 
Gaussianity in the curvature perturbation. For example, with x °^^ ~ 0.07 and C — 1, one 
finds /nl — —100, which is within reach of current and future experiments. However, in 
this case, if the curvaton is the only field responsible for the curvature perturbations, the 
spectral index of scalar fluctuations will be largely blue. This is in disagreement with current 
observations and is therefore excluded. 

Let us finally consider the special case in which the curvaton is not responsible for the 
linear fluctuations observed by CMB and large-scale structure probes. Using the same po- 
tential as in the example above, this case constrains the fraction r in Eq. (12.1821) to be small, 
as the contribution to the curvature power spectrum is controlled by this parameter (see 
Eq. (12.1831) ). On the other hand, a can take large values without violating constraints on 
the power spectrum, which is dominated by the inflaton perturbations. As for second order 
perturbations, /nl is still given by the formula (12.1921) and can be large due to the freedom 
in the mass and the small value of C. This happens at all scales and the non-Gaussianity is 
induced through the evolution of fluctuations on superhorizon scales. 

It is important to note that, in this case, the blue spectrum of the curvaton perturbations 
may dominate at small scales. If this happens, and if the perturbations have enough power 
on small scales, a significant amount of PBHs would be produced. In this special case, as 
in many other versions of inflation, an important constraint on the model comes from the 
probability of PBH formation as we will study in the following chapters. 



Chapter 3 

Statistics of non-Gaussian fluctuations 



3.1 Introduction 



In the inflationary paradigm, the prediction that the spectrum of fluctuations should exhibit 
Gaussian statistics has recently been challenged. This prediction follows from the fact that 
the curvature perturbation, which commonly refers to the comoving curvature perturbation 
defined in Eq. (12.431) . is treated as a free field during inflation, 



d 3 k 
(2^ 



ikx 



(3.1) 



where there is no coupling between the lZ(t, k) for different k. With this understanding, 
Eq. (13.11) means that 1Z does not interact with either itself or any other particle species in 
the universe. The real-space field 1Z(t, x) is obtained by summing an infinite number of 
independent, identically distributed, uncorrected oscillators. U nder these circumstan ces the 
Gaussianity of lZ(t,x) follows from the central limit theorem [|Bardeen et all Il986|l . given 
reasonable assumptions about the individual distributions of the lZ{t, k). The exact form of 
the distributions of the 1Z(t, k) is mostly irrelevant for the inflationary density perturbations. 

In conventional quantum field theory, all details of 1Z and its interactions are encoded in 
the n-point correlation functions of 71, written as (out|7£(£i, xi) • • -Tl{t n , x n )|in). Working 
in the Heisenberg picture, where the operators carry time dependence but the states 
{|in), |out)} do not, these functions express the amplitude for the early-time vacuum |in) to 
evolve into the late-time vacuum |out) in the presence of the fields 7 t{tj, x^). Given the re- 
point functions for all n at arbitrary x and t, one can determine lZ(t, x) [|Streater & Wightman . 
2000], at least in scattering theory. In the context of the inflationary density perturbations, 
these vacuum evolution amplitudes are not directly relevant. Instead, one is interested in 
the equal time expectation values (in\1Z(t, Xi) ■ • - TZ(t, x n )|in), which can be used to mea- 
sure gravitational particle creation out of the time-independent early vacuum |in) during 
inflation. These expectation value s are calculated us ing the s o-called 'closed - time-path for- 
malism', which was in trod uced by Schwingerl [|196l|l ; see also ICalzetta & Hul [|1987|l ; Ijordan 
II 1 9 8 611 : IPeWittl [1200311 and lHajicekl 111 97911 . In this formalism there is a doubling of degrees of 
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freed om, which is also manifest in finite te mperature calculations IILe Bellad 



2005 



200ft ISloth . 



2006 



2000; Rivers, 



Seeryl . l2008ll 



198811 ■ This method has recently been used [|WeinbergL 
to extend the computation of the correlation functions of 1Z to beyond tree-level. 

Knowledge of the expectation values of 1Z in the state |in) is sufficient to predict a large 
number of cosmological observ ables, including the power spectrum of the density perturba- 



tions generated during inflation [|Guth & Pi 



1982 



Hawkin 



point functions of the CMB temperature anisotropics IIHu & Sugiyamal 



Komatsu & Spergel 



ZaldarriagJ 



200 ll:lKogo & Komatsu . 



h, 



2006 



1 98 211, and the two - and three 



Okamoto & Hu, 



6: 



2002 



1995 



Hu 



2001 



Babich et al 



OoHifci 



2004 



Babich & 



2004; 



Babic 



20051 : iLiguori et al 



200> 



Cabella et al 



2006]. Because they are defined as expectation values in the quantum vacuum, these observ- 
ables all have the interpretation of ensemble averages, as will be discussed in more detail 
below. 

On the other hand, one sometimes needs to know the prob ability that fluctuat i ons of 



197 



4; 



some given magnitude occur in the curvature p erturbation 1Z [[Press & Schechterl . 
Bardeen et al. . 1986 : IPeacock & Heavens . 1990 1. This is not a question about ensemble 
averages, but about the probability measure on the ensemble itself. As a result, such in- 
formation cannot easily be obtained from inspection or simple manipulation of the n-point 
functions. 

For example, if we know by some a priori means that 1Z is free, then the argument given 
above, based on the central limit theorem, implies that at any position x, the probability 
density of fluctuations in 1Z of amplitude g must be 



P(f?) 



1 



2ixo 



exp 



2a< 



(3.2) 



where the variance in 1Z is 



a 2 = (JZ(t,x) 2 ) = / dlnkVnik) 



(3.3) 



The quantity Vn(k) is the dimensionless power spectrum, which is defined in terms of the 
two-point function of 1Z, calculated from the quantum field theory in-vacuum: 



(in|7^(t, k x )7e.(t, k 2 )|in> = (27r) 3 «5(k 1 



2tt 2 



(3.4) 



This is the only relevant observable, because it is a standard property of free fields that 
all other non-vanishing correlation functions can be expressed in terms of the two-point 
function (13.41) . and hence the power spectrum. Those expressions can be achieved through 
a gen eralisation of Wick's theorem using an equal-time normal-ordering ||Luo & Schramm . 
199311 . In practice, in order to give a precise meaning to (13.21) . it would be necessary to 
specify what it means for 1Z to develop fluctuations of amplitude g, and whether it is the 
fluctuations in the microphysical field 1Z or some smoothed field TZ which are measured. 
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These details affect the exact expression (13.31) for the variance of g. 

The average in Eq. (13.41) . denoted by (in| • • • |in), is the expectation value in the quantum 
in- vacuum. To relate this abstract expectation value to real- world measurement probabilities, 
one introduces a notional ensemble of possible universes, of which the presen t universe an d 
the density fluctuations that we observe are only one possible realisation (e.g. iLythl [2006D). 
However, for ergodic processes, we may freely trade ensemble averages for volume averages. 
The ergodicity of a system refers to that property of processes by which the average value of 
a process characteristic measured over time is the same as the average value measured over 
the ensemble. 

If we make the common supposition that the inflationary density perturbation is indeed 
ergodic, then we expect the volume average of the density fluctuation to behave like the 
ensemble average: the universe may contain regions where the fluctuation is atypical, but 
with high probability most regions contain fluctuations with root-mean-square amplitude 
close to a. Therefore the probability distribution on the ensemble, which is encoded in 
Eq. (13.41) . translates to a probability distribution on smoothed regions of a determined size 
within our own universe. 

In order to apply the above analysis, it is necessary to know in advance that 1Z is a free 
field. This knowledge allows us to use the central limit theorem to connect the correlation 
functions of 1Z with the probability distribution (13.21) . The situation in the real universe is not 
so simple. In particular, the assumption that during inflation 1Z behaves as a free field, and 
therefore that the oscillators 1Z(k) are uncorrected and independently distributed, is only 
approximately correct. In fact, 1Z is subject to self-interactions and interactions with the 
other constituents of the universe, which mix k-modes. Consequently, the oscillators lZ(k) 
acquire some phase correlation and are no longer independently distributed. In this situa- 
tion the central limit theorem gives only approximate information concerning the probability 
distribution of 1Z(x), and it is necessary to use a different method to connect the correlation 
functions of 1Z with its probability distribution function (PDF). 

In this chapter we give a new derivation of the PDF of the amplitude of fluctuations in 
1Z which directly connects P(p) and the correlation functions {lZ{k\) ■ ■ - lZ{k n )), without 
intermediate steps which invoke the central limit theorem or other statistical results. When 
the inflaton is treated as a free field, our method reproduces the familiar prediction (13.21) 
of Gaussian statistics. When the inflaton is 'not' treated as a free field, the very signifi- 
cant advantage of our technique is that it is possible to directly calculate the corrections to 
P(#). Specifically, the interactions of 1Z can be measured by the departure of the correlation 
functions from the form they would take if 1Z were free. Therefore, the first corrections to 
the free-field approximation are contained in the three-point function, which is exactly zero 
when there are no interactions. 

The three-point function for single-field, slow-roll inflation has been calculated by Maldacena 
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[2003|], whose result can be expressed in the form llSeery & LidseyL 12005 ah 



Vn 



(^(k 1 )^(k 2 )^(k 3 )) = 4vr 4 (2vr) 3 ( 5( £ k,) ^^A(h, k 2 , k 3 ), 



(3.5) 



where A is Maldacena's ^4-function divided by two llMaldacenal . I2003ll Fl VfC measures the 



amplitude of the spectrum when the kj crossed the ho r izon. (For ear l ier wo r k on the deriva 
tionof t he three-point-function, see 



11996H 



FalketaL 



lll993h : lGangui et all 



1 Jl994ll : 



Pvne & Carroll 



Acquaviva et all Q2003D.) This result has since been extended to cov er the non 



Gaussianity produced during s low-roll inflation with an arbitrary number of fi elds [ Mald acenal . 
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Lvth & Rodriguez . 
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2006Q, preheating IIEnqvist et all 



Lvth & Zaballa, 
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bEE 



200611 . models where the dominant non-Gaussianity is produced by a 



light scalar which is a spectator during inflation IIBoubekeur & Lythl 120061 : lAlabidi & Lythl 



2006bl : lLyth . 



2006 J1, and alternative models involving a small speed of sound for the inflaton 



pertur bation llSeerv & Lidsevl l2005bl : lAlishahiha et all 120041 : ICalcagnil 120051 : 



Arkani-Hamed et al 



2004 



CreminelliLl2003ll 



For single-field, slow-roll inflation, the self-interactions of 1Z are suppressed by powers 
of the slow-roll parameters. This means that the correction to Gaussian statistics is not large. 
In terms of the ^.-parametrised three-point function (13.51) . this is most commonly expressed 
by writing, in an equivalent form to £ in Eq. (12.1311) 



ft(t,x) =K 1 (t,x) - -/nl (Kl(t,x) - (nl(t,x))) 



(3.6) 



with 



/nl — — - ^ ^3 — 0(esR, rjsB.) 
giving the rel ative contribution of the non-Gaussian 
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piece in 1Z and 1Z\ being a Gaussian 
Verde et all I2OOO1I . (Note that there are differing 



20060, here we stick to that used by the WMAP 



team.) In models with mor e degrees of freedom, much larger non-Gaussianit i es are e xp ected , 
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Rigopoulos & Shellard, _ 
Lvth & Rodriguez! . l2005bl -1 Vernizzi & Wands, 



2007; 



200' 



2006bll; 

3. If the infla- 



tionary perturbation has a speed of sound differe n t from unity, th e n larg e non-Gaussianities 



may also appear (e.g. ISeery & Lidsevl ll2005bh : 



Lo Verde et al. 



Il2008ln . although in this 



case it is difficult to simultaneously achieve scale invariance. The current observational 
constraint, as we mentioned in Chapter [2l is of order |/ NL | < 100. In the absence of 
a detection, the forthcoming PLANCK mission may tighten this constraint to |/nl| ^ 3 



In Maldacena's normalisation, the numerical prefactor in Eq. J3 . 5 1 is not consistent with the square of the two-point 
function, Eq. l!3.4t. We choose A s o that the prefactor becomes 47r 2 (27r) 3 . This normalisation of Eq. J3 . 5 b was also 
employed bv lSeerv & Lidsevl ll2005bl lall. although the distinction from Maldacena's A was not pointed out explicitly. 
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Komatsu & Spergelll200ll : lLiguori et all 1200611 . 

Non-Gaussian PDFs have been studied previously by several authors. The closest analy- 
sis to the method developed in this chapter comes from lMatarrese et all B2000J], who worked 



with a path integral expression for the de nsity fluctuation smoothed on a s cale R (which 
they denoted by '<5r'). Also, the analysis of 



Bernardeau & Uzan 



1120021 1200311 has some fea- 



tures in common with our own, being based on the cumulant generating function. Moreover, 
the expression for the probability density in those papers is expressed as a Laplace trans- 
form. Our final expression, Eq. (13.701) . can be interpreted as a Fourier integral, viz (13.151) . 
which (loosely speaking) can be related to a Laplace integral via a Wick transformation. De- 
spite these similarities, the corres pondence between the two analyses is complicated because 



Bernardeau & Uzan 



II2002L 1200311 work in a multiple- field picture and calculate a probability 
density only for the isocurvature field '5s', which acquires its non-Gaussianity via a mixing 
of isocurvature and adiabatic modes long after horizon exit, of which particular cases were 
presented in the Chapter[2l This contrasts with the situation in the present chapter, where we 
restrict ourselves to a single-field scenario and compute the PDF for the adiabatic mode 1Z. 
This would be orthogonal to 5s in field space and its non-Gaussianity is generated exactly at 
horizon exit. 

In the older literature it is more common to deal with the density fluctuation 5 P measured 
on comoving slices, rather than the curvature perturbation 1Z. For slowly varying fields, on 
sca les larger th an the horizon, 1Z and 5 P can be related in the comoving gauge via Eq. (25) 



of IILvth . 



1985 



aH 



5 n = - 



w 



(3.8) 



to first order in cosmological perturbation theory for a barotropic fluid. (One may use the 
SN formalism to go beyond leading order as in Chapter [2l but to ob tain results valid on 
sub-hor izon scales one must use the full Einstein equations; see, e.g.. Langlois & Vernizzi 
Il2005b|]an .) For fluctuations on the Hubble scale (k ~ aH), this means \1Z\ ~ 8 p , so 1Z pro- 
vides a useful measure of the density fluctuation on such scales. By virtue of this relationship 
with the density fluctuation, the probability distribution F(g) is an important theoretical tool, 
especially in studies of stru cture formation. For exa mple, it is the principal object in the 
Press-Schechter formalism llPress & Schechten . 1197411 . As a result, there are important rea- 
sons why knowledge of the detailed form of the PDF of g, and not merely the approximate 
answer provided by the central limit theorem, is important. 

Firstly, large amplitude collapsed objects, s uch a s primordial black holes (PBHs) natu- 



1975 



Carr&HawkingLll974ll . Such 



rally form in the high-£ tail of the distribution HCarri . 
large fluctuations are extremely rare. This means that a small change in the probability den- 
sity for | £ | 3> ca n make a large difference in the mass f raction of the universe which 
collapses into PBHs [|Bullock & PrimackLll997l : llvanovLll998|1 . Thus one may hope to probe 
the form of the PDF for g using well-known and extremely stringent constraints on PBH 
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formation in the early universe JCan, 



Zaballa et all 



2007 



Josan et al. 



2003 



Carr & Lidsev , 



1993 



Green & Lid dle. 



1997 



2009]. The corrections calculated in this chapter are there- 



fore not merely of theoretical interest, but relate directly to observations, and have the poten- 
tial to sharply discriminate between models of inflation. 

Secondly, as described above, although the non-Gaussianities produced by single-field, 
slow-roll inflation are small, this is not mandatory. In models where non-Gaussianities are 
large, it will be very important to account for the effect of non-Gaussian fluctuations on struc- 
ture f ormation | Verde et al. . 2000 : Matarrese et alll200ol : Verde et al. . 2001 : Verde & Heavens 



200 111 . The formalism presented in this c hapter provides a system atic way to obtain such pre 



dictions, extending the analysis given by 



Matarrese et al. 



[2000]. 



The outline of this chapter is as follows. In Section l3~2l we obtain the probability measure 
on the ensemble of possible fluctuations. This step depends on the correlation functions of 
1Z. In Section 1331 we discuss the decomposition of 1Z into harmonics. This is a technical 
step, which is necessary in order to write down a path integral for F(g). First, we Fourier 
decompose 1Z. Then we write the path integral measure, and finally we give a precise speci- 
fication of g, which measures the size of fluctuations. We distinguish two interesting cases: 
a 'total fluctuation' g, which corresponds to 1Z (or approximately 5 P ) smoothed over regions 
the size of the Hubble volume; and the 'spectrum' V g (k), which describes the contributions 
to q from regions of the primordial power spectrum around the scale described by wavenum- 
ber k. In Section [3~4l we evaluate F(g). We give the calculation for the Gaussian case first, in 
order to clearly explain our method with a minimum of technical detail. This is followed by 
the same calculation but including non-Gaussian corrections which follow from a non-zero 
three-point function. In Section l331 we calculate F[V e (k)]. Finally, we summarise our results 
in Section [3761 



3.2 The probability measure on the ensemble of 1Z 



Our method is to compute the probability measure P*[R] on the ensemble of realisations of 
the curvature perturbation R(x), which we define to be the value of 1Z(t,x) at some fixed 
time t. This probability measure is a natural object in the Schrodinger approach to quantum 
field theory, where the elementary quantity is the wavefunctional \&t[R], which is related to 
Pt[R] by the usual rule of quantum mechanics, that Pt[R] oc |\I/ t [R]| 2 - Once the measure 
Pt[R] is known, we can directly calculate (for example) F t (g) by integrating over all R that 
produce fluctuations of amplitude g. Although the concept of a probability measure on R 
may seem rather formal, the Schrodinger representation of quantum field theory is entirely 
equivalent to the more familiar fo rmulation in term s of a Fock space. This representation is 
briefly discussed, for example, by lPolchinskil Iil998ll and lVisserl Iil996ll. A brief in troduction 
to infinite-dimensional probability measure s is given by Albeverio et aP 11997]. Indeed, 
a similar procedure has been discussed by llvanovl 11199811 . who calculated the probability 
measure on a stochastic metric variable ai s (x) which can be related to our 7Z(x). Although 
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the approaches are conceptually similar, our method is substantially different in detail. In 
particular, the present calculation is exact in the sense that we make no reference to the 
stochastic approach to inflation, and therefore are not obliged to introduce a coarse-graining 
approximation. Moreover, Ivanov's analysis appeared before the complete non-Gaussianity 
arising from 7^.-fie ld interactions around the time of horizon crossing had been calculated 
||Maldacenal . 1200311 • and therefore did not include this effect. 



3.2.1 The generating functional of correlation functions 

The expectation values (7£(xi) • ■ ■1Z(x.2)) m the vacuum |in) at some fixed time t can be 
expressed in terms of a Schwinger-Keldysh path integral^ 



[diz^diz + ] ^+ {t ' x)=7M< ' x) 7e(t, Xl ) 



H(t, x n ) exp (i J (n) [TZ+] - i J (n) [n. 



(3.9) 



Here [d1Z+] is the integrand of the path integral over 1Z. In cosmology we are generally 
interested in TZ evaluated at different spatial positions on the same t-slice, so we have set all 
the t equal in (13.91) . The path integral is taken over all fields TZ which begin in a configuration 
corresponding to the vacuum |in) at past infinity. The correlator is equal to the expectation 
value of three copies of 1Z at time t so we require two path integrals: the first integral [d7Z + ] 
evolves the vacuum state |in) from past infinity to the state 7Z + (t,x) at time t where the 
n copies of the field 1Z are averaged, and a second path integral [dlZJ\ which will project 
back the average to the vacuum state through a second functional integral . \TZ] is the 
action for the fluctuation 7Z, which is computed perturbatively to order n in 1Z when we 
want to com pute correlations o f the same order. For example, 1^[TZ] i s given to third 
order in 1Z by iMaldacenal (1200311 in the context of slow-roll inflation and by ISeery & Lidsey 
||2005allb|l in the inflationary models where the kinetic energy is not negligible. The action 
[R] is time ordered for the argument 1Z_ and anti-time o rdered for 1Z+. (For d e tails of 
the Sc h winger-Ke 



dysh o r 'closed t i me pa t h' forma l ism, s e e ICa lzetta & Hul 11987b ; 



ril986h : IWeinbergl l!2005h : lLe Bellad l!2000h : lHajicekl I11979ll : 



Rivers 



[1988].) 



Jordan 



An expression equivalent to Eq. (13.91) can be given in terms of the 'equal time' generating 
functional 



Z t [q] = I [dR] / [d7^ + £f x)=R(x) 



(3.10) 



exp J (n) [K+] - i X (n) [KJ\ + i J d 3 x R(x)g(x) J 



where q is some arbitrary source field, also known in the theory of special functions as 



6 Henceforth, we use the notation (• ■ ■ ) to mean expectation values in the in- vacuum, and no longer write in) explicitly 
where this is unambiguous. 
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the formal argument of the generating function. £ t is a spatial slice at coordinate time 
t. The equal-time correlation functions (TZ(t, xi) • • • lZ(t, x n )) are recovered from Z t [q] by 
functional differentiation, 



(^,xO--^(t,xJ)=I^...^ln^] 



(3.11) 

g=0 



Up to normalisation, this is merely the rule for functional Taylor coefficients, so it is straight- 
forward to invert Eq. (13.1 II) for Z t [q\. We obtain 

Zt[q] = exp{j2^ J •• y d 3 xi---d 3 x n , g(x 1 )---g(x n )(^(t,x 1 )---^(t,x n ))}. 



(3.12) 

Eq. (13.101) for the generating functional can be rewritten in a suggestive way. We define 
the wavefunctional at time t as 

%[R] = j [d7^ x)=R(x) exp (i J<»>[ft]) . (3.13) 

This definition is simply the functional generalisation of the familiar quantum-mechanical 
wavefunction. It expresses the amplitude for the field 1Z(t, x) to have the spatial configura- 
tion R(x) at time t, given the boundary condition that 1Z started in the vacuum state in the 
far past. In terms of \Pt[R], the generating functional can be rewritten as 

Z t [q] = J [dR] ^[R] f ^[R] exp M j d 3 x R(x) 9 (x)J = |t£[R]j 2 cc P[R], (3.14) 

where a tilde denotes a (functional) Fourier transform, and f denotes Hermitian conjugation. 
Eq. (13.141) implies tha t Z t \q\ is the complementary function for the probability distribution 
Pt[R] PAlbeverio et allll997|l . which can formally be obtained by inversion of Z t [q}. Hence, 



up to an overall normalisation, 



P t [R] oc J [dq] exp f-i J d 3 x R(x)g(x)J Z t [q\. (3.15) 

The normalisation is not determined by this procedure. We will fix the R-independent pref- 
actor, which correctly normalises the PDF, by requiring Jd^P(^) = 1 at the end of the 
calculation. For this reason, we systematically drop all field-independent prefactors in the 
calculation that follows. 



3.2.2 The probability density on the ensemble 

So far, all our considerations have been exact, and apply for any quantum field 1Z(t,x). 
For any such field, Eq. (13.151) gives the probability density for a spatial configuration R at 
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time t, and implies that to obtain F t [R] we need to know all such functions for all n-point 
correlations and at all spatial positions x. In practice, some simplifications occur when 1Z is 
identified as the inflationary curvature perturbation. 

The most important simplification is the possibility of a perturbative evaluation. The 
dominant mode of the CMB fluctuation is constrained to be Gaussian to high accuracy, so the 
non-Gaussian corrections to the leading order cannot be large. Moreover, the amplitude of 
its spectrum is constrained by CMB observations. Specifically, as mentioned in Section [L2l 
in the range of wavenumbers probed by the CMB, the spectrum has amplitude V^( 2 ~ 1CT 5 . 
Therefore, each higher-order correlation function is suppressed by an increasing number of 
copies of the spectrum, Vn(k), as we have shown in Chapter[2l 

Provided the amplitude of V-r, is small, it might seem reasonable to truncate the expo- 
nential in Eq. (13.151) for a given n and to work with a perturbation series in V. However, this 
simple approach is too nai've, because the integrals over q eventually make any given term in 
the series large, and this invalidates simple perturbative arguments based on power-counting 
in Vn- The perturbation series can only be justified a posteriori, a point to which we will 
return in Section l3.4.2l 

We work to first-order in the three-point correlation, that is, we consider non-vanishing 
two and three-point correlators in Eq. (13.121) . 

Z t [q]=ex V {-^ J y"d 3 xid 3 x 2 g(x 1 )g(x 2 )(^(t,x 1 )^(t,x 2 ) (3.16) 

~\J J y'd 3 x 1 d 3 a; 2 d 3 X3g(x 1 )g(x 2 )g(x3)(^(t,x 1 )^(t,x 2 )^(t,X3))}. 

This generating functional is introduced in the expression derived for the probability, Eq. (13.151) . 
for which we expand the exponential third order term in q(x) in power series to lowest order. 
We finally arrive to the product 



P t [R]oc J[dq]r t [q]LU t [q;R], (3.17) 
where T[q] and uj[q; R] are defined by 

T t [q] = (l - ~J d " kl ^f h qikMkM^intMWMWM)) , 



(3.18) 



and 



r ~! ( [ d 3 k 1 d 3 k 2 o(ki)g(k 2 ) , Sm/ , s . f d 3 k ., A 

u t [q;R] = expl-J 1 ^ l -f 2 ^ (n(t, kOftfo k 2 )) - i J — 9 (k)R(k)J 



(3.19) 



The expression for uo t gives rise to the Gaussian part of the PDF. T is of the form 1 plus a 
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correction which is small when the perturbative analysis is valid. Higher-order perturbative 
corrections in V-R can be accommodated if desired by retaining higher-order terms in the 
power series expansion of the exponential in (13.151) . Therefore our method is not restricted 
to corrections arising from non-Gaussianities described by three-point correlations, but can 
account for non-Gaussianities which enter at any order in the correlations of TZ, limited only 
by the computational complexity. However, in this chapter, we work only with the three- 
point non-Gaussianity. 

We now complete the square for uj t [q; R] in (13.191) and make the finite field redefinition 

g (k)» W = 9 (k) + (2^i ( ^ k ^_ k)> , , (3.20) 

where the prime in (lZ(t, \t)7Z(t, — k))' indicates that the momentum-conservation 5-function 
is omitted. The measure [dq] is formally invariant under this shift, giving J[dq] = f[dq], 
whereas u t [q; R] can be split into an R-dependent piece, which we call T t [R], and a piece that 
depends only on q but not R, 

u t [q; R] » r,[R]exp (~ J g(k 1 )g(k 2 )(^(t, k^f, k 2 )>) , (3.21) 

where T t [R] is a Gaussian in R, 

r,[R] = exp (-1 / d%d% mMmM)> m ^^_ K)) ) • (3.22) 
Eq. (13.171) for the probability density becomes 

P t [R] oc r 4 [R] J [dq] T t [g]exp J g(k 1 )g(k 2 )(^(t, ^(t, k 2 )}) , 

(3.23) 

One can easily verify that this is the correct expression, since if we ignore the three-point 
contribution (thus setting T t = 1), one recovers (after applying a correct normalisation) 



[dR] R(k 1 )R(k 2 )T t [R] = (^kO^ka)). (3.24) 



The remaining task is to carry out the q integrations in T t . The only terms which contribute 
are those containing an even power of q, since any odd function integrated against e~ q 
vanishes identically. In the expansion of f| i q(ki) in terms of q, there are two such terms: 
one which is quadratic in q, and one which is independent of q. These are accompanied by 
linear and cubic terms which do not contribute to P<[R]. For any symmetric kernel K and 
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vectors p, q e M m , one has the general results [IRiversl 1 1 9 8 811 

1 



[df] exp( -- / d m xd m y/(x)/(y)K(x, y ; 



= (det Kr 1/2 , 

[d/] /( P )/(q)exp f-1 1 d m xd m |//(x)/(y)K(x,y)j = K^faq) (det K) 



(3.25) 

1/2 _ 

(3.26) 



These rules allow us to evaluate the q integrals in Eq. (13.231) . giving 

p t [R]ocr\[R] (i + r! 0) [R] + rf[R]), 



(3.27) 



where 



T (o )[R] 



Tj 2) [R] 



d 3 ^ d 3 k 2 d 3 k 3 (TZ(t, kjKit, k 2 )TZ(t, k 3 )> 



R(k 1 )R(k 2 )R(k 3 ) 

(3.28) 



d 3 ^ d 3 £; 2 d 3 k 3 (R(t, ki)ft(t, k 2 )K(t, k 3 )) x (3.29) 
R(k 1 )^(k 2 + k 3 ) 

X n ¥3 (^,k i )^,-k i )>' +{perms} - 

In the last expression we include the possible permutations of the labels {1, 2, 3} since these 
give rise to distinct integrands. 

(2) 

In fact, T t is negligible. This happens because the three-point function contains a 
momentum-conservation 5-function, 5(k 1 + k 2 + k 3 ), which requires that the vectors kj sum 
to zero in momentum space. [For this reason, it is often known as the "triangle condition", 
and we will usually abbreviate it schematically as 5(A).] In combination with the 5-function, 
<5(k 2 + k 3 ), the effect is to constrain two of the momenta (in this example k 2 and k 3 ) to be 
equal and opposite, and the other momentum (in this example, ) to be zero. This corre 



2003 



Creminelli & Zaldarriaga . 



sponds to the extreme lo cal or 'squeezed' limit ||Maldacenal . 
20041 : lAUen et all 1200611 . in which the bispectrum reduces to the power spectrum evaluated 
on a perturbed background, which is sourced by the zero-momentum mode. Written explic- 
itly, Tj 2 ^ behaves like 



tJ 2) [R] 



6 



d 3 h d 3 k 2 
(2tt) 3 



aR(k 1 )5(k 1 ) + {perms}, 



(3.30) 



where we have written lim A = ak 3 ,, for some known constant a. In particular, Eq. (13.301) 

vanishes, provided R(k) approaches zero as k — > 0. This condition is typically satisfied, 
since by construction R(k) should not contain a zero mode. Indeed, any zero mode, if 
present, would constitute part of the zero-momentum background, and not a part of the 
perturbation R. Accordingly, Eqs. (I3.27I) - (I3.28I) with = give ¥ t [R] explicitly in terms 
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of the two- and three-point correlation functions. 



3.2.3 The smoothed curvature perturbation 

The probability density P t [R] oc (1 + Xj 0) [R])T t [R] expressed in Eq. (I3T23T) relates to the 
microphysical field 1Z(t, x) which appeared in the quantum field theory Lagrangian. A given 
k-mode of this field begins in the vacuum state at t — > — oo. At early times, the mode is far 
inside the horizon (k ^> aH). In this ('subhorizon') regime, the k-mode cannot explore the 
curvature of spacetime and is immune to the fact that it is living in a de Sitter universe. It 
behaves like a Minkowski space oscillator. At late times, the mode is far outside the horizon 
(k <C aH). In this ('superhorizon') regime, the k-mode asymptotes to a constant ampli 



tude, provided that only one field is dynamically relevant during inflation IILyth et all 12005 : 
Wands et all boodl H If we restrict attention to tree-level diagrams, then under reasonable 



conditions the integrals which define the expectation values of 1Z are typically dominated by 
the intermediate ('horiz on crossing') regime, where 7£(k) is exiting the horizon (k ~ aH) 



HWeinbergL I2005L 1200611 . As a result, the correlation functions generally depend only on the 
Hubble and slow-roll parameters around the time of horizon exit. 

The simple superhorizon behaviour of 1Z means that we can treat the power spectrum as 
constant outside the horizon. As has been described, its value depends only on the Hubble 
parameter and the slow-roll parameters around the time that the mode corresponding to k 
exited the horizon. For this reason, the time t at which we evaluate the wavefunctional 
\P*[R], the generating functional Z t [q] and the PDF P t [R] is irrelevant, provided it is taken to 
be late enough that the curvature perturbation on interesting cosmological scales has already 
been generated and settled down to its final value. Indeed, we have implicitly been assuming 
that t is the time evaluated in comoving slices, so that observers on slices of constant t see 
no net momentum flux. Because 1Z is gauge-invariant and constant outside the horizon, 
our formalism is independent of how we choose to label the spatial slices. The evolution 
of 1Z outside the horizon is the principal obstacle involved in extending our analysis to the 
multiple- field scenario. 

When calculating the statistics of density fluctuations on some given length scale 2-K/kn, 
one should smooth the perturbation field over wavenumbers larger than kjj. To take ac- 
count of this, we introduce a smoothed field R which is related to R via the rule R(k) = 
W(k, fc#)R(k), where W is some window function. The probabilities we wish to calculate 
and compare to the real universe relate to R rather than R. The exact choice of filter W is 
mostly arbitrary. For the purpose of analytical calculations, it is simplest to pick a sharp 
cutoff in k-space, which removes all modes with k < k H . Such window function is given by 

W(k,k H ) = e(k-k H ), (3.31) 



7 Where multiple fields are present, there will typically be an isocurvature perturbation between them: hype rsurfaces of 
constant pressure and density will not coincide. Under these circumstances 1Z will evolve IWands et all 2000]. We do not 
consider the evolving case in this chapter, but rather restrict our attention to the single-field case where the superhorizon 
behaviour of 1Z is simple. 
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where, Q(x) is the Heaviside step function defined in Eq. (12.1521) . This choice of window 
function has the disadvantage that it is non-local and oscillatory in real space, which makes 
physical interpretations difficult. The most common alternative choices, which do not suffer 
from such drawbacks, are (i) a Gaussian or (ii) the so-called 'top hat', which has a sharp 
cutoff in real space. We allow for a completely general choice of C° function W, subject 
to the restriction that W ^ except at k = oo and possibly at an isolated set of points 
elsewhere (we will work with more specific forms of the window function in the following 
Chapters |4] and [5]). This restriction is made so that there is a one-to-one relationship between 
R and R. If this were not the case, it would be necessary to coarse-grain over classes of 
microphysical fields R which would give rise to the same smoothed field R. 

In addition to this smoothing procedure, the path integral must be regulated before car- 
rying out the calculation in the next Section. This is achieved by artificially compactifying 
momentum space, so that the range of available wavenumbers is restricted to k < A, where 
A is an auxiliary hard cutoff or 'regulator'. At the end of the calculation we take A — > oo. 
Some care is necessary in carrying out this compactification. We set R = for k > A. In 
order to maintain continuity at k = A, we introduce a 1 -parameter family of functions Wa- 
These functions are supposed to satisfy the matching condition lim Wa(/c) = W(k), and 
are subject to the restriction Wa(A) = 0. (These conditions could perhaps be relaxed.) The 
relationship between R and R becomes 

R(k) = 0(A - k)W A (k; fcff)R(k) (3.32) 

To minimise unnecessary clutter in equations, we frequently suppress the A and kn depen- 
dences in W, writing only W(k) with the smoothing scale k H and hard cutoff A left implicit. 
Both the Gaussian and the 'top-hat' window functions approach zero as k — > oo, and are 
compatible with (13.321) in the A — > oo limit. In this limit, the final result is independent of 
the exact choice of W\(k, kn)- 

We are inter ested in the probability of observing a given fil tered field R. One can express 



this via the rule [|Matarrese et all 



20001 : lTaylor& Watts . 



200 



0] 



p t [R] = J [dR] P t [R]<5[R = 0(A - fc)WR]. (3.33) 
3.3 Harmonic decomposition of the curvature perturbation 

In the previous Section, we obtained the probability density for a given smoothed spatial 
configuration of the curvature perturbations. Given this probability density, the probability 
P that the configuration exhibits some characteristic of 1Z, such as fluctuations of amplitude 
q or a 'fluctuation spectrum' of the form V g (k), is formally obtained by integrating over all 



configurations of R which exhibit the criteria which define g IIMatarrese et all 1200011 . In this 



section, we give a precise specification of these criteria. Before doing so, however, we exploit 



3.3: Harmonic decomposition of the curvature perturbation 



77 



the compactification of momentum space introduced in (13.321) to define a complete set of 
partial waves. The smoothed field R can be written as a superposition of these partial waves 
with arbitrary coefficients. Moreover, the path integral measure can form ally be written as a 



product of conventional integrals over these coefficients [|HawkingLll977|1 . 

In the following we assemble the necessary formulae for the partial-wave decomposition. 
In particular, we will obtain expressions for the decomposition of R, for the characteristics g 
and V g (k), and a precise specification of the path integral measure. 

3.3.1 Harmonic expansion of R 

We expand R(k) in harmonics on the unit sphere and along the radial k — |k| coordinate: 

oo I oo 

R ( k ) = EEE R2|»*W*> <P)Mk)- (3-34) 

£=0 m=-l 71=1 

The Y im (9,(f)) are the standard spherical harmonics on the unit 2-sphere, while the *jj n (k) 
are any complete, orthogonal set of functions on the finite interval [0, A]. These harmonics 
should satisfy the following conditions§: 

1. ipn(k) — > smoothly as k — > 0, so that power is cut off on very large scales, and the 
universe remains asymptotically FRW with the zero-mode a(t), which was used when 
computing the expectation values (1Z • ■ ■ 1Z); 

2. ip n (k) — > smoothly as k — > A, so that the resulting R is compatible with Eq. (13.321) : 

3. ipn(k) should have dimension [M~ 3 ], in order that Eq. (13.341) is dimensionally correct; 

4. thetp n (k) should be orthogonal in the measure f Q A dk k 5 V^ 1 (k)W~ 2 (k). 

In addition, there is a constraint on the coefficients R^ n , because R(k) should be real in 
configuration space and therefore must obey the Fourier reality condition R(k)* = R(— k), 
where an asterisk denotes complex conjugation. The R™ n are generically complex, so it is 
useful to separate the real and imaginary parts by writing R™ n = a^ n + i6™ n . The condition 
that R is real in configuration space implies 

a~™ = (-lY +m a% n , (3.35) 
b~: = (-lY +m+ \ n . (3.36) 



8 When expanding functions on R 3 in terms of polar coordinates, a more familiar expansion involves spherical waves 
Zi m \k <x Jt(kr)Yi m (9, (f), where Je is a spherical Bessel function. These waves are eigenfunctions of the Laplacian in 
polar coordinates, viz, V 2 Z lm \ k = —k 2 Z lm \ k . An arbitrary function on R 3 can be written in terms of spherical waves, 
which is equivalent to a Fourier expansion. We do not choose spherical waves as an appropriate complete, orthogonal 
set of basis functions here because we do not wish to expand 'arbitrary' functions, but rather functions obeying particular 
boundary conditions at k — 0. The spherical waves for low £ behave improperly at small k for this purpose. Moreover, it 
is not possible to easily impose the boundary condition R(fc) — » as k — > A. 
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These conditions halve the number of independent coefficients, since the a and b coefficients 
with strictly negative m are related to those with strictly positive m, whereas for the m = 
modes, the b coefficients vanish if £ is even and the a coefficients vanish if i is odd. 

Condition Q] is made because, in the absence of this constraint, R could develop un- 
bounded fluctuations on extremely large scales, which would renormalise a(i). Therefore, 
condition Q] can be interpreted as a consistency requirement, since the inflationary two- and 
three-point functions are calculated using perturbation theory on a FRW background with 
some given a(t), which must be recovered asymptotically as |x| — > oo. It will later be 
necessary to sharpen this condition to include constraints on the behaviour of V-jz(k) near 
k = beyond the weak requirement that a 2 = J Vn{k) din k is finite. Condition |4] is 
a technical requirement made for future convenience. Any other choice of normalisation 
would work just as well, but this choice is natural, given the /^-dependence in the Gaussian 
kernel G[R]. Indeed, with this condition, the Gaussian prefactor in P(g) will reduce to the 
exponential of the sum of the squares of the a™ n and b^L. Condition [3] ensures that the inner 
product of two ip n (k) in the measure dk k 5 V^{k)W~ 2 {k) is dimensionless. Condition[2] 
has less fundamental significance. It follows from the condition Wa(A) = and the artifi- 
cial compactification of mo mentum space. However, as in the usual Sturm-Liouville theory 



flMorse & Feshbachi 1195311 . the precise choice of boundary condition is immaterial when 
A — »■ oo, so this does not affect the final answer. 

To demonstrate the existence of a suitable set of ip n (k), we can adopt the definition 

/, m V^2 Vn(k)W(k) T ( n k\ 

Mk) = - T TT-o J v \ a »T » ( 3 - 37 ) 



J v+1 {p%) Ak 2 \ u A 

where J v (z) is the Bessel function of the first kind and of order v, which is regular at the 
origin, and a" is its n-th zero. The order v is arbitrary, except that in order to obey conditionQ] 
above, we must have k u ~ 2 Vu(k) — > as k — > 0. This assumes that W(k) — >• 1 as k — > 0, 
as is usual for a volume-normalised window function. The i^ n {k) obey the orthonormality 
condition 

dk V n (k)W 2 {k) ^ n(yk ^ m(yk ^ = 6mn ' (3 ' 38) 
where 5 mn is the Kronecker delta. The completeness relation can be written 

k 5 

S(k - feo)|fce[ ,A] = p.g(fc )vv2(fc ) ^^( fc )^(fc ), (3.39) 

where the range of k is restricted to the compact interval [0, A]. 

Although we have given an explicit form for the ip n in order to demonstrate existence, 
the argument does not depend in detail on Eq. (13.371) . The only important properties are 
Eqs. (I338l) - (l339l) . which follow from condition H 
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3.3.2 The path integral measure 



Since any real C° function R obeying the boundary conditions R(k) — > as k — >• and 
R(k) — > as k — > A can be expanded in the form (13.341) . one can formally integrate over 
all such R by integrating over the coefficients R^ n . This prescription has been widely used 
for obtaining ex plicit results from path integral calculations. (For a textbook treatment, see 
Kleinertl [2004J].) In the present case, one should include in the integral only those R(x) 
which are real and so correspond to a physical curvature perturbation in the universe. Since 
the Yi m are complex, this means that instead of integrating unrestrictedly over the R^ ra , the 
reality conditions (13.351) must be respected. A simple way to achieve this is to integrate only 
over those a^ n or b^ n with m > 0. The m = modes must be treated separately since the a 
and b coefficients vanish for odd and even £, respectively. 
The integral over real R can now be written as 



|dR] 



nnn^ 

£=0 m=l ?i=l 



db% n 



oo oo 

nn^ 

r=0 s=l 
r even 



(3.40) 



where the subscript R on the integral indicates schematically that only real R(x) are included. 
The constants \i and ji account for the Jacobian determinant which arises in writing J[dR] in 
terms of the harmonic coefficients R™ n . Their precise form is of no importance in the present 
calculation as they will be absorbed by the final normalisation factor. 

As noted above, the detailed form of the measure (13.401) is not absolutely necessary for 
our argument. The important point is that each a or b integral can be carried out indepen- 
dently for m > 0. For this purpose, it is sufficient that the spectrum of partial waves be 
discrete, which follows from the (artificial) compactness of momentum space. However, al- 
though it is necessary to adopt some 'regulator' A in order to write the path integral measure 
in a concrete form such as (13.401) . we expect the answer to be independent of the specific reg- 
ulator which is chosen. In the present context, this means that our final expressions should 
not depend on A, so that the passage to the A — > oo limit becomes trivial. 



3.3.3 The total fluctuation g and the spectrum V e (k) 

There are at least two useful ways in which one might attempt to measure the amplitude of 
fluctuations in R. The first is the 'total smoothed fluctuation' at a given point x = x . By a 
suitable choice of coordinates, we can always arrange that x is the origin, so the parameter 
becomes g = R(0). When R is smoothed on scales of order the horizon size this gives a 



ume, since distan c es of l ess than a horizon size 
Shibata & Sasaki 1 1999h have proposed that g 



measure of the fluctuation in each Hubble vo 
no longer have any meaning. For example, 
defined in this way represents a useful criterion for the formation of PBHs, with form ation 
occurring whenever g exceeds a threshold value £> t h of order unity IIGreen et all 1200411 . This 
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measure of the fluctuation is non-local in momentum space. Making use of the relation 

/ dQ(9, cj)) Yt m (6, 4>) = \/47r^,o^m,o f° r me homogeneous mode of the spherical harmonics, 
one can characterise the amplitude as 

/d 3 k - \/Atx f 

^-pR(k)e ik - x | x=0 = J dk k 2 £ <\nMk). (3-41) 



On the other hand, one might be interested in contributions to the total smoothed fluctua- 
tion in each Hubble volume which arise from features in the spectrum near some characteris- 
tic scale of wavenumber k. For this reason, we consider a second measure of the fluctuation, 
which we call the 'fluctuation spectrum', defined by 

(Thus the total smoothed fluctuation can be obtained by integrating its spectrum according 
to the usual rule, viz, g — J V e (k) d In k.) This condition is local in k-space. Differentiating 
(13.411) . one can characterise V g (k) by the functional constraint 



V e {k) = ^-rj J2 a °o\nk 3 Mk). (3.43) 

^ ' n=l 

We will calculate the statistics of both the total fluctuation g and the spectrum V g (k). In each 
case, the calculation is easily adapted to other observables which are non-local or local in 
momentum space§. Indeed, both the non-local g and the local V g (k) are members of a large 
class of observables, which we can collectively denote by and which all share nearly- 
Gaussian statistics. Specifically, Eqs. (13.411) and (13.431) can be written in a unified manner in 
the form 

;>>° |n £ n (£;) = ^?(£0, (3.44) 



/47T 

n=l v 

where 

'/„ A d« 2 iW (* = <?); 345 

k 3 Mk) ($ = V e (k)). 

Note that, in the first case, the S n are independent of k. Any characteristic which can be put 
in this form, coupling only to the real zero-modes a^ n of R, will necessarily develop nearly- 
Gaussian (i.e., weakly non-Gaussian) statistics. More general choices of characteristic are 
possible, which cannot be cast in the form (13.441) . For example, one can consider charac- 
teristics which depend non-linearly on the a^ n . Such characteristics will generally lead to 



'The local and non-local variables defined here should not be confused with the local and equilateral triangulations 



which are specifically defined for the bispectrum in Chapter|2] 
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strongly non-Gaussian probabilities. The Gaussianity of the final PDF depends on the ge- 
ometry of the constraint surface in an analogous way to the decoupling of the Fadeev-Popov 
ghost fields in gauge field theory IIWeinbergL 1200511 . These non-Gaussian choices of char- 
acteristic can also be handled by generalising our technique, but we do not consider them 
here. 



3.4 The probability density function for g 



We first calculate the probability density for the non-local constraint g, given by Eq. (13.411) . 
The expression is 



oc / [dR] F[R]6 



(2tt) s 



n=l 



4ir 



(3.46) 



To obtain this density, one treats g as a collective coordinate parameterising part of R. The 
remaining degrees of freedom, which are orthogonal to g, are denoted by R ± . Therefore the 
functional measure [dR] can be broken into [dR- 1 ] and dg. After integrating the functional 
density P[R] [dR] over R- 1 , the quantity which is left is the probability density ¥(g) dg. In this 
case, the integration over the orthogonal degrees of freedom R- 1 is accomplished via the 5- 
function, which filters out only those members of the ensemble which satisfy Eq. (13.411) . We 
emphasise that this is a conventional 5-function, not a ^-functional. There is no need to take 
account of a Fadeev-Popov type factor because the Jacobian associated with the constraint 



(13.441) is field-independent, in virtue of the linearity of Eq. (13.411) in a 



On' 



3.4.1 The Gaussian case 

We first give the calculation in the approximation that only the two-point function is retained. 
In this approximation, the PDF of g will turn out to be purely Gaussian, which allows us to 
develop our method without the extra technical difficulties introduced by including non- 
Gaussian effects. 

If all correlation functions of order three and higher are set to zero, then we are in a 
Gaussian regime and hence P[R] oc G[R]. Using (13.41) . one can write 



G[R] = exp ( - i J dQ J k 2 dk 



1 



(27r) 3 27r 2 V n (k)W 2 (k) 



(3.47) 



£i,mi,m ^2, "12 ,ri2 



The harmonics Y em and ip n integrate out of this expression entirely, using the orthonormality 
relation (13.381) and the spherical harmonic completeness relation 



(3.48) 



3.4: The probability density function for g 



82 



Moreover, after rewriting the a and b coefficients with m < in terms of the m > coeffi- 
cients, we obtain 



oo I oo 



G[R] = exp ( - EEE 0«%| 2 + I^Tnl 2 ] " (3-49) 

v ' £=o m=l n=l 

. oo oo 

EED4J 2 +i&Wi 2 ])- 



47T 2 (27T) 3 

V > 1=0 71=1 

£ even 



The 5-function in (13.461) constrains one of the a[L (e.g. a^ ) in terms of g and the other 
coefficients. It would then be possible to evaluate f(g) by integrating out the 5-function 
immediately. However, this does not turn out to be a convenient procedure. Instead, we 
introduce the Fourier representation of the 5-function and rewrite (I3.46|) as 



¥(g) oc / [dR] ! dz G[R] exp ( iz 



\3 



(2tt 



v 47T 

n=l 



(3.50) 



where the functional measure is understood to be Eq. (13.401) . The final answer is obtained by 
integrating out z together with all of the a and b coefficients. In order to achieve this, it is 
necessary to separate a^ n , z and g from each other by successively completing the square in 
a°| and z. Working with a^ Q first, we find 



exp 



( i i 00 00 A 

"4^ (2^ E KJ 2 + a °o\nZ n (3.51) 

\ n=l n=l / 

f 1 1 fv 

^ 4vr 2 (2vr)3 ^ la °'" 



= 6XP I ~4^ (2% E(°o|n - 2^(2vr)^S„) 2 - ( 27 r) WE 2 

where we have introduced a function £ 2 = Y^n=i ^n- ^ n me fi na l PDF, £ 2 will turn out 
to be the variance of g. From Eq. (13.511) . it is clear that making the transformation a^ n i— > 
a^ n + 27r 2 (27r) 3 iz£ n suffices to separate a^ n from z. The measure, Eq. (13.401) . is formally 
invariant under this transformation. Exactly the same procedure can now be applied to z and 
g, giving 



exp ^-(27r) 3 7rV£ 2 - ^Li^ = exp -(27r) 3 7r 2 £ 2 (^z + 



10 Q 2 



2tt 2 v / 4^S 2 J 2S 2 
(3.52) 

As before, the finite shift z > 2 — i^/27r 2 v^47rS 2 leaves the measure intact and decouples z 
and g. The a, 6 and z integrals can be done independently, but since they do not involve g, 
they contribute only an irrelevant normalisation to P(f?). Thus, we obtain Gaussian statistics 
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for g: 



It remains to evaluate the variance £ 2 . In the present case, we have E n = J Q dk k 2 ifj n (k). 
From the completeness relation Eq. (13.391) . it follows that 

£ klUko)k 2 Mk) = ^f^^ - *b). (3-54) 

n 

£ 2 is now obtained by integrating term-by-term under the summation. The result coincides 
with the 'smoothed' conventional variance (cf. Eq. (13.31) ). 

Zl(k H )= / d\nkW 2 (k;k H )V n (k). (3.55) 
Jo 

Thus, as expected, Eq. (13.531) reproduces the Gaussian distribution (13.21) which was derived 
on the basis of the central limit theorem, with the proviso that parameters (such as S 2 ) 
describing the distribution of g are associated with the smoothed field R rather than the mi- 
crophysical field R. S 2 is therefore implicitly a function of scale, with the scale-dependence 
entering through the window function. Note that it was only necessary to use the complete- 
ness relation to obtain this result, which follows from condition |4] in Section l3-3.il 

3.4.2 The non-Gaussian case 

The non-Gaussian case is a reasonably straightforward extension of the calculation described 
in the preceding section, with the term in Eq. (13.271) now being included. However, some 
parts of the calculation become algebraically long, and there are subtleties connected with 
the appearance of the bi spectrum. 

The inclusion of T^ ^ corrects the pure Gaussian statistics by a quantity proportional to 
the three-point function, (IZIZIZ), which is given in Eq. (13.51) . This correction is written in 
terms of the representative spectrum V-r , which prescribes when the slow-ro ll prefactor, 



given by the amplitude of the spectrum, should be evaluated IIMaldacena . 



200311 . For modes 



which cross the horizon almost simultaneously, with size k\ ~ k 2 ~ k$, this prefactor should 



be V-r 2 = Vn(k) 2 , where k is the common magnitude of the k { . In the alternative case, where 



r-T 2 



one k-mode crosses appreciably before the other two, Vn should be roughly given by 

Vn = P^(max/c i )'P- R .(min/c i ). (3.56) 

Since the difference between this expression and the expression when all the k are of the same 
magnitude is very small, it is reasonable to adopt Eq. (13.561) as our definition of Vy 2 . We 



stress that this prescription relies on the conservation of 1Z outside the horizon jAllen et al. . 
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200611 . and it would therefore become more complicated if extended to a multiple- field sce- 
nario. 

With this parametrization, the probability measure on the ensemble is obtained by com- 
bining dH, (1527b . (1X281) and (1331) : 



1 f d^d^d 3 ^ gftU R(k 1 )R(k 2 )R(k 3j 

P[R] oc G[R] | 1 - - J (27r)627r2 8(A) vn[h) w{kimk2)w{h) 



(3.57) 



This expression should be integrated with the constraint (13.411) and measure (13.401) to obtain 
the probability F(g). At first this appears to lead to an undesirable consequence, since the 
integral of any odd function of R multiplied by G[R] must be zero. It may therefore seem that 
the non-Gaussian corrections we are trying to obtain will cancel out. This would certainly be 
correct if the integral were unconstrained. However, the presence of the 5-function constraint 
means that the shifts of a^ n and z which are necessary to decouple the integration variables 
give rise to a non- vanishing correction. 

The finite shift necessary to decouple a°, n and z is not changed by the presence of non- 
Gaussian corrections, since it only depends on the argument of the exponential term. This 
is the same in the Gaussian and non-Gaussian cases. After making this shift, which again 
leaves the measure invariant, the integration becomes 



¥(q) oc f [dR] f dz G[R]exp ( -(27r) 3 7r 2 £V 



(27r) 3 : 



47T 



izg) (I-J0-J2), (3.58) 



where 



Jo 



d 3 h d 3 k 2 d 3 k, 



27T 4 (27r)^ 



5(A) 



3(4vr) 3 / 2 



ni,n 2 ,n :j 



UiVn(h) 

^1(^1)^2(^2)^3(^3) 
W(k l )W(k 2 )W(h) ' 



(3.59) 



and 



J2 



d 3 fc 1 d 3 fc 2 d 3 fc 3j(A) VnA ^ ^2 



6(27r) 3 V47r 



(3.60) 



ni t2,rn,2,n2 I3,m,3,n2 



X izTi r 



^ni(k 



1 w(fci) R S™2 R Si«3 F ^ m 2^ 2 ' ^2)^3(^3, ^)^j: 



^2(^2)^3(^3) 



2)W(k 3 ) 

1^21 + 11^31. 



The symbol [1 ^ 2] represents the expression in square brackets with the labels 1 and 2 
exchanged, and similarly for [1^3]. The range of the m 2 and m 3 summations is from — £ 2 
to £ 2 an d —£3 to £ 3 , respectively. In addition, the shift of a^, n generates other terms linear 



and cubic in the R^ n , but these terms do not contribute to 



and we have omitted them 
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from (13381) . 

After shifting z to decouple z and g, the integrals Jo and J 2 develop terms proportional 
to z°, z, z 2 and z 3 . Of these, only the z° and z 2 survive the final z integration. Consequently, 
we suppress terms linear and cubic in z from the following expressions. The integral Jo 
becomes 



Jo 



d 3 hd 3 k 2 d 3 k 3 7r2(2?r) 



3(4vr) 



1 g 3 o^Q\ mm _ 



S ni E n2 E n3 W(A;i)n ; (A . 2)n ; (A . 3) 



ni,n2,n3 



while J72 simplifies to 

d 3 /ci d 3 /c 2 d 3 k 3 £> 



J2 



487r 3 (27r) 3 S 2 



E E E 

ni ^2,m2,n2 ^3, 1713,713 



x E 



"0m (A; 



■hi 



+R7? R?f K 



|n 2 V 3 |n 3 -'< ; 2m2> 



7713^3, <?3 



>n 2 (fc 2 )^n 3 (^3) 
W(k 2 )W(k 3 ) 

+ ri ^ 21 + ri - 



(3.61) 



(3.62) 



the m summations being over the same range as before. Thus J contains corrections pro- 
portional to q and g 3 , whereas J2 only contains corrections proportional to g. 

The a, b and z integrations can now be performed, with the integrand written entirely in 
terms of the a™ n and b^ n with m > 0. There are no a or b integrations in J . There are no z 
integrations in J 2 but the a and b integrations involved in the product R^ n2 R^f n3 A x ^2 = ^3, 
m 2 = m 3 and n 2 = n 3 . One then uses the spherical harmonic completeness relation, 



Yl Y tJPiM Y Ufi*<h) = S(tf>i - 2 )5(cos^ - cos£ 2 ) 

£=0 m=-l 

and the equivalent relationship for the ip -harmonics, Eq. (13.391) . to obtain 



(3.63) 



0J 2 



d 3 ki d 3 k 2 d 3 k 3 g 



VnA V n (h)W 2 (k 2 ) 



24vr S2 ^ 'UiVMWfc 



1c 3 



y^E re n (fei)^(k 2 + k 3 



[1 - 21 + ri ^ 31. 



(3.64) 



The terms with 1 exchanged with 2 and 3 generate the same integral as the first term and can 
be absorbed into an overall factor of 3. 
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J Q involves only z integrations. It can be written as 

f d% d% d% ( ff< , g\„ A , VkA 

x S ril E ri2 S ri 3?/; ni (/ci)?/; n2 (/c2)^n3(^3)- 

ni ,n 2 ,n 3 

To simplify these expressions further, it is necessary to obtain the value of the sum Yl^Li ^n^n 
Reasoning as before from the completeness relation Eq. (13.39k it follows that 

2_^ll n il) n {k) = — . (3.66) 

n=l 

From this, it is straightforward to show that 

f d 3 h d 3 k 2 d 3 k 3 - 2 / f? 3 e \ 



967r 3 ni*iW -1 (^) V s6 
where E 2 is the smoothed variance, Eq. (13.551) . On the other hand J 2 becomes 

j2 = j 2^/3 ^ 6{A) iMkV) wiki)A m ■ ( } 

After integrating out k 3 and the angular part of ki and k 2 , this gives 

J 2 = 2, J dfc 2 k\ / d*. S(k ,)±W(k l)l? ±_ lim 4^1, (3.69) 
where we have used the fact that k\ is constrained to zero by the 5-function to evaluate the 



bispectrum A in the 'squeezed' limit where one of the m omenta goes to zero IIMaldacenal . 



20031 : 1 Allen et all 120061 : ICreminelli & ZaldarriagaL l2004f| . In this limit, min(A; i ) = k\ and 



2 

max(/cj) = k 2 = k 3 , so it is possible to expand Vn unambiguously. Moreover, lim A = 

ak 3 ,, so J 2 = if Vfi{k)/k 3 — > as k — > 0. This more stringent condition on how strongly 
large-scale power is suppressed was anticipated in Section l3.3.1l It requires that Vn(k) falls 
at small k faster than k 3 . If this does not occur, then the integral diverges. (There is a 
marginal case when Vn(k)/k 3 tends to a finite limit as k approaches zero. We assume that 
this is not physically relevant.) 

The J2 integral contains a 5-function <5(k 2 + k 3 ). It can therefore be interpreted as 
counting contributions to the bispectrum which come from a correlation between the modes 
k 2 and k 3 , in a background created by ki, which exited the horizon in the asymptotic past. 
As we have already argued, modes of this sort are included in the FRW background around 
which we perturb to obtain the correlation functions of 1Z, so we can anticipate that its 
contribution should be zero, as the above analysis shows explicitly. In this interpretation, the 
condition Vn{k)/k 3 — > as k — > ensures that the perturbation does not destroy the FRW 
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background. Indeed, fluctuations on very large scales in effect describe transitions from one 
FRW world to another via a shift in the zero-momentum modes of the background metric. 
In this case, there is only one such mode, which is the scale factor a(t). These transitions are 
rather like changing the vacuum state in a quantum field theory. As a result, fluctuations of a 
large volume of the universe between one FRW state and another are strongly suppressed. 
For fluctuations on the Hubble scale, therefore, the PDF should be 



P(f?) 



exp I ) > (3-70) 



where we have used the fact that the corrections are odd in g and therefore do not contribute 
to the overall normalisation of F(g). The (dimensionless) coefficient J is 

f d'hd 3 k,<l>k 3 „. w -2, ..... 

3 = .1 9^n,ktw-'(k</ [A)Vk x <3 ' 71) 

This explicit expression is remarkably simple. Although it is preferable for calculation, it 
can be recast directly as the integrated bispectrum with respect to W: 



J = 48(27r) 3 (27r2) 3 J d 3 hd 3 k 2 d 3 k 3 (^(k 1 )7e(k 2 )^(k 3 ))>V(fc 1 )>V(fc 2 )>V(A ;3 ). 

(3.72) 

As a consistency check, we note that the expectation of g, defined as E(g) = J g¥(g) dg, 
is zero. This is certainly necessary, since the universe must contain as many underdense 
regions as overdense ones, but it is a non-trivial restriction, since both the g and g 3 corrections 
to ¥(q) do not separately average to zero. The particular combination of coefficients in (13.701) 
is the unique correction [up to 0(g 3 ), containing only odd powers of g] which maintains 
E(f?) = 0. 

Finally, we note that Eqs. (13.70l) - (13.71l) do not explicitly involve the cut-off A, except as 
a limit of integration in quantities such as £ A and Wa which possess a well-defined, finite 
limit at large A. As a result, there is no obstruction to taking the A — > oo limit to remove the 
regulator entirely. 



3.4.3 When is perturbation theory valid? 

It is known from explicit calculation that the bispectrum is of order V\ multiplied by the 
quantity, /nl, which is predicted to be small when slow-roll is valid. It is therefore rea- 
sonable to suppose that whenever the window functions W are peaked around some probe 
wavenumber k+, one has the order of magnitude relations S 2 ~ V± and J ~ V%, where 
represents the spectrum evaluated at k = Since the g 3 correction dominates for g > \/3E, 
this means that for g not too large, , the perturbative correction we have calculated 

will be small. As g increases, so that g ^> V+ 3 ^ 2 , perturbation theory breaks down and the 
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power series in g needs resummation. In any case, at such large values of g, the calculation 
described above ought to be supplemented by new physics which can be expected to become 
important at high energy densities. The details of these corrections presumably do not mat- 
ter too much, because at any finite order, the fast-decaying exponential piece suppresses any 
contributions from large values of g. 

At some value of g, corrections coming from the trispectrum can be expected to become 
comparable to those coming from the bispectrum that we have computed. We do not know 
preciselly which are the dominant contributions of the correction from the trispectrum. Such 
corrections to the non-Gaussian PDF are to be explored in the future. 



3.5 The probability density function for V Q (k) 

The probability density function for V e (k) can be obtained by a reasonably straightforward 
modification of the above argument, taking account of the fact that the constraint, Eq. (13.431) 
is now a functional constraint. This means that, when splitting the functional measure [dR] 
into a product of \dV s (k)] and the orthogonal degrees of freedom [dR- 1 ], the result after 
integrating out the R 1 - coordinates gives a functional probability density in \dV g (k)}. In 
particular, with the definition (13.451) . the 5-function in Eq. (13.501) is now represented as 



[dz] exp 



(2tt) ; 



i j( dk z(k) f J2<\ n k 3 Mk) - ^= V o(k) 



(3.73) 



In order to carry out this calculation, we write z{k) formally as 



*(*) = E 



z n ip n (k). 



(3.74) 



The integration measure J[dz] becomes Y[ n ft dz n , where, as before, ft is a field-independent 
Jacobian representing the change of variables from z(k) ^ z n . Its value is not relevant to 
the present calculation. In addition, we introduce a set of coefficients V e to describe V g (k), 



k 3 



oo 



n=l 



V en Mk). 



(3.75) 



The V Qn can be calculated using the rule V 6n = J Q A dk k 2 V^(k) W~ 2 (k)V e {k)ip n (k). Note 
that, in order to do so, we have made the implicit assumption that V e (k) /k 3 — > as k — ► 0, to 
ensure that (13.751) is compatible with the boundary conditions for the ip n (k). In other words, 
we make the ansatz of a suppression of power in modes with low k. 
With these choices, the 5-function constraint becomes 



n 



pL I dz n exp 



m=l 



(3.76) 
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In contrast to the nonlocal case of g, where a single extra integration over z coupled to g, 
we now have a situation where a countably infinite tower of integrations over z n couple to 
the the coefficients V e . In all other respects, however, the calculation is much the same as 
the nonlocal one, and can be carried out in the same way. The shift of variables necessary to 
decouple a^ n and z n is 

a° h^ag, +i27r 2 (27r) 3 z n ; (3.77) 



and the shift necessary to decouple the z n and V g is 

z n » z n - - en . (3.78) 
27r 2 v47r 



When only the two-point function is included, we obtain a Gaussian in the V Bn , 
PpV*)]«expf~^A (3-79) 

The sum over the V Qn can be carried out using Eq. (13.751) and the completeness and orthog- 
onality relations for the 4> n {k): 

Using this expression, and integrating over all V e (k) which give rise to a fluctuation of 
amplitude g, one recovers the Gaussian probability profile Eq. (13.531) with variance given by 
Eq. (13.551) . This serves as a consistency check for Eqs. (13.791) and (13.531) . 

When the non-Gaussian correction T^ ^ is included, one again generates a probability 
density of the form 



n / 



nV e (k)} <x{l-K - K 2 ) exp -- , (3.81) 



where K 2 has the same form as J%, and therefore vanishes for the same reasons, and 

K ° J iiwPi ((A)Ps <3 ' 82) 



X 



P P ( 2 + 3) " 11 



The first term contains a 5-function which squeezes k\ into the asymptotic past. It formally 
vanishes by virtue of our assumption about the behaviour of V g (k) near k = 0, which is 
implicit in Eq. (13.751) . As a result, the total probability density for the fluctuation spectrum 
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can be written as 



F[V e (k)] oc (l-K)exp 



( 



1 



/ 



din A; 



) 



(3.83) 



2 



Vn(k)W 2 (k) 



where 



K 



1 



d 3 h d s k 2 d 3 k 3 



(3.84) 



96?r 2 YliVnihWiki) 



As before, one can show that this expression is consistent with Eqs. (13 .701) — (13 .711) by integrat- 
ing over all V g (k) which reproduce a total fluctuation of amplitude g, after dropping another 
term which is squeezed into the asymptotic past owing to the presence of a 5-function. This 
is a non-trivial consistency check of Eqs. (13 .831) — (13 .84-1) . 

As in the local case, Eqs. (I3.83I) - (I3.84I) are entirely independent of A (except as a limit 
of integration), so the regulator can be freely removed by setting A = oo. 

3.6 Summary of results 

In this chapter we have obtained the connection between the n-point correlation functions of 
the primordial curvature perturbation, evaluated at some time t, (7Z(ki) • • • 7Z(k n )), and the 
PDF of fluctuations in the spatial configuration of 1Z. We have obtained an explicit expres- 
sion for the PDF of a fluctuation of amplitude g when 1Z is smoothed over regions of order 
the horizon size. This is a probability density in the conventional sense. In addition, we 
have obtained an expression for the probability that g has a spectrum V e (k). This is given by 
/ d In k V e {k) = q, although mapping V g (k) i— > g is many-to-one. This is a functional prob- 
ability density, and can potentially be used to identify features in the fluctuation spectrum 
near some specific wavenumber k ~ k*. Our result is independent of statistical reasoning 
based on the central limit theorem and provides a direct route to incorporate non-Gaussian 
information from the correlators of the effective quantum field theory of the inflaton into 
theories of structure formation. 

Both these probabilities are Gaussian in the limit where 1Z only possesses a two-point 
connected correlation function. If there are higher-order connected correlation functions, 
then 1Z exhibits deviations from Gaussian statistics, which we have explicitly calculated 
using determinations of the inflationary three-point function during an epoch of slow-roll 
inflation. Our method can be extended to incorporate corrections from higher connected n- 
point functions to any finite order in n. We have not computed these higher corrections, since 
we anticipate that their contribution is subdominant to the three-point correction (which is 
already small). 

Our argument is based on a formal decomposition of the spatial configuration of the 
curvature perturbation in k-space into spherical harmonics, together with harmonics along 
the radial k direction. However, we have emphasised that our results do not depend on the 
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details of this construction, but require only a minimal set of assumptions or conditions. 
The first assumption is that the power spectrum Vn(k) goes to zero sufficiently fast on large 
scales, specifically Vn(k)/k 3 — > as k — > 0. (In addition, in the case of the fluctuation 
spectrum, we require V g (k)/k 3 — ► as k — > 0.) Such a condition is certainly consistent 
with our understanding of large-scale structure in the universe and, within the perturbative 
approach we are using, we have argued that in fact it describes a self-consistency condition 
which prevents perturbative fluctuations from destroying the background FRW spacetime. 
Our second assumption is that the spatial configuration R can be smoothed to R via a window 
function W to obtain a configuration for which R — > as k — »■ oo. In this case, it is fair to 
compare R to the primordial power spectrum. 

In addition to these fundamental assumptions, which relate to the behaviour of real phys- 
ical quantities, a large part of the calculation has relied on an auxiliary technical construc- 
tion. This construction is based on an artificial compactification of momentum space, im- 
plemented by a hard cutoff A. There is an associated boundary condition on R at k = A 
which discretises the harmonics (partial waves) in k. However, in both the non-local (to- 
tal fluctuation g) and local (fluctuation spectrum V e (k)) cases, the final probability density is 
independent of both the details of the partial wave construction and A (except as a limit of in- 
tegration). It is also independent of the choice of the family of window functions W\(k; ku), 
and depends only on the limit lim W^{k] k H ) = W(k; k H ). Therefore the regulator can be 

A^oo 

removed by taking the limit A — > oo. Moreover, the boundary condition at k = A becomes ir- 
relevant in this limit, which is a familiar result from the theory of Sturm-Liouville operators. 
As a consistency check, one can integrate P[P £ ,(fc)] with the condition J d In k V e {k) = g in 
order to obtain F(g). 

In Chapters ID and [5] we present applications of this method, and the PDF obtained, to 
improve the estimation of the probability of PBH formation. 



Chapter 4 

Probability of primordial black hole 
formation 



4.1 Introduction 

Primordial Black Holes (PBHs) are a unique tool to probe inhomogeneities in the early uni 
verse. The probability of PBH formation is extensively studied because it is use ful in con 



straini n g the amplitude of primordial inhomogeneities generated by inflation (e.gJCarr et al 



[ 1 994] ; L iddle & Greenl 11998] ; 



Sendouda et al. 



[2006]; 



Zaballa et al. 



D2007D; 



Bugaev & Klimai 



(20Q6D). What makes PBHs a unique tool in cosmology is the range of scales that can be 
probed by their formation. The anisotropics probed by the CMB data cover the range of 
wavenumbers 7 x 10~ 4 < fc/Mpc -1 < 0.021. Equivalently these modes enter the horizon 
when the cosmological horizon or Hubble mass is between 10 19 < M/M Q < 10 23 while the 
overdensities forming galactic haloes have associated masses 10 8 < M/M Q < 10 12 . The 
inhomogeneities forming PBHs are much smaller and they can span the range of wavenum- 
bers 10 3 < fc/Mpc" 1 < 10 16 which correspond to 10" 24 < M/M Q < 10 6 , a set of values 
that can change with the model of inflation and it's reheating scale. In any case, this is the 
largest range of scales probed by any single observable in the universe. 

Another advantage of studying PBH statistics is that, in a radiation background, the grav- 
itational collapse of fluctuations takes place shortly after horizon crossing. Consequently, 
PBH statistics do not suffer the bias problem or the late-time nonlinear evolution that signif- 
icantly modifies the mass and statistics of other bound objects. 

The absence of direct detections of PBHs has prompted studies of processes that could be 
influenced by the gravitational effects of PBHs or their evaporations. Some of the processes 
and observations that limit the abundance of PBHs are the following: 

1 . If the number of PBHs is large enough, they could constitute a significant fraction of 
the dark matter. The current density of PBHs therefore cann ot exceed the observe d 
density of dark matter, i.e., ^pbh(M > 10 15 g) < tl DM = 0.28 IIKomatsu et allboOSn . 



2. The Hawking radiation from PBHs liHawkingL 1197411 can be the source of the back 
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ground radiation at various waveleng ths in our universe JCanl. 1 19761: IPage & Hawkins, 



1976 



Bugaev & KonishchevL 1200 1Q and cosmic rays flBugaev & KonishchevL 1200 1Q . 
As mentioned in Chapter [H PBHs of mass M cvap = 5 x 10 14 g should be evapo- 



rating tod ay and observations of the gamma-ray background imply f2 PR 



5 x 



0" 8 llPage & HawkingL 



1976] 



Carri . 



1976 



MacGibbon & Can 



1991 



Kimetal 



1 999|1 . This is the tightest constraint on the density of P BHs although future ob serva 



tions of the 21 cm radiation might impose a tighter limit IIMack & WesleyL 12008]. 



3. Black holes with mass M < M cvap have already evaporated and the decay products 
should not spoil the well understood chemical history of the universe. Indeed, limits 
on /3p BH (M) can be obtained in the mass range 10 9 < M/g < 10 12 by looking at the 
effects of hadrons and neutrinos emitted by PBHs on the Big Bang nucleosynthesis of 
helium and deuterium I Miyama & Sato . 19781 : Novikov et al. . 



19791. 



A com plete list of numerical bounds can be found in Table I, as compiled by 



Green & Liddle 



111 99711 . All these bounds have been used to probe early universe fluctuations llCarr et al 



1994 : lLiddle & GreenLI 19981 : ISendouda et all l200ft IZaballa et all 120071 : Bugaev & Klimai . 



2006]. They can be translated into limits on the root-mean-square amplitude of density or 
curvature perturbations 7£rms on scales inaccessible to the CMB. 

Here we explore how the bounds to 7£rms can be modified in view of the considera- 
tion of a non-Gaussian probability distribution. We use the PDF derived in Chapter [3] and 
calculate the mass fraction of PBHs with the aid of the Press-Schechter formalism. The ef- 
fects of non-Gaussian perturbations on PBHs have already bee n studied for specific models 



Bullock & PrimackL 1 19971: llvanov . 



1998; 



Pina Avelino 



2005] but a precise quantification 



of the non-Gaussian effects is still required. Indeed, it is only now, with a much better 
understanding of the effects of higher order perturbations, that we are able to describe the 
general effects on PBHs. This discussion is crucial in the light of recent claims that only ex- 
otic extensions of the canonical slow-roll inflationary potentials can produce a n appreciable 
number of PBHs llChongchitnan & Efstathioul 120071 : Bugaev & KlimaiL 12008(1 (see however 
[IPeiris & Eastherl 1200811 where it's argued that a large number of PBHs can be formed even 
within the slow-roll regime). Here we explore whether the consideration of non-Gaussian 
perturbations in inflationary models could increase the mass fraction of PBHs significantly. 
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Table I Constraints on the mass fraction /3 PBH (M) of the universe going into PBHs 



CONSTRAINT 


MASS RANGE (g) 


NATURE 


1.25 x 1U y 10 n g ) 


^ in 11 


entropy of the universe 


/ \ 1/2 
4 1 X 10" 3 1 -Mr- 1 


10 9 - 10 11 


pair-production at nucleosynthesis 


4.9xlO- 7 ( 10 A 4j 


10 10 - 10 11 


Deuterium destruction 


/ \ 7 / 2 

6.5x10 ^ 1011g J 


1 nil t nl3 
1U — 111 


Helium-4 spallation 


io- 18 (^) _1 


IQll - 10 13 


CMB distortion 


3.1 x 1CT 27 


3.6 x 10 14 - 10 15 


7-rays from evaporating PBHs 




> 10 15 


fipBH(to) < 0.47 



4.2 The non-Gaussian PDF 

Let us introduce the elements of the non-Gaussian distribution of probabilities for the cur- 
vature perturbation field 1Z. We first describe how the Gaussian PDF is constructed in the 
context of the linear theory. The amplitude of the curvature perturbations 1Z is derived by 
solving the perturbed Einstein equations to linear order. Statistically, the mean amplitude is 
written in terms of the two-point correlation function as 

(7e G (k 1 )7e G (k 2 )) = (2 7 r) 3 5(k 1 + k 2 )\n RMS (k)\ 2 , (4.1) 

where, as before, 7^c(k) are Gaussian perturbations in Fourier space. 

The two-point correlator defines the dimensionless power spectrum V(k) through the 
relation 

2tt 2 

(^ G (k 1 )^ G (k 2 )) = (2 7 r) 3 5(k 1 + k^—Vih). (4.2) 

As discussed in Chapter[3l the perturbations are smoothed over a given mass scale ku- Here 
we choose a truncated Gaussian window function 

W M (k) = 6(A; max - k) exp f-^rj , (4-3) 

were 6 represents the Heaviside function and the fiducial scale k mSLX is introduced to avoid 
ultraviolet divergences. The smoothing scale k M is defined by 



k M = 2nH M = M/2, 



(4.4) 
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where M is the Hubble mass at the time the scale k M enters the horizon. 
The variance of the smoothed field is related to the power spectrum by 

^WUk)V(k). (4.5) 

The power spectrum encodes important information about the underlying cosmological model. 
For example, in the case of perturbations deriving from the quantum fluctuations of a single 
inflationary fi eld <j> with a potential V dominating the cosmological dynamics, the explicit 



expression is llStewart & LyfhL 1 199311 



Hi V[ 

(2vr) 2 02 m 2 ~ [dV/d^lm] 



Wk) = ln „l ; « ,*„A^2 > (4-6) 



Here an asterisk denotes values at the time when the relevant perturbation mode exits the 
cosmological horizon, k = a*if* = a(t*)H(t*). 

The tilt of the power spectrum is parametrised with a second observable, the spectral 
index, which is defined as 

"■- 1 = ^i lnp w- (4 - 7) 

If n s < 1, the root-mean- square amplitude TZ^ms increases on larger scales, corresponding to 
a red spectrum. Conversely, n s > 1 indicates larger power on smaller scales and corresponds 
to a blue spectrum. 

The power spectrum and the tilt are derived directly from linear perturbations as reviewed 
in Section l2~4l of Chapter[2l In observations of the CMB, it is possible to determine with great 
accuracy the numerical values of the power spectrum and its tilt on scales larger than the 
horizon at the time of last scattering, that is (k < k\ s = 1.7 x 10~ 3 Mpc~ 1 ). On such scales, 
the five-year results o f WMAP, combined w ith the galaxy counts, give T(k\ s ) = 2.4 x 1CT 9 



and n s = 0.95 ±0.1 IlKomatsu et al. 



1200811 . 

In linear perturbation theory one makes use of the central limit theorem to construct the 
PDF. To first order, the perturbation modes are independent of each other. If we assume 
the field of linear perturbations TZ has zero spatial average, then the central limit theorem 
indicates that the PDF of TZ is a normal distribution which depends only on the variance E^, 

i / n 2 \ 

P G (ft) = -=— exp -— TJ — . (4.8) 



2^ n ^\ 2Y? n (M) 

A successful linear theory of structure formation will predict this probability distribution and 
match the numerical values at the relevant observational scales. Higher order correlations of 
the perturbation field 1Z offer an exciting way to distinguish between cosmological models 
with common properties at linear order. As discussed in Chapter [Hand Chapter [2l the devi- 
ations from Gaussianity are described to lowest order by the nonlinear parameter / NL . This 
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parameter appears in the expansion (e.g. iLyth & Rodriguez! I12005ah ) 



n(k) = n G (k) - -hL(K G *n G (k) - (n 2 G )), 



(4.9) 



where a star denotes the convolution of two copies of the field. The interaction of Fourier 
modes does not admit the use of the central limit theorem and the non-Gaussian probability 
distribution must be constructed by other means. 

In Chapter [3] we have provided a method to calculate the correction to the Gaussian PDF, 
and to derive a new PDF which includes the linear order contribution from the 3 -point func- 
tion. Such a correlator can be derived through a second order expansion of the perturbations 



2005 



Seery et al. . 



in the Einstein equations HBartolo et all 120041 ; iRigopoulos & ShellardL 
2008]. Alternatively, an explicit expression for the three-point correlator can be obtained 
from the third-order quantu m perturbations to t he E instein-Hilbert action. Pi oneering works 



using this method come fr omlMa ldacen al 000311 and 



Seery & Lidseyl ||2005bl]a|l . Here we use 



the expression derived by ILyth & Rodriguez! ll2005ah for the correlator in Fourier space. At 
tree-level this reduces to 



(K{h)n{k 2 )K(k 3 )) 



4vr 4 ^ NL 
5 



V(h)V(k 2 ) 

1,31,3 



{perms} 



(4.10) 



Current observations provide numerical bounds for /nl through the three-point correlation 
of the temper ature fluctuation mode s. The WMAP satellite gives the constraints —151 < 
/nl" 11 < 253 BKomatsu et all 1200811 fo r an equilateral triangulation of the bispectrum and 



_4 < /i°£ aI < 80 llSmith et alTEoOoll for a local triangulation (the local and equilateral 



triangulations have been defined in Sec. ll.3l) . Both of these values are determined at the 95% 
confidence level and consider an invariant value at all scales probed by the CMB. 

In the following, the basic components of the non-Gaussian PDF derived in Chapter [3] 
are presented. The amplitude of the perturbation is characterised by its value at the centre of 
the configuration 



0r 



ft(x = 0). 



(4.11) 



This specification is necessary to construct an explicit expression of the PDF. The param- 



eter # n is particu 



llShibata & Sasaki! . 



a rly useful to discriminate the relevant inhomogeneities forming PBHs 



1999 



Green et all I2004TI . The non-Gaussian probability distribution 



function for a perturbation with central amplitude i3 , derived in Eq. (13.701) . is 



Png(^c 



2ttS 



n L 



?9 3 
_o_ 

Y 3 



3tf \ J 



Y 3 



cxp 



l) 2 



?Y 2 



(4.12) 
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where the factor J encodes the non-Gaussian contribution to the PDF: 

J= l_j dkl ^ ka W M (A:i)>VM(A:2)WM(A: 3 )(^(k 1 )7e(k 2 )7e(k3)), (4.13) 
1 I" dki dk 2 dk 3 \ 

5 J n.mwj'it) v ; ' 



MM + {perms} 
kik 2 



(4.14) 



This last equation is valid at tree level in the expansion of (IZIZIZ). It is justified as long 
as the loop contributions to the three-point function, generated by the convolution of 1Z- 
modes, are sub-dominant. This requirement is met when the second order contribution to 1Z 
in Eq. (14.91 ) does not exceed the linear contribution. This is requirement is met if we demand 
that 



/nl ^ l/VW)- (4-15) 

The complete derivation of the PDF in Eq. (14.121) was already provided in Chapter [3] 
Here it is sufficient to say that the time-dependence of this probability is eliminated when 
the averaging scale is ku < a(t)H(t) providing the growing mode of the perturbation 1Z is 
constant on superhorizon scales. This is true in particular for perturbations 1Z considered in 
the radiation era, when the PBHs considered here are formed (see Chapter [I]). 

In order to adapt the PDF in Eq. (14.141) to the computation of PBH formation proba- 
bilities, this expression is integrated between the limits fc m i n and fc max defined to cover the 
relevant perturbation modes for PBH formation. PBHs are formed long before today, so 
in the large-box (small wavenumber) limit of integral (14.141) , t he pre sent Hubble horizon 
k min = H is a reasonable lower limit for PBH formation llLvth , [l 992 1 . At the other end of 



the spectrum, the smallest PBHs have the size of the Hubble horizon at the end of inflation. 
A suitable upper limit in this case is the wavenumber associated with the comoving horizon 
at the end of inflation, k max = a(t cn d)H en d. It is important to mention that, even though the 
integral in Eq. (14.141) should include all fc-modes, finite limits are imposed to avoid loga- 
rithmic divergences. Due to the window function factors Wm(A;), the dominant part of the 
integral is independent of the choice of integration limits as long as they remain finite. 

The integral (14.141) is considered only at the limit of equilateral configurations of the 
three-point correlator, that is, considering correlations for which ki = k 2 = k 3 . This is not 
merely a computational simplification. In the integral, each perturbation mode has a filter 
factor Wm(^) which, upon integration, picks dominant contributions from the smoothing 
scale k M common to all perturbation modes. In this case J can be written in the suggestive 
way: 

J=-o ~T~ lWM(k)V(k)} 2 — /nl • (4.16) 



ky. 



k 1 y ' v /J \5 



With the complete non-Gaussian PDF at hand, it is possible to characterise its effects 
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on the probability of PBH formation. In the next section, J is computed numerically for 
inflationary perturbations generated in a single-field slow-roll inflationary epoch. The results 
in this case are shown to be consistent with previous works on non-Gaussian computations of 
the probability of PBH formation. In Section l4~4l the non-Gaussian PDF is generated for the 
case of constant / NL . This will be used to test the magnitude of the effects of non-Gaussianity 
on the probability of PBH formation. 



4.3 Non-Gaussian modifications to the probability of PBH formation 



The simplest models of structure formation within the inflationary paradigm are those where 
a single scalar field drives the accelerated expansion of the spacetime and its quantum fluc- 
tuations evolve into the observed structure in subsequent stages of the universe. Although 
small in magnitude, the non-Gaussianity of the fluctuations generated in this simple model 
provide a qualitative hint to the consequences that non-Gaussianity has for the probability of 
PBH formation. 

In fact, for single-field inflationary models, the effect s of non-Gaussianity onPBHs have 
been explored in the past but with inconclusive results. iBullock & Primackl Ill997ll studied 
the probability of formation of PBHs numerically for non-Gaussian perturbations with a blue 
spectrum (n s > 0). The motivation for this was that any inflationary model with a constant 
tilt and a normalisation consistent with the perturbatio ns at the CMB scale must have a 
blue s pectrum to produce a significant number of PBHs llCarr et all Il994t iGreen & Liddlei 
1997(1 . Their analysis is based on the stochastic generation of perturbations on superhorizon 
scales, together with a Langevin equation for computing the PDF. For all the cases tested, 
the non-Gaussian PDF is skewed towards small fluctuations. In consequence, the probability 
of PBH formation, which integrates the high amplitude tail, is suppressed with respect to the 



Gaussian case. An example of the kind of potential studied by lBullock & Primackl [|1997|l is 



1 + arctan ( — 
1 + (4xl(F) U 



21 



>0, 
< 0. 



(4.17) 



where Vq is the amplitude of the potential at <\> = 0. This potential features a plateau for 
4> < 0. This produces an increase in the power of matter fluctuations corresponding to the 
production of PBHs of mass 10 32 gr. 

Another way of generating large perturbations in the inflationary scenario is to consider 
localised features in the potential dominating the dynamics regardless of the tilt of the spec- 
trum. As one can see from Eq. (14.61) . an abrupt change in the potential would generate a 
spike in the spectrum of perturbations. This is valid as long as we avoid a 'fiat' or 'static' 
potential in which dV/dtp = 0. In such case = and Eq. (14.61) i s invalid (for a treatment 
of this particular case, also known as 'ultra-slow roll' inflation, see lKinneyl Q2005D). 

The description of a model of inflation with large amplitude in the power spectrum is 
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incomplete if we do not take on account the effects of nonlinear fluctuations. The effects of 
non-Gaussianity for an inf lationary mod el producing features in an otherwise red spectrum 



(n s < 0) were explored by 



Ivanov 



, using the toy model 



Vo 



A(<j) 2 - 



for < <f>i, 
for 2 > (j) > 
for <fi > 02- 



(4.18) 



where A and A are coupling constants. Through a stochastic computation of the PDF, Ivanov 
found that the non-Gaussian P DF is skewed toward s large perturbations. This result goes in 
the opposite direction to that of Bullock & Primackl [1997]. 

To understand this difference and generalise the effects of non-Gaussianity, it is conve- 
nient to look at the fractional difference of the Gaussian and non-Gaussian PDFs: 



P 



NG 



P 



G 



Pf 



4- -3—1 4' 



(4.19) 



Both 



Bullock & Primackl 1199711 and 



Ivanov 



ll 199811 use perturbations generated in a piece- 
wise slow-roll inflationary potential for which inflation is controlled by keeping the slow-roll 
parameters, defined in Eq. (12.931 ), smaller than one. Here the slow-roll approximation is used 
to explore the qualitative effects of Eq. (14.191) . 

To linear orde r, there is a straightfo rward expression for the spectral index in terms of 



these parameters jStewart & Lyfhlll993|l . 



2fa 



SR 



3e SR ). (4.20) 

On the other hand, by using a first order expansion in slow-roll parameters, Maldacenal 1 2003 ] 
provi des an expression for the nonlinear factor /nl in terms of these parameters jMaldacenal . 
2003ll : 



5 5 
/nl =— (n s + F{k)n t ) = - (r) SK - 3e SR + 2J r (k)e S R) 
LZ 



(4.21) 



where n t = 2esR is the scalar-tensor perturbation tilt and F(k) is a number depending on the 
triangulation used. For the case of equilateral configurations, when JF = 5/6, 



NL 



^SR — g^SR 



(4.22) 



cq 



This last expression is used to evaluate the integral (14.161) for J. The non-Gaussian effect on 
the PDF is illustrated in Fig. 14.1 I for the potentials given by Eqs. (14.171) and (14.181) in terms of 
the fractional difference (14.191) . This difference represents the skewness of the non-Gaussian 
PDF. The non-Gaussian contribution encoded in the factor J is the integral of /nl over all 
scales relevant for PBH formation. Consequently the sign of /nl is what determines the 
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enhancement or suppression of the probability for large amplitudes $0 in the non-Gaussian 
PDF. For the two cases illustrated, the scalar tilt n s dominates over the tensor tilt n t , so that 
the sign of / NL coincides with that of n s . This result is illustrated in Fig. 14.11 
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Figure 4. 1 : The fractional departure from the Gaussian PDF is plotted for two types of non-Gaussian distri- 
butions Png> as defined in Eq. ( 14.12b . For the potential in Eq. (14.171 i, /nl > and the departure is plotted 
with a solid line. For the potential in Eq. (14.18b . /nl < and the departure is shown by a dashed line. 



4.4 Constraints on non-Gaussian perturbations of PBH range 



A stand ard practise in calculating the PBH mass fraction is to use the Press-Schechter for- 
malism [ Press & Schechter . il 974 1. As described in Chapter [T] this involves integrating the 
probability of PB H formation over the relevant matter perturbation amplitudes, 5, measured 
at horizon epoch llcarr . 1975 1 and gives 



/W > M) = 2 / F(6 P (M)) d6 p (M). 
Js th 

For the large ratio S^fEp this can be approximated as 



(4.23) 



/W > M) 



E p (M) 



exp 



S 2 



2£2(M)J 



(4.24) 



where Y? S (M) is the variance corresponding to the mass scale M and <5 t h is the threshold 
amplitude of the perturbation necessary to form a PBH. When the relevant amplitudes of 
a smoothed perturbation are integrated, /?pbh represents the mass fraction of PBHs with 
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M > w 3 ' 2 M H « w 3 / 2 k M /(2n) llCarrl . 1 197511 . where w is the equation of state at the time of 
formation. Note that the approximation (14.241) is valid only for a Gaussian PDF. 

The integral /3pbh establishes a direct relation between the mass fraction of PBHs and the 
variance of perturbations. The set of observational constraints on the abundance of PBHs is 
listed in Table I and has be en used to place a bound to the mean amp li tude 5 in a variety o f 
cosmological models ( e.g. 



Carr et al. 



Sendouda et al. 



119941 1: iGreen & Liddld 1 119971 1: Idancy et all 1200311 : 
(2006J]). The Press-Schechter formula has also been tested against oth er 
estimations of the probability of PBH formation, such as peaks theory | Green et al. . 2004|] . 



The threshold value <5 th us ed in Eq. (14.231) has been modified with the improvement of 



gravi t ational collapse studies |Carr|, 



1975 



Niemeyer & JedamzikL 



1998 



Shibata & Sasak 



i, 



19991 : lHawke & Stewart! 1200211 . A more appropriate approach has be en noted recently, where 
simulations have addressed the problem using curvature fluctuations flShibata & Sasaki Il999l : 



Musco et al 



2005 



: |Polnarev & Muscol 1200711 . The co rresponding threshold value of the cur 



vature perturbation can be deduced from the relation [|Liddle & Lythi 12000(1 



S k (t) 



2(1 



w] 



5 + 3w 



k 



(4.25) 



which at horizon-crossing during the radiation-dominated era gives, TZtb = 1.01 for 5 



0.3. 



his value has 



3een also confirmed in the numerical simulations of 



tb 



Shibata & Sasaki 



lll999h . lGreen et all [1200411 and Musco et all B2005I1 . We will make use of it throughout. 

The threshold value 1Z th indicates the minimum amplitude of an inhomogeneity required 
to form a PBH. Consequently, the probability of PBH formation is best described by a non- 
linear treatment and this is the major motivation for our analysis. In the following we adapt 
the Press-Schechter formalism to derive the non-Gaussian abundance of PBHs. The use 
of the Press-Schech ter integral for distributions of the curvature perturbation is not new. 



Zaballa et al. 



ll2007ll use it to estimate the PBH formation from the curvature perturbations 
which never exit the cosmological horizon. We apply the integral formula in Eq. (14.231) 
to the non-Gaussian probability distribution (14.121) . The result of the integral is the sum of 
incomplete Gamma functions r inc and an exponential: 
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(4.26) 



The Taylor series expansion of these functions around the limit S^/^o = gives 
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For the mass fraction shown in Eq. (14.271) . the observational limits of Table I could in 
principle constrain the values of the variance Y? n and of /nl- However, when the mean 
amplitude of perturbations, E^, is normalised to the value at CMB scales, the obtained 
limits for /nl are of order 10 4 . This is inconsistent with the analysis presented here because 
in such regime higher order contributions are expected to dominate non-Gaussianity. In fact, 
the expansion in Eq. (14.91) shows that when 



\h 



NL 



< 



5 1 



3TZ 



RMS 



3E 



(4.28) 



the quadratic term of Eq. (14.91) dominates over the linear term, and in the computation of the 
three-point function Eq. (14.101) . the loop contributions to the correlators become dominant. 
For the values of 7£rms required to form a significant number of PBHs, the limit on |/nl| 
is of order 10. The computation of non-Gaussianities in this case goes beyond the scope of 
the present work . ( For discussions on the loop corrections to t he correlation functions, see 



Weinbergl H2QQ5I1 . 



Zaballa et al 



H2006H . Byrnes etalJ H2007H and lSeervi I2008IU 
It is interesting to look at the values allowed for /nl from WMAP and test the modifi- 
cations that large non-Gaussianities bring to the amplitude of 1Z at the PBH scale. Fig. 14.21 
presents the set of bounds on the initial mass fraction of PBHs listed in Table I. The corre- 
sponding limits to the variance of the curvature E^ are shown in Fig. 14.31 for the Gaussian 
and non-Gaussian cases. Independently of the model of cosmological perturbations adopted, 
one can use the observational limits on / NL to modify the bounds for on small wave- 
lengths. For the non-Gaussian case we choose to plot the central value of the present limits 
to f^ lL = 51 komatsu et all l2008ll and the limit value / NL = 5/(3E^) » -66 men- 
tioned in Eq. (14.281) . The tightest constraints on E^ come from perturbations of initial mass 
M « 10 15 g. With the non-Gaussian modification the limit is log (E^) < —1.2, compared to 
the Gaussian case log (E^) < —1.15. As shown in Fig. 14.31 the modification to E^ cannot 
be much larger if instead the limit value of Eq. (14.281) is used. 



4.5 Closing remarks 

The present chapter shows, to lowest order in the contribution of the bispectrum, the ef- 
fects of non-Gaussian perturbations on PBHs formation. Using curvature perturbations with 
a non-vanishing three-point correlation, an explicit form of the non-Gaussian PDF is pre- 
sented, which features a direct contribution from the non-Gaussian parameter / NL . Further- 
more, it is shown how the sign of this parameter determines the enhancement or suppression 
of probability for large- amplitude perturbations. Using the simple slow-roll expression for 
/nl in the context of single field inflation, a previous discrepancy in the literature regarding 
effects of non-Gaussianity on the abundance of PBHs has been solved. 

As a second application of the non-Gaussian PDF presented here is to use the Press- 
Schechter formalism of structure formation to determine the non-Gaussian effects on PBH 
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Figure 4.2: The constraints on /3pbh in Table I are plotted together with the smallest value considered for each 
mass. 
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Figure 4.3 : A subset of the constraints on from overproduction of PBHs is plotted for a Gaussian and non- 
Gaussian correspondence between (3 and E^, Eqs. (I4.24l i and ( 14.271 ) respectively. The dashed line assumes a 
constant /ml = 51 and the dotted line a value /nl = — V^tc ~ — 66- The solid line represents the constraints 
for in the Gaussian case 



abundance. In section 14.41 it is shown how the PBH constraints on the amplitude of per- 
turbations can be modified when a non-Gaussian distribution is considered. As an exam- 
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pie, it is shown that the limit S^(M = 10 15 g) < 6.3 x 10 2 is reached for t he marginal 
value /nl = —66.35 modifying the known bounds for 1Z at the end of inflation | Carr et all 



199411 . This limit is, however, much larger than the observed amplitude at CMB scales, 
where E-r. « 4.8 x 10~ 5 . The order of magnitude gap between the mean amplitude observed 
in cosmological scales and that required for significant PBH formation remains almost intact 
and, as a consequence, non-Gaussian perturbations do not modify significantly the standard 
picture of formation of PBHs. 



Chapter 5 

Curvature profiles of large overdensities 



5.1 Introduction 

As mentioned in Chapter [H previous studies of PBH formation take the amplitude of the mat- 
ter density or curvature inhomogeneities as the only parameter determining the probability 
density of PBH formation. Also the mass fraction in the form of PBHs is usually calculated 
with the aid of the Press-Schechter formula. Here we argue that this rough estimation is 
incomplete and that a different approach should be taken to evaluate the threshold value S th , 
or the equivalent curvature inhomogeneity lZ th , in the investigations of PBH formation. 

From the first numerical simulations of PBH formation, it was evident that the process 
of PBH formation depends on the pressu re gradients in the collapsing configuration as well 



as their amplitude. iNadezhin et al.l ll 197811 found that such pressure gradients can modify the 
value of 5 t h significantly. This has been confirmed in more recent works, which describe 
the configuration in terms of the curvature inhomogeneity 1Z(r) (note that in this chapter 
we work with spherical coordinates {r} and not any set of coordinates {x}). As we will 
show below, the Einstein equations relate the curvature profiles directly with the internal 
pressure gradients. This is the main motivation for considering the probability of curvature 
configurations. 

We extend here the Press-Schechter formalism to consider a two-parameter probability. 
We include here for the first time a parameter related to the slope of curvature profile at the 
edge of the configuration. We start by calculating the probability of finding a spherically 
symmetric curvature configuration with a given r a dial p rofile. We can justify the sphericity 



assumption using the argument of IZabotin et alJ tl987|]: PBH formation takes place only 



from nearly spherical configurations. In our analysis, we describe the radial profiles by 
introducing two parameters: the central amplitude of the curvature inhomogeneity lZ(r = 0) 
and the central second radial derivative 71" (r = 0). The introduction of these parameters 
is a first step towards the full parametrisation of profiles in terms of all even derivatives at 
the centre of configurations. (The odd derivatives are all zero due to the assumed spherical 
symmetry.) 

The method presented to derive a multiple-parameter probability enables us to compute 
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the probability of any number of parameters describing the curvature profile. However, only 
families of curvature profiles described by two parameters are currently available, so we limit 
ourselves to a two-parametric description. More accurate future codes will simulate PBHs 
formation with a larger number of parameters. The number of parameters required for the 
complete description of these profiles and their probability distribution will be the same0. 

The central amplitude 71(0) has been used in previous calculations of gravitational col- 
lapse and the probability of PBH formation, as illustrated in the previous chapter. Here we 
compute the probability to find a given configuration as a function of the two parameters 
[JZ(0),TZ"(0)]. We subsequently illustrate how this two-parametric probability is used to 
correct the probability of PBH formati on. For this purpose we us e the results of the latest 
numerical simulation of PBH formation IPolnarev & Mused Il2007n . Such an exercise shows 



how the corrections to /?pbh are potentially significant and they will be considered in more 
detail in future studies of PBH formation. 



5.2 Probability of profile parameters of cosmological perturbations 



Formally, the high amplitude inhomogeneous profiles describing configurations which col- 
lapse into PBHs are not perturbations. However, such regions are included in the statistics of 
random primordial curvature perturbations in the sense that the statistics of random fields can 
be used to estimate the probability of finding high- amplitude inhomogeneities. To describe 
su ch inhomogeneities, we consider the nonlinear curvature field lZ(t,r), as first described 
by ISalopek & Bond! 119901]. The nonlinear curvature 7Z(t, r), defined in terms of the met- 
ric in the following equation, represents the relative expansion of a given local patch of the 
universe with respect to its neighboring patches. It is described by the metric 



ds 2 = -N 2 (t, r) dt 2 + a 2 (t)e 2n ^ r) ^ tj (dr l + N l (t, r) dt)(dr J + N\t, r) dt), (5.1) 



where a(t) and 7 are the usual scale factor and the intrinsic metric of the spatial hypersur- 
faces. The gauge-dependent functions N and N l are the lapse function and shift vector, 
respectively. These variables are determined by algebraic constraint equations in terms of 
the matter density p, pressure p and metric variables 71, a and 7^. 

Here we consider the nonlinear configurations which correspond to large 7Z inside some 
restricted volume and zero 7Z outside, where the expansion of the universe follows the back- 



ground FRW solution. There are several advantages of working with metric (15.11) . First, as 
shown in Chapter El 7Z is defi ned as a gauge-invariant combination of the metric and matter 
variables ll Wands et all 1200011 . Second, with the aid of the gradient expansion of the metric 
quant i ties, 1Z(r, t) appe a rs in t h e Einstein equations in a no n -perturbative wav lIStarobinskyL 



19861 : ISalopek & Bond . 



1990 



Deruelle & Langlois , 



1995 



Rigopoulos & Shellard . 1200511 . 



In the context of dark matter haloes, the profile of the initial inhomogeneity is effectively irrelevant because galaxies 
are formed from pressureless configurations. The density profiles and shapes of virialised haloes result from the evolution 
of the initial high peaks and are not linked to the profile of initial configurations that we investigate here 
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Th ird, 7Z does not d epen d on time for scales larger tha n the cosmological horizon, as proved 



by lLyth et al. 



and 



Langlois & Vernizzil 12005 afl . In the present chapter we work in the 



superhorizon regime, where the field 7Z(t) can be assumed to be time-independent. 

The primordial field of random perturbations we consider presents a Gaussian probability 
distribution. The expressions for the PDF of the parameter 71(0) in Chapter [3] are recovered 
here. For convenience we use a different notation, replacing g in Section 13.3.11 with the 
amplitude i3 = 71(0) and the variance with S( 2 ) = The PDF for the central amplitude 
$ is identical to that in Eq. (l3.53l) : 



1 



2tt£ 



exp 
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(5.2) 



We now derive the density of the probability for the central second derivative to have 
amplitude 



tf 2 = 71" (0) 



d 2 



(5.3) 



J r=0 



In order to compute the probability of a specific property of 7Z(r), we integrate the original 
PDF, which encodes all the information about the field, weighted with the Dirac 5-functions 
of relevant arguments. Hereafter we assume 71(0) and 71" (0) as statistically independent 
parameters. The validity of this assumption is not explored here but is left for future inves- 
tigations. Following this assumption, the probability of having 71" (0) = $ 2 is given by the 
integral 



[dK]F(7l) 8 [K"(0) - t? 2 ] 



(5.4) 



where \d7V\ indicates integration over all possible configurations 7t(k) in Fourier space. In 
order to compute this integral, we expand the smoothed curvature perturbation profile 7Z(r) 
in terms of spherical harmonic functions: 



n(r) 



d 3 k 



7Z{k) exp (ik ■ r), 



(5.5) 



with 



^( k )=E E J2 n ^ m (0A)Mk)- 



(5.6) 



e=0 m=-l n=l 



Here Y^ m are the usual spherical harmonics on the unit 2-sphere and ip n (k) are a complete 
and orthogonal set of functions in an arbitrary finite interval < k < A. (The explicit 
expression for ip(k) and the value of A are given by Eq. (13.371) of Chapter[3l) The coefficients 
in the expansion are genetically complex, so we separate the real and imaginary parts by 
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introducing TZ™ n = a^ n + ib^ n . The reality condition for the curvature field, 1Z* (k) 
TZ(— k), is met when 



tin 



\l+m+l ira 



(5.7) 
(5.8) 



In particular, the m = modes require a°i n and b^ n to be zero for odd and even £, respec- 
tively. To integrate (15.41) we use the Fourier expansion (15.51) so that 



d 3 k 

W) 1 



1Z{k){\k) 2 exp (ik-r)| r=0 . 



K"(0) 



Furthermore, we use (15.61) and the orthogonality of the spherical harmonics 



(5.9) 
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(5.10) 



to obtain 

oo 
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5° 
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i; n (k)k A = tf 2 . 



(5.11) 



To compute the probability (15.41) we proceed by integrating over all configurations in 
Fourier space. With the aid of the expansion (15.61) we can express the measure of the integral 
in terms of the expansion coefficients satisfying the reality conditions (15.71) and (15.81) . as 
shown in Eq. (13.341) this is 
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p=0g=l 



(5.12) 



for any given functional ^ of 7£(k). The constants ji and /2 are weight factors to be included 
in the final normalisation of the joint probability. 

The Gaussian PDF we are restricted to is written in terms of the spherical harmonic 
coefficients as (cf. Eq. (13.491) ) 



oo £ oo 



m = exp [ - W + W 2 ] 



e=0 m=0 n=l 



oo oo 

EE 
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(5.13) 



^7° I 2 
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l^2p+l| g | 2 ] 



p=Q q=l 

In order to obtain the probability in Eq. (15.41) . we use the standard representation of the Dirac 
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5-function 

/oo 
dz exp[izx]. 
-oo 



(5.14) 



This allows us to write the (^-function in Eq. (15.41) in terms of the spherical harmonic coeffi- 
cients as 



6[K"(0)-& 2 ] = J dzexp 



ME 



'On 



"2\n 



,n=l 



dk 



(5.15) 



We now have all the elements required to integrate probability density of finding 72." (0) with 
amplitude i? 2 - Substituting expressions (15.131) and (15.151) into Eq. (15.41) . we perform the 
functional integral with the aid of the decomposition (15.121) . For this case we have 



P(tf 2 ) oc / [dK] / dzP[72]exp 
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3(2vr) 3 tf 2 



47T 
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where we have simplified the expression by defining the factor 



4 4 > = / dkk 4 Mk). 



(5.16) 



(5.17) 



In the process of integration, we discard all the Gaussian integrals because they contribute to 
the probability only with a multiplicative constant which will be included in the final normal- 
isation. On the other hand, the Dirac 5-function contributes to the integral with exponential 
functions of ajL and a% n . The integrals of these parameters are computed by completing 
squares of the exponential arguments. First, we collect the terms of the integral with factors 
of a° |n , that is, 



exp 



^ oo oo 



47T 2 (27T 

Completing the squares, this last expression becomes 
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In the same way we can complete the squares for the expansion factors a% n : 
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(5.20) 
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Finally we can complete the squares for the terms containing the variable z, these being 
independent of a^ n and a°| n : 
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where for simplification we have written 
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(5.22) 
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So by making the change of variables 
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we can perform all the integrals and eliminate the Gaussian ones which contribute only up to 
an overall numerical factor subsequently absorbed by normalisation. The remaining factor 
expresses the probability of finding a perturbation 1Z with a central second derivative of value 

tf 2 : 



F[K"(t = 0) = 2 ] oc exp 



2S (4) 



(5.23) 



The quantity E 2 4 ^ represents the 'variance' of the PDF for 72." (0). To evaluate this variance 
we integrate Eq. (15.221) and use the property (13.381) in Chapter [3] to integrate the complete 
sum and obtain 



Y 2 

^(4) 



dhikW 2 (k,k H )V(k) k 4 



(5.24) 
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The final probability density for the pair of parameters 72.(0) and 72" (0) is the product of 
Eqs. (Q and (15331) 



P ( 72(0) = O , 72"(0) = 2 ) = A exp ( 



r) 2 
?Y 2 

^(2) 



/2j (4) 



(5.25) 



Here £( 2 ) and E( 4 ) are the dispersion of the amplitude and the second derivative respectively, 
and A is a normalisation factor obtained from the condition that the integral of the joint PDF 
over all possible values of the two independent parameters equals unity. The final normalised 
joint probability density is 



P(0o,0 2 
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It is worth mentioning that the standard PDF containing only amplitudes O , Eq. (15.21) . is 
recovered from Eq. (15.251) when we set all gradients in the Hubble scale equal to zero, i.e. 
V72| r=rH = 0. The Fourier transform of this expression demands \k H \ — »■ oo. Using this in 
Eq. (15.241) means that Ea) — > oo and the argument 2 goes to zero in the probability density 
of Eqs. and OS- 



According to the Press-Schechter formalism llPress & Schechterl . 1197411 . the PDF is inte- 
grated over all perturbations which collapse to form the astrophysical objects under con- 
sideration. In this way we calculate the mass fraction of the universe in the form of such 
objects. To apply this formalism and calculate the probability of PBH formation and inte- 
grate the PDF (15.261) . we require the range of values 72.(0) and 72." (0) which correspond to 
PBH formation. In the next section we will obtain this range w ith the help of the results of 
numerical computations presented by lPolnarev & Mused Il2007ll . 



5.3 The link between perturbation parameters and the curvature pro- 
files used in numerical calculations 

5.3.1 Initial conditions 



As demonstrated by the first numerical simulations of PBH formation IINadezhin et al J . 1 1 97 811 . 
whether or not an initial configuration with given curvature profile leads to PBH formation 
predominantly depends on two factors: 

• The ratio of the size of the initial configuration r to the size of the extrapolated closed 
universe r k = a(t) dr/y/1 — r 2 , which is a measure of the strength of gravitational field 
within the configuration. 

• The smoothness of the transition from the region of high curvature to the spatially flat 
FRW universe, which is characterised by the width of the transition region at the edge of 
the initial configuration and is inversely proportional to the pressure gradients there, strong 
pressure gradients inhibiting PBH formation (This is an argument beyond the Jeans' stability 
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criterion and applies to configurations beyon d the linear regime) 



The numerical computations presented in lPolnarev & Mused 11200711 (hereafter PM) give 
the time evolution of the configurations with initial curvature profiles accounting for the 
above-mentioned factors. In that paper the initial conditions are obtained with the help of a 



quasi-homogeneous asymptotic solut ion valid in the limit t — » 



stein e quat ions was first introduced b y llifshitz & K halatniko villi 96311: s ee also 



Nadezhin et al. 



his solution to the Ein- 



Zeldovich & Novikov 



119781], PM used this asymp- 



119831] and lLandau & Lifshitzllll975n . Following!] 
totic solution to set self-consistent initial conditions for curvature inhomogeneities, the ini- 
tial curvature inhomogeneity being described by the spherically symmetric curvature profile 
K(f). This sets the initial conditions for the process of black hole formation. Asymptotically, 
the metric can be presented in terms of K(f) as 



ds 2 



-dr) 2 + s 2 (fi) 



1 - K(f)f 2 



dr 2 + f 2 [d9 2 + sin 2 6 1 



(5.27) 



where 77 is the conformal time, s(i]) is the scale factor for this metric. As we will show in the 
paragraph after Eq. (15.421) . this is identical to the usual scale factor 0(77) of a fiat Friedmann 
universe, only here we use a different notation to distinguish between metrics. Also, we 
write f for the radial coordinate to distinguish it from the coordinate of the metric (15.11) . 

An advantage of working with this metric is that it contains the curvature profile K(f) 
explicitly. We choose a set of coordinates with the origin at the centre of spherical symmetry 
and fix K(0) = 1. The condition that K(r) is a local inhomogeneity requires that K(f) = for 
radii f larger than the scale f where the metric matches the homogeneous FRW background. 

In PM the profiles K(r) are presented in two forms, one of which is characterised by two 
independent parameters a and A as 



K(r) 



1 + a 



2A 2 



exp 



2A 2 



(5.28) 



The parameter A describes the width of the Gaussian profile, while a parametrises linear 
deviations from this profile. The results of the numerical simulations in PM indicate that 
PBHs are formed in the region of the parameter space [a, A] shown in Fig. 15.1a . 



5.3.2 Physical criteria for the identification of parameters 

We proceed to find the correspondence between the two sets of parameters, [TZ(0),TZ"(0)] 
and [a, A], both of which describe the initial curvature profiles. First let us note that the sets 
of coordinates {t, r} and {77, f } are those of the metrics (15.11) and (15.271) . respectively. Thus 
we require a relationship between these set of coordinates too. Assuming that the size of the 
configuration, r , is much larger than the Hubble radius, r H = H" 1 , we can use the gradient 
expansion of the functions in metrics (15.11) and (15.271) . In this case, the time derivative of any 
function f{t, r) is of order f/t ~ Hf and significantly exceeds the spatial gradient which is 
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Figure 5.1: (a) The top plot shows the parameter values for initial configurations which collapse to form black 
holes according to iPolnarev & Mused 1200711 . (b) In the [R.(0), 1Z" (0)] plane three regions of integration are 
considered to compute the probability of PBH formation. Area I is the region enclosed by the solid curves 
and corresponds to the area denoted by BH in Fig. la. Area II is the region to the right of the grey dotted 
fine, representing the area of integration considered in previous studies where only the amplitude is taken 
into account. Area III is the region above the solid line and between the dashed lines. This contains those 
configurations which have a smooth profile in the centre and present the amplitudes 71(0) that are found to 
form PBHs in llPolnarev & Muscol l2007Tl . The complete description of the physical characteristics of profiles 
with values in this region is given in Section l531 



of order f /r . Hence the small parameter in the gradient expansion is 

6 = = — , (5.29) 

r ah 

where k is the wave-number corresponding to the scale of the configuration. 

For the metric (|5.1I) . using the coordinate freedom to set N l = and ignoring any tensor 
contributions (i.e. 7^ = Sij), the expansion of the Einstein equation G ° = 87rGT ° to order 
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e 2 can be written as0 

1 /6c 

2 W 2 



\(% + (3) R - ~ir(N - 1) ) f 0(e 4 ) = 8nG (p + Sp) + 



(5.30) 



where ^ R is the spatial curvature, or the Ricci scalar for the spatial metric . To zero order 
in e, we have 



— = 8vrGp , 



(5.31) 



which corresponds to the homogeneous part of Eq. (15.301) . The time-slicing can be taken to 
be the uniform expansion gauge in which 



N 



l + 3w 
1 + w 



S p + 0(e 4 



(5.32) 



where w is the equation of state UShibata & Asadal . 



1995 



Shibata & Sasaki 



1999l : lTanaka & Sasaki 



20070 ■ Using (I5.30I) . (I5.31I) and (15.321) . we find the equivalence between the spatial curvature 
and the matter overdensity: 



7 + 3w 



3 r V3 + 3w / 

In consequence, the gradients if this quantity relates to the pressure gradient: 

7 + 3w 



(5.33) 



V (3) R = — - - V (Sp) = — - 

6 6{W + 1) 6 



3w (w 



(5.34) 



where V = (g rT ) _1 / 2 d/ 'dr. Hence, subject to the two physical conditions at the edge of the 
configuration listed at the beginning of Section [5731 we relate the profiles lZ(r) and K(r) by 
equating the spatial curvature and its gradient for the metrics (15.11) and (15.271) . That is, 



(3) R 



27e"(r) + (n'(r)f exp(-2TL(r)) = 3K(f) + rK'(f) 



(5.35) 



and 



1 d 



(®R) 



^/g^dr 
[7e , 7e ,, + 7e w ]exp(-37e(r)) 



"1 - Kf 2 " 


1/2 ^ 







2fK'(f) + ^f 2 K"(f; 



(5.36) 



By definition, the 3-curvature must vanish at the edge of the configuration, so Eq. (15.351) 



11 Fo r the complete second order expansion of the metric quantities, see iLvth et all l2005ll and lLanglois & Vernizzil 
l2005bll . 
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implies 

2n"(r ) + (n'(r )) 2 = (5.37) 

and 

3K(f ) + foK'(fo) = 0. (5.38) 
Thus the gradient relation (15.361) can be written as 



[JZ'(r ) 3 - 2K'"(r )] exp(-3^(r )) 



1 — Kri ' 



12K(f )+f 2 K"(f )]. (5.39) 



This establishes a relation between 1Z(r) and K(f) at the edge points r and f . The config- 
uration K(f) is parametrised by [a, A], as shown in Eq. (15.281) . As follows from condition 
(15.381) . the radius r can be written in terms of those parameters as 



5a - 2 + y/(5a - 2) 2 - 24a 



2a 



A 2 . (5.40) 



Then we use two more equations obtained from the conformal transformation of coordinates 
at zero order in e: 

a 2 (r)e 2 ^ dr 2 = s 2 ( V ) g— (5.41) 

and 

a 2 (r) e 27 ^ r 2 drf = s 2 ( V ) f 2 dtf. (5.42) 

Asymptotically, in the limit [r, f] — * oo, the homogeneous Einstein equations are identical in 
both metrics, therefore, the homogeneous scale factors a(r) and s(r)) can be identified. Thus 
we find a relation between the radial coordinates, 

e n(r) r = ~ ? (5 43) 

and an integral relation between the configurations, 

e n ^dx= f , dX (5.44) 
o Jo y/1 - K(x)x 2 

One can verify that Eqs. (15.351) . (15.361) and (15.441) are not independent. For example, Eq. 
(15361) follows from (15351) and (f5~44b . 

In the previous section we have developed a method to account for the probability of any 
set of parameters describing the curvature profile. For simplicity we have chosen the pair 
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[11(0), K"(Q)]. We now illustrate how to relate [71(0), K"(0)] and [a, A] by considering the 
parabolic profile 

TZ(r) = 72(0) + -K"(0) r 2 . (5.45) 

This parametrisation meets the minimal requirement of covering the [a, A] parameter space 
in Fig. 0(a). 

Eqs. (15.371) . (15.441) and (15.431) are now reduced to the following system of algebraic equa- 
tions: 

2 



„2 _ 

n"(oy 



< = -^TTTTTV, (5-46) 



2 f f '° dx \ 

= _ 2 exp(2R(0) - 2) 

where f is given in terms of [a, A] by Eq. (15.401) . 
5.3.3 Parameter values leading to PBH formation 

The numerical computations of PM, which used the parametrisation (15.281) . show that PBHs 
are formed in the [a, A] region shown in Fig. l5.lf a). Eqs. (15.471) and (15.481) map this region 
to Area I in the [71(0), 72." (0)] plane shown in Fig. 15. lb . The Jacobian of the transformation 
corresponding to this mapping is non-vanishing, which guarantees a one-to-one correspon- 
dence of the 'BFT region in Fig. 15.1a with Area I in Fig. 15.1b . Each point here corresponds 
to a parabolic profile which leads to PBH formation. 

For each one of these parabolic profiles, there is a family of non-parabolic profiles with 
the same central amplitude 71(0), the same configuration size r , and the same behaviour 
near the edge, as shown in Fig. 15.21 In this figure, the profiles lying below the parabola 
correspond to larger absolute magnitudes of 72" (0) and do not form PBHs because they have 
lower average gravitational field strength and higher average pressure gradient. The non- 
parabolic profiles which lie above the parabolic one (with smaller absolute magnitude 72" (0)) 
should also collapse to form PBHs because they correspond to higher average gravitational 
field strength and lower pressure gradient. 

In the parameter space [72(0), 7Z"(0)}, this last set of profiles corresponds to Area III in 
Fig. I5. lb . This region will be included in the calculation of the probability of PBH formation 
in the next section. 
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Figure 5.2: The curvature profile for three different families of configurations with common central amplitude 
TZ(0) = 1. The configurations shown by the dashed lines have values of 7Z"(0) larger in absolute magnitude 
than the parabolic one shown in black. The configurations shown by the dotted lines have values of 1Z"(0) 
smaller than the parabolic one. All profiles satisfy conditions (15.37b and ( 15.391 l. 



5.4 Two-parametric probability of PBH formation 



To calculate the probability of PBH formation, which is equivalent to the mass fraction 
of the universe going into PBHs of given mas s, it is customary to use the standard Press- 
Schechter formalism [jPress & Schechten . ll974|] . This has bee n widely used in previous cal 



culations of the one parametric probability of PBH formation | Carrl 



Liddle & GreenL 1 19981: ICan. 



1975 



Carr et al 



Zaballa et al.1 



1994 



2007]. 



200.4 Ichongchitnan & Efstathioui 12007 
When the probability depends on a single amplitude parameter, this method reduces to the 
integration of the corresponding PDF over the relevant pe rturbation a mplitudes. The final 



integral is equivalent to the mass fraction of PBHs of mass llCarrl 1197511 



M ~ w 3/2 M H « w 3/2 k M /(2n) 



(5.49) 



with the equation of state w measured at their formation time. Here we extend the standard 
Press-Schechter formalism to include for the first time an additional parameter accounting 
for the radial pressure in the initial configuration. When the [TZ'(0), 1Z"(0)} area is a square 
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[Hi < 11(0) < 1Z 2 , K" < 11" (0) < fc'2), the integrated two-parametric probability is 



TL 2 



(3 PBH (M) =2 / d& Q / d0 2 P(0o,0 2 ) 




(5.50) 



We use this result to integrate numerically over a mesh of small squares covering each of the 
areas of the plane [7^(0), 71" (0)] shown in Fig. 15.1b . The results of this integration for two 
different power-law spectra Vn{k) oc /c n_1 are shown in Fig. 15.31 

From that figure we note that the probability function Ppbh(M) has a maximum at a 
value of M max that changes with the of the spectral index. This can be easily derived by 
computing the solution of dppBn/dM = 0. We find that this equation provides a formula 
for the value M max which indeed depends sensitively on the spectral index n s . Assuming 
n, > 1 we have: 



M max = M eq 



3 7 P< 



eq 



2 V 2 



cxp 



n„ - 1 



(5.51) 



where M eq and P cq are the Hubble mass and the power spectrum at the time of matter- 
radiation equivalence (A; cq = 8.9 x 10 _2 Mpc~ 1 ), and 7 is a factor of order unity that changes 
slightly with the value of n s . 

We contrast the case of parabolic profiles described by Eq. (15.451) with the non-parabolic 
set presented in Fig. I5.2l by plotting the probability /?pbh for different values of V-r,- This is 
presented in Fig. 15.41 The figure shows that the probability of PBH formation can be larger 
th an the one-paramete r probability computed in previous studies from the integration of Area 
II ||Green et all 12004(1 . This important result requires confirmation from more detailed nu- 
merical simulations of PBH formation in this parameter area. The uncertainty is explained 
by the fact that the two-parametric calculation of the probability of PBH formation is still 
incomplete. This should be complemented in the future by the introduction of all relevant 
higher-order derivative parameters and the higher-order correlations in the PDF. 



5.5 Discussion 

We have developed a method for calculating the two-parametric probability of PBH forma- 
tion, taking into account the radial profiles of nonlinear curvature cosmological inhomo- 
geneities. This is the first step towards calculating the iV-parametric probability, which takes 
into account the radial profiles more precisely than studies using the amplitude as the only 
relevant parameter. We have incorporated the derived contribution to the total probability of 
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Figure 5.3: The logarithmic probability of PBHs for two tilts in the power spectrum (n s = 1.23 on the top 
figure, n s = 1.47 on the bottom figure), integrated for the three different regions sketched in Fig. 15.11 The 
integrals over Areas I and II correspond to the dashed and solid lines, respectively. The probability integrated 
over Area III is represented by the dotted lines in both figures. 



PBH formation by considering the range of values of 1Z"(ti) that will form PB Hs, using the 
results of the numerical computations presented by IPolnarev & Mused Q2007[| . Finally, we 
have provided an example of the consequences of this probability for the statistics of PBHs. 

The results obtained show that, if we restrict ourselves to PBH formation from parabolic 
profiles (as described in Section 15.31) . then the total PBH probability is orders of magni- 
tude below previous estimates! On the other hand, if non-parabolic configurations are also 
included (see Fig. 15.21) . the total probability of PBH formation is higher than the single- 
parametric probability estimated in previous works. In this case, we can impose new bounds 
on the power spectrum in the scales relevant for PBH formation. Analysing the uncertainty of 
our results, we have demonstrated how much we still have to understand about the formation 
and statistics of PBHs. The physical arguments supporting our results should be verified by 
numerical hydrodynamical simulations of PBH formation, which would provide a valuable 
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Figure 5.4: The grey dashed line shows the ratio of the total probability /?pbh which results from integrating 
over Area I on the [72.(0), 1Z" (0)] parameter space of Fig. lb to the probability which results from the integrating 
over Area II. The black line is the ratio of the probability integrated over Area III to the probability integrated 
over Area II. 



feedback to the initial motivation of this work. 

The main argument of this chapter is that the amplitude of initial inhomogeneities is 
not the only parameter which determines the probability of PBH formation. The ultimate 
solution of the problem requires a greater set of parameters and a larger range of their values 
to determine all high curvature configurations that form PBHs. This is a huge task for future 
research. In the meantime, we have a method to operate with the statistics of all these 
parameters. 



Chapter 6 

Conclusions and future work 



In this thesis we have presented a study of large inhomogeneities in the early universe. Such 
large concentrations of matter may collapse to form primordial black holes (PBHs). The 
number of PBHs in our universe is calculated by integrating the probability distribution func- 
tion (PDF) of primordial inhomogeneities, this encoding all the statistical information of pri- 
mordial inhomogeneities. The main objective of this thesis is to quantify the probability of 
PBH formation in the context of nonlinear perturbation theory. This represents a significant 
improvement in the study of large- amplitude inhomogeneities since, by definition, these are 
nonlinear. 

The statistics of inhomogeneities are the point of contact between theory and observa- 
tions. In theoretical studies the statistics of primordial fluctuations are studied in the frame- 
work of cosmological perturbation theory. Until recently, perturbation theory was restricted 
to consider only linear departures from the homogeneous background. Linear perturbations 
are Gaussian due to the independence the perturbation modes. This is an excellent approxi- 
mation to describe the structures observed in the universe. Indeed, observationally, only the 
variance, or second statistical moment, has been measured. However, the detailed observa- 
tions of large-scale structure (LSS) and the cosmic microwave background (CMB) now allow 
us to test for corrections to the linear approximation. This motivates the study of extensions 
of linear perturbation theory. In particular, the non-Gaussianity of curvature fluctuations has 
been a subject of intense investigation. 

6.1 Summary of results 

In Chapter [2] we presented a brief introduction to cosmological perturbation theory within 
general relativity. We reviewed the basic results of this theory for the cosmological inflation 
paradigm. From the evolution equations, we identified the conditions under which curvature 
fluctuations can grow significantly at superhorizon scales. As shown in Eqs. (12.871) and (12.881) . 
these conditions are mainly the presence of a non-adiabatic component in the matter field 
fluctuations. This has motivated several previous studies of non-Gaussianity resulting from 
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the second-order perturbations of an isocurvature (non-adiabatic) field x- The conditions for 
inflation show that a non-adiabatic field is only a subdominant component of the total matter 
during inflation. The cosmological model in which x iS responsible for the curvature pertur- 
bations is called the curvaton model. 



In Section [231 we calculated non-Gaussian correlators for some special cases of the cur- 
vaton model. The lowest order non-Gaussian signature is a non-vanishing skewness or third 
moment of the PDF. In perturbation theory this is equivalent to the correlation of three copies 
of the curvature perturbation field. This correlator itself is only present when we consider 
nonlinear perturbations. In order to derive the three-point correlator we have considered non- 
linear field fluctuations Sx- We have calculated the second order perturbations of a single 
isocurvature field during inflation and radiation domination. We have done this by solving 
the Klein-Gordon equation of the perturbation 5x to second order. For simplicity we consider 
only the matter fluctuations, assuming a large contribution from the third derivative of the 
potential d 3 W(x)/dx 3 - We find that an effectively massless field does not generate a large 
nonlinear contribution to the perturbation 5x- Conversely, a slightly massive field allows an 
exponential growth of the nonlinear perturbation. With the aid of a new method to compute 
non-Gaussian correlators, we derived the field bispectrum -F(kj) given by Eq. (12.1691) . We 
then derived the curvature perturbation bispectrum -B(kj), considering a dominant contri- 
bution from the field bispectrum F(ki). Equation (12.1881) expresses the the non-Gaussian 
parameter / NL in terms of the elements of the potential W(x)- Chapter [2] closed with a brief 
discussion of the observational limits to the curvaton. 

One of the main objectives of this thesis is to present the modified probability of struc- 
ture formation from the non-Gaussian PDF. From the central limit theorem, we know that 
the non-Gaussian PDF produces non-trivial moments of order higher than two. To determine 
the shape of the distribution uniquely, one requires the a priori knowledge of all moments. 
From studies of non-Gaussianity in perturbation theory, however, we only know the skew- 
ness (third moment) and kurtosis (fourth moment) of some models of structure formation. 
Finding a PDF which encodes the contribution of only these two higher order moments is not 
trivial. In Chapter [3] we have constructed, in the context of quantum field theory, the general 
non-Gaussian PDF for curvature perturbations 1Z. Formally, this is a probability functional 
for the ensemble of realisations of 7£(x) at some specified time t. We refer to this probabil- 
ity as Pt[7£]. We first derived a mathematical expression of the above statement by writing 
f t [JZ] in terms of the n-point correlation functions (see Eqs. (13.121) and (13.151) ). We then 
constructed an explicit expression for the PDF using only the first three statistical moments. 

We found that, in order to calculate the PDF, it is necessary to consider a regularised 
function R(k). We must therefore consider a field sufficiently smooth on small scales, with 
a smoothing scale customarily set as the horizon scale at time t. We also require an upper 
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limit for the A;-numbers in order to avoid divergences in the integrations required to construct 
the PDF. This is achieved by artificially compactifying the momentum space over a scale 
A > k. This regularisation is a common requirement of calculations in the theory of pertur- 
bations and in the statistics of LSS. The results of this chapter are important because the final 
expression is an explicit functional probability P t [7£(k)]. This means that the probability of 
any parameter appearing in the function 1Z(k), or equivalently 7£(x), can be retrieved from 
this PDF. We rely on this property to study two important modifications of the probability of 
PBH formation in the subsequent Chapter|4]and Chapter[5] 



In the last two chapters of this thesis we revisit the calculation of the probability of PBH 
formation, taking into account two important effects which are characteristic of nonlinear 
inhomogeneities. In Chapter |4] we calculate the probability of PBHs using a non-Gaussian 
PDF we considered the non-Gaussian PDF of curvature perturbations 7Z. The featured PDF 
includes a linear contribution from the three-point correlation function as derived in Chap- 
ter [3] In Section I4~2l this PDF was adapted to curvature configurations 7Z(r) that give rise 
to PBHs. As previous works show, the amplitude at the centre of the curvature configuration 
7Z(r = 0) is a good parameter to determine the formation of PBHs. Eq. (I4.12I) gives the non- 
Gaussian PDF for the mentioned parameter. With the aid of this PDF we have reproduced 
qualitatively the effects of the non-Gaussian contribution considered in two previous works. 
We first identified the source of inconsistencies in previous works studying non-Gaussian 
effects in the probability of PBH form ation. We showed that the fundame ntal dif f erence in 
the inflationary models considered by iBullock & Primackl jl997|l and by llvanovl ||1998|l is 
the spectral index. In the first work, perturbations involve a blue power spectrum (for which 
the spectral index accomplishes n s — 1 > 0), while the power spectrum is red (n s — 1 < 0) in 
the second. Noting that, in single-field inflation, the non-Gaussian parameter / NL is directly 
related to the spectral index n s — 1, we have shown the source of the discrepancy. The main 
effect on the non-Gaussian PDF is the respective suppression and enhancement of the prob- 
ability for large values of 71. 



Chapter |4] also presented the non-Gaussian modifications to the probability of PBH for- 
mation (or the mass fraction of PBHs). In Section l4~4l we have shown how the new PDF 
can modify the bounds to the variance of curvature fluctuations This comes from the 
observational limits to the abundance of PBH for each mass scale. Such modifications are 
illustrated in Fig. 15.41 Note that in this figure we have used the maximum value of the pa- 
rameter / NL allowed by perturbation theory. In the future, greater values could be considered 
by constructing PDFs with the techniques described here. 

In Chapter [5] we have studied the probability of configurations 7£(x) from another per- 
spective: We compute the probability of a parameter describing the curvature profile in ad- 
dition to the probability of 71(0). Specifically, we compute the probability of the second 
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radial derivative at the centre of the configuration, 11" (0) = d 2 lZ/dr 2 \ r=0 . As studied in that 
chapter, the consideration of additional parameters describing curvature profiles is a signifi- 
cant improvement in the study of gravitational collapse. In other words, the choice of initial 
configurations collapsing to form PBHs relies on two sets of parameters. Parameters of the 
profile lZ(r) are required in addition to the amplitude parameters customarily used. We used 
the results of the latest simulations of PBH formation to integrate all the allowed configu- 
rations parametrised with the pair [1Z(0),1Z"(0)]. The result shows, heuristically, how the 
probability of PBH formation can be drastically changed by considering curvature profiles 
in the PDF. 



6.2 Future research 

The non-Gaussian signatures of cosmological inhomogeneities offer good prospects for model 
discrimination. An example of this is given at the end of Chapter [2], where we were able to 
limit special cases of the curvaton model with the observed constraints on the parameter /nl- 
A number of extensions to this work are possible. First, the curvaton m odel can be adapted , 
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with pertinent modifications, to describe models of modulated reheating |Zaldarriaga 
In such models, the auxiliary field modifies the expansion in different patches of the universe 
during the reheating process. The present work can be extended to cover such models by 
modifying the scales of the mass and expectation values of the auxiliary field \- I n this wa Y 
one can search for feasible models of modulated reheating which satisfy the observational 
limits of non-Gaussianity. 

Another application of our study is to compute higher-order correlations from the derived 
solutions to the nonlinear Klein-Gordon equation. Future probes of non-Gaussianity could 
detect the 'trispectrum' of curvature perturbations, which is a higher order discriminator be- 
tween models of inflation. Computing the corresponding four-point function is thus crucial 
for a characterisation of the hypothetical detection of non-Gaussianity at this level. 

An important complement of the work presented in Chapter [2] is the computation of the 
solutions to the Klein-Gordon equation allowing for metric perturbations. This has been ig- 
nored here because we assumed that the field fluctuations dominate over all other sources, 
as applies in the slow roll limit of the Klein-Gordon equation. As mentioned in Section [231 
however, the curvature perturbation contribution (often called backreaction) may entail im- 
portant corrections for b\ [see e.g. Malikl [2007]]. It is important to compute such contri- 



butions because, on the one hand, they could be the dominating component in the growth of 
the fluctuations and, on the other hand, the curvature back-reaction could cancel large /nl 
values. This could prompt reconsideration of models previously excluded by observations. 
Which case applies is an open question that should be addressed in the near future. 
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There is another set of problems where the methods of Chapter [3] find an important ap- 
plication. This is the determination stochastic sources in the evolution equations of classical 
fields. A functional probability is written in terms of products of the n-point correlation 
functions with n copies of the field configuration in Fourier space (see Eq. (l3. 151) for the case 
of curvature perturbations). This can be used, in particular, to derive and extend the stochas- 
tic equations of inflation by Istarobinsky & Yokoyamia Il994 1. This seminal work presents 
a Fokker-Planck equation for the probability of the configuration 0(x) on small scales and 
for a single-field inflationary field. A first c onnection betw een the stochastic framework and 
the work presented here has been given by [SeeryJ [|2OO90. In that paper, a wave-functional 
similar to Eq. (13.131) is considered and the governing Hamiltonian operator for the scalar field 
is recovered from its action. The Fokker-Planck equation suggested by Starobinsky can 
be deduced easily from the Schrodinger equation for that wave-functional. Such a method 
can be extended to calculate PDFs of multi-scalar or non-canonical models of inflation. The 
construction of the probability distribution for configurations 0(x) and other possible fields 
allows for the consideration of full non-Gaussian distributions. This is clearly the way to go 
beyond approximations like the one considered in Chapter [3] 



Regarding the probability of PBH formation, the results of Chapter [5] cannot be conclu- 
sive because we do not have at hand the complete set of collapsing configurations. Deter- 
mining the set of all configurations collapsing to form PBHs is a huge task to be explored 
elsewhere. We can assert, however, that if PBHs are to be used as a tool for cosmology, the 
curvature profile parameters have to be taken into account in the derivation of the PDF. A less 
ambitious task is to have an estimate of how severe the modifications to the single-parameter 
approximation are. This would require the determination of more appropriate parameters 
describing curvature profiles, a topic currently under investigation. 

From the results of Chapter |4] we are able to set constraints on inflationary models. 
Specifically one can look at models with an enhancement of the power spectrum at small 
scales. In such cases, the constraints from PBHs can be more or less stringent, depending on 
the values and the sign of the non-Gaussian parameter / NL . Here we have provided a tool for 
testing those models. Such tool can also be improved as more constraints are derived from 
observational tests to the abundance of PBHs. 
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